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ABSTRACT 


High  Resolution  Measurements 
O'f  OH  In-frared  AirglOM  Structure 

by 

Parris  Cornel  Neal •  Doctor  04  Philosophy 
Utah  State  University!  1985 

Major  Processors:  Or.  Doran  J.  Baker  and  Dr.  Kay  O.  Baker 
Department:  Electrical  Engineering 

Disturbances  in  the  normally  calm  atmospheric  airgloM 
layer,  which  cause  bright  and  dark  bands  or  '^str i pes*^  to 
appear,  have  been  observed.  These  disturbances  are 

attributed  to  '"^gravity  waves*  propagating  through  the 
atmosphere.  An  instrument  capable  oC  resolving  the 
temporal,  spatial,  and  spectral  attributes  o4  OH  infrared 
emissions  was  designed  to  gather  quantitative  data  on 
airglow  structure. 

An  optically/compensated  interferometer  /4pect  rometer  was 

chosen  as  the  basic  instrument  to  measure  this  phenomenon. 

2 

This  high-throughput  instrument  (0.285  cm  sr )  is  an  order 
of  magnitude  more  sensitive  than  more  conventional 
spectrometers  having  a  noise  equivalent  spectral  radiance  of 
16  R/cm  ^  at  1.5  Mm.  A  spectral  resolution  of  2  cm  ^  was 
obtained.  The  high-throughput  of  the  optically-compensated 


i  nt»r-f  •ram* 


makes  possibla  tsmporal  resolution  as  short 
as  30  seconds. — '^Spatial  data  were  obtained  by  matching  the 
inter-ferometer  ’  s  higlv^throughput  to  a  unique  optical  system 
which  includes  a  50-cm  diameter  telescope.  Th i s -r^at.-4vaiy' 
large  diameter  Dall<^irkham  telescope  maintains  the  large 
throughput  o-f  the  inter-ferometer  but  narrows  the  instrument 
-field  of  view  to  less  than  a  degree.  The  spatial  resolution 
of  the  system  is  14  mi  1 1  iradians.y 

The  interferometer  was  operated  in  conjunction  with  a 
low-light  level  infrared  imaging  isocon  camera  system 
provided  by  the  University  of  Southampton,  England.  The 
camera  was  co-aligned  with  the  telescope  to  provided  an 
infrared  video  "eye"  for  the  interferometer. 

A  bright  OH  Meinel  airglow  structure  event  was  recorded/ 
on  June  19,  1983  from  an  observation  site  located  at 
Sacramento  Peak,  New  Mexico.'^  The  structures  were  measured 
at  elevation  angles  near  the  horizon.  Apparent  wavelengths, 

f  •  i  ■ 

periods,  and  phase  velocities,  of  24±;1  km,  14f'l  minutes,  and 
28^2  meters/second  respectively,  were  calculated  for  the 
recorded  structure.  The  interferometer  data  show  intensity 
modulations  of  20^40%  within  the  structure.  Rotational 
temperatures  were  also  calculated  using  the  interferometer 
spectral  data.  A  mean  rotational  tempe.-ature  of  165  *K  was 
calculated  and  temperature  modulations  of  5-10  *K  were 
recorded  in  phase  with  the  intensity  modulations. 
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=  zenith,  day  166,  3:30-6:45  hrs.  UT  .... 
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C-24. 

OH  (7,4)  band  smoothad  rotational 
tamparatura  and  standard  daviation, 
viewing  angla  »  zanith, 

day  166,  3:30-6:45  hrs.  UT  . 
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C-2S. 

OH  (4,2)  band  relative  intensity  and 
standard  daviation,  viawing  angla  =  17*  El. 

328*  Az.,  day  166,  7:30-8:30  hrs.  UT  .... 
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C-26. 

OH  (4,2)  band  rotational  tamparatura  and 
standard  daviation,  viawing  angla  >  17*  El. 

328*  Az.,  day  166,  7:30-8:30  hrs.  UT  .... 
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C-27. 

OH  (4,2)  band  smoothad  rotational 

tamparatura  and  standard  daviation, 

viawing  angla  >  17*  El.  328*  Az. , 

day  166,  7:30-8:30  hrs.  UT  . 
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C-28. 

OH  (3,1)  band  relative  intensity  and 
standard  daviation,  viawing  angla 
*  17*  El.  328*  Az., 

day  166,  7:30-8:30  hrs.  UT  . 
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C-29. 

OH  (3,1)  band  rotational  tamparatura  and 
standard  daviation,  viawing  angle 
*  17*  El.  328*  Az., 

day  166,  7:30-8:30  hrs.  UT  . 
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C-30. 

OH  (3,1)  band  smoothed  rotational 

tamparatura  and  standard  daviation, 

viawing  angla  •  17*  El.  328*  Az., 

day  166,  7:30-8:30  hrs.  UT  . 
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OH  (8,5)  band  rclativ*  Intensity  and 
standard  davlatian,  viawing  angla 
«  17*  El.  328*  Az., 

day  166,  7:30-8:30  hrs.  UT  . 

OH  (8,5)  band  rotational  tamparatura  and 
standard  daviation,  viawing  angle 
■  17*  El.  328»  Az., 

day  166,  7:30-8:30  hrs.  UT  . 

OH  (8,5)  band  smoothed  rotational 
temperature  and  standard  daviation, 
viewing  angle  «  17*  El.  328*  Az., 
day  166,  7:30-8:30  hrs.  UT  . 

OH  (7,4)  band  relative  intensity  and 
standard  deviation,  viewing  angle 

*  17*  El.  328*  Az., 

day  166,  7:30-8:30  hrs.  UT  . 

OH  (7,4)  band  rotational  temperature  and 
standard  deviation,  viewing  angle 
»  17*  El.  328*  Az., 

day  166,  7:30-8:30  hrs.  UT  . 

OH  (7,4)  band  smoothed  rotational 
temperature  and  standard  deviation, 
viewing  angle  •  17*  El.  328*  Az., 
day  166,  7:30-8:30  hrs.  UT  . 

OH  (4,2)  band  relative  intensity  and 
standard  deviation,  viewing  angle 

*  15.5*  El.  340*  Az., 

day  166,  8:30-9:15  hrs.  UT  . 

OH  (4,2)  band  rotational  temperature  and 
standard  deviation,  viewing  angle 
«  15.5*  El.  340*  AZ., 

day  166,  8:30-9:15  hrs.  UT  . 

OH  (4,2)  band  smoothed  rotational 
temperature  and  standard  deviation, 
viewing  angle  «  15.5*  El.  340*  Az. , 
day  166,  8:30-9:15  hrs.  UT  . 

OH  (3,1)  band  relative  intensity  and 
standard  deviation,  viewing  angle 

*  15.5*  El.  340*  Az., 

day  166,  8:30-9:15  hrs.  UT  . 


OH  (3,1)  band  rotational  temperature  and 
standard  deviation,  viewing  angle 
>  15.5*  El.  340*  Az., 

day  166,  8:30-9:15  hrs.  UT  . 
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C-42 


OH  (3,1)  band  amaothad  rotational 
tamparatur*  and  standard  daviation, 
viawing  an9la  »  IS. 5*  El.  340*  Az. , 
day  166,  8:30-9:13  hrs.  UT  . 

C-43.  OH  (8,S)  band  ralativa  intanaity  and 
standard  daviation,  viawing  an^le 
«  13.5*  El.  340*  Az., 

day  166,  8:30-9:13  hrs.  UT  . 

C-44.  OH  (8,3)  band  rotational  tamparatura  and 
standard  daviation,  viawing  angla 
=  IS. 3*  El.  340*  Az., 

day  166,  8:30-9:13  hrs.  UT  . 

C-4S.  OH  (8,3)  band  smoothad  rotational 

tamparatura  and  standard  daviation, 
viawing  angla  «  IS. 3*  El.  340*  Az. , 
day  166,  8:30-9:13  hrs.  UT  . 

C-46.  OH  (7,4)  band  ralativa  intansity  and 
standard  daviation,  viawing  angla 
«  IS. 3*  El.  340*  Az., 

day  166,  8:30-9:13  hrs.  UT  . 

C-47.  OH  (7,4)  band  rotational  tamparatura  and 
standard  daviation,  viawing  angla 
=  13.3*  El.  340*  Az., 

day  166,  8:30-9:13  hrs.  UT  . 

C-48.  OH  (7,4)  band  smoothad  rotational 

tamparatura  and  standard  daviation, 
viawing  angla  »  IS. 3*  El.  340*  Az., 
day  166,  8:30-9:13  hrs.  UT  . 

C-49.  OH  (4,2)  band  ralativa  intansity  and 
standard  daviation,  viawing  angle 
*  13.3*  El.  309*  Az., 

day  166,  9:13-10:13  hrs.  UT  . 

C-SO.  OH  (4,2)  band  rotational  tamparatura  and 
standard  daviation,  viawing  angla 
■  IS. 3*  El.  309*  Az., 

day  166,  9:13-10:13  hrs.  UT  . . 

C-Sl.  OH  (4,2)  band  smoothad  rotational 

tamparatura  and  standard  daviation, 
viawing  angla  «  IS. 3*  El.  309*  Az., 
day  166,  9:13-10:13  hrs.  UT  . 


C-S2. 


OH  (3,1)  band  ralativa  intansity  and 
standard  daviation,  viawing  angla 
«  IS. 3*  El.  309*  Az., 

day  166,  9:13-10:13  hrs.  UT  .  .  .  .  . 


C-S3 


OH  (3,1)  band  rotational  temparature  and 
standard  daviation,  viawing  angia 
»  1S.5*  El.  309*  Az., 

day  166,  9:15-10:15  hrs.  UT .  194 

C-54.  OH  (3,1)  band  smoothad  rotational 

tamparatura  and  standard  daviation, 
viawing  angla  »  15.5*  El.  309*  Az . , 

day  166,  9:15-10:15  hrs.  UT .  195 

C-55.  OH  (8,5)  band  ralativa  intansity  and 
standard  daviation,  viawing  angla 

-  15.5*  El.  309*  Az., 

day  166,  9:15-10:15  hrs.  UT .  196 

C-56.  OH  (8,5)  band  rotational  tamparatura  and 
standard  daviation,  viawing  angla 

■  15.5*  El.  •'09*  Az., 

day  166,  9:  Ib-lO:  15  hrs.  UT .  197 

C-57.  OH  (8,5)  band  smoothad  rotational 

tamparatura  and  standard  daviation, 
viawing  angla  >  15.5*  El.  309*  Az., 

day  166,  9:15-10:15  hrs.  UT .  198 

C-58.  OH  (7,4)  band  ralativa  intansity  and 
standard  daviation,  viawing  angla 

■  15.5*  El.  309*  Az., 

day  166,  9:15-10:15  hrs.  UT .  199 

C-59.  OH  (7,4)  band  rotational  tamparatura  and 
standard  daviation,  viawing  angla 

-  15.5*  El.  309*  Az., 

day  166,  9:15-10:15  hrs.  UT .  200 

C-60.  OH  (7,4)  band  smoothad  rotational 

tamparatura  and  standard  daviation, 
viawing  angla  >  15.5*  El.  309*  Az., 

day  166,  9:15-10:15  hrs.  UT .  201 


CHAPTER  I 
INTRODUCTION 

The  mesosphere  is  the  interface  region  between  the 
earths’  inner  and  outer  atmospheres.  Occurrences  at  the 
mesopause  include:  the  temperature  gradient  makes  a  sign 
change,  the  atmospheric  pressure,  density,  and  mean 
molecular  weight  all  have  an  inflection  point  in  their 
respective  curves  CBanks  and  Kockarts  19731.  The  atmosphere 
makes  a  transition  from  a  fluid  to  free  molecular  flow  in 
this  region  which  accounts  for  these  changes.  The  unique 
properties  of  the  mesospheric  region  are  of  great  interest 
in  understanding  the  middle  atmosphere  and  its  influence  on 
the  energy  budget  of  the  earth. 

The  80  to  100  km  region  (mesosphere)  is  a  difficult 
region  to  study  because  the  altitude  is  too  low  for  direct 
satellite  observations  and  too  high  for  direct  balloon  or 
airplane  measurements.  Ground  based  studies  are  hampered  by 
the  intervening  atmosphere.  There  are  some  relatively 
transparent  atmospheric  "windows”  in  the  near  infrared, 
however.  Hydroxyl  radicals  (OH)  reside  in  this  mesospheric 
region  in  sufficient  concentrations  to  radiate  a  large 
amount  of  energy  at  the  red  and  near  infrared  wavelengths 
[Baker  19781. 


Th»  •xcitation  o-f  OH  is  caussd  by  various  solar  and 
chemical  processes.  The  excited  OH  radical  is  a  complex 
vibrat ional -rotat ional  system  which  emits  radiation 
spectrally.  The  spectral  radiation  distribution  is  near 
Maxwel  1 -Bol tzmann  in  nature  CBaker  19783  1  there-fore,  the  OH 
spectral  radiation  can  be  measured,  the  Boltzmann 
distribution  determined,  and  a  rotational  temperature 
calculated  CUare  19803.  This  temperature  can  then  be  used 


to  help  understand  the  chemistry  and  physics  of  the  entire 
region.  Banks  and  Kockarts  C 19733  show  that  mesospheric 
temperatures  during  the  summer  months  at  a  mid-latitude  site 
can  be  expected  to  be  between  150  *K  and  190  *K. 

Recent  studies  have  shown  that  mesospheric  optical 
radiation  called  airglow  has,  at  times,  exhibited  some 
"wavelike*  structure  CTaylor  et  al.  19803.  These  waves  have 
been  studied  using  photographic  and  photometric  methods. 
The  object  o-f  this  study  was  to  develop  and  utilize  an 
instrument  to  provide  quantitative  data  o-f  OH  rotational 
temperature  and  intensity  variations  during  these  airglow 
structure  events  . 


Airglow  Structure  Measurement  Background 


The  atmospheric  airglow  layer  has  been  observed  and 
studied  -for  many  years  using  a  variety  of  methods.  The 
studies  conducted  during  the  decades  from  1930  to  1970, 
however,  failed  to  recognize  the  nature  of  the  airglow 
structure  phenomenon.  Rayleigh  (19313  was  among  the  first 


to  rocosniz*  th*  dif-foronc*  botwoon  airglow  and  aurora.  Ha 
ra'farrad  to  tha  anhancad  airflow  as  "non-polar  aurorae." 
Photographs  oi  tha  airslow  structure  ware  presented  in  1952 
by  Ho4f mai star  1195231  but  asain,  neither  the  identity  nor 
tha  source  of  the  airflow  structure  was  understood. 
Chamberlain  C19613  brie-fly  outlines  the  historical  e-fforts 
in  airslow  studies  up  until  1961.  During  the  decade  of  the 
60s  the  techniques  of  photometry  were  perfected  and  most 
optical  atmospheric  study  efforts  were  centered  around  these 
methods. 

Kieffaber  C 19731  presented  photographic  evidence  in  1972 
of  apparent  airglow  "waves*  and  "cells’"  in  the  750-900  nm 
wavelength  region  using  infrared  film  and  a  35-mm  camera. 
She  proposed  that  the  airglow  stripes  originated  from  a 
disturbance  in  the  OH  layer.  Again  in  late  1972,  Peterson 
and  Kieffaber  119733  recorded  more  occurrences  of  structure 


at  their  mid-latitude  site  near  Albuquerque,  New  Mexico. 
The  photographically  recorded  events  were  also  tracked  with 
infrared  photometers  <1.65  and  2.2  lim)  and  shown  to  be 
moving  between  20  and  40  meters/second. 

In  1975,  Crawford  et  al.  119753  flew  an  image 
intensified  isocon  television  system  on  board  NASA’s  CV990 
aircraft.  Peterson  and  Kieffaber  119753  obser .ed  on  the 
same  flight  with  their  cameras  and  photometers.  Both  groups 
recorded  "cloud-like"  airglow  structures.  Again  using  35-mm 
cameras,  Moreels  and  Herse  119773  measured  extensive  OH 
airglow  structure  over  Europe.  Their  findings  were  similar 


to  those  of  Peterson  and  Kieffaber.  Waves  on  the  order  of 
40-km  spatial  wavelensth  appeared  to  be  travel in9  at 
horizontal  speeds  from  15  to  20  meters /second .  Peterson 
C 19791  was  able  to  record  numerous  occurrences  from  1975 
through  1978  with  some  of  the  events  being  enhanced  enough 
to  see  the  structure  with  the  naked  eye.  The  University  of 
Southampton  Atmospheric  Physics  group  CTaylor  et  al .  1980 
and  Taylor  1984]  recorded  many  structure  events  with  the 
image  isocon  television  cameras  from  1975  through  1983. 
Using  radiometric  techniques,  Huppi  and  Baker  C 19761  also 
recorded  OH  intensity  variations.  Takeuchi  and  Misawa 
C19811  record  some  short-period  waves  of  OH  intensity  and 
rotational  temperature  using  a  tilting  filter  photometer. 
This  study  used  a  narrow  field  of  view,  fast  scan  rate 
instrument.  The  spectral  resolution  was  rather  coarse, 
however,  (unable  to  resolve  the  base  line  between  adjacent 
band  lines)  and  the  measurements  were  taken  without  the  aid 
of  any  photographic  or  video  equipment  making  it  difficult 
to  identify  what  was  being  observed.  Airglow  structure  was 
also  recorded  by  Peterson  and  Adams  C 19831  during  a  total 
lunar  eclipse  in  the  summer  of  1982.  In  this  lunar  eclipse 
study  extensive  use  was  made  of  photographic  equipment  as 
well  as  a  vertical  sounding  radar. 

During  the  decade  of  the  1970’s,  i nter f erometr i c - 
spectroscopy  applied  to  middle  atmospheric  research  matured 
as  a  measurement  science.  Inter ferograms  processed  using 


computer-based  fast  Fourier  transform  (FFT)  methods  yielded 


hi9h-r»solution  OH  airglow  spectra  ti'om  which  intensity  and 
rotational  temperatures  could  be  extracted.  Baker  C 19783 
presents  an  excellent  summary  o-f  the  studies  conducted  in 
this  area.  However,  because  the  instruments  used  had  wide 
-fields  o-f  view  (there-fore  integration  over  large  areas),  low 
throughput  (there-fore  long  integration  times),  or  low 
spectral  resolution  the  small  spatial  airglow  variations 
were  not  able  to  be  spectroscopically  measured  at  high 
resolution. 

The  re-ferences  cited  and  many  others  have  recorded  OH 
airglow  structure  events.  However,  none  o-f  the  researchers 
were  able  to  provide  a  measure  o-f  high-resolution  spectral 
changes  and  there-fore  calculate  the  di-f -ferences  in  OH 
rotational  temperature  o-f  the  dark  and  bright  band  "waves." 
This  dissertation  gives  the  design,  development,  and 
operation  o-f  a  special  high-throughput,  narrow  -field  o-f 
view,  fast  scan  interferometer-spectrometer  which  can 
spectrally,  spatially,  and  temporally  resolve  the  OH  airglow 
emission  structure. 


Airglow  Structure  and  the  Theory 


of  Atmospheric  Gravity  Waves 


In  a  landmark  paper,  Hines  [19601  suggested  that  under 
certain  conditions  the  atmosphere  could  be  disturbed  by  a 
"gravity  wave."  Later,  Hines  [19653  hypothesized  that  the 
passage  of  internal  gravity  waves  (IGW’s)  through  the 
atmosphere  would  cause  some  reversible,  adiabatic  heating 


from  5  to  2S  *K  with  the  high  vibrational  level  transitions 


appearing  hotter. 


Interferometers  and  Fourier 
Transform  Spectroscopy 

The  Michelson  interferometer  was  invented  in  the  1880’ s 
by  Albert  Abraham  Michelson  CShankland  19741.  In  the 
Michelson-Mor ley  experiment,  the  instrument  was  used  in  an 
[  attempt  to  measure  the  earth’s  movement  through  an  "ether." 

I  Michelson  also  used  his  interferometer  to  determine  the 

I  exact  length  of  the  standard  meter  and  to  measure  the 

diameter  of  celestial  bodies.  He  also  discovered  the 
spectral  fine  structure  of  hydrogen,  mercury,  and  thallium 
I  CShankland  19741.  This  pointed  out  the  potential  for  what 

would  be  later  be  called  "Fourier  spectroscopy." 

Opt i cal ly-sensi t i ve  detectors  were  subsequently  used  in 
'  conjunction  with  Michelson  interferometers,  in  which  one 

mirror  was  mechanically  displaced  at  a  constant  rate,  to 
produce  an  electrical  i nterf erogram .  The  i nterf erogram  was 

i 

then  inverted  using  Fourier  transform  techniques  to  yield 
direct  spectral  data.  In  1910,  Ruben  and  Woods  CConnes 

19631  obtained  the  first  far  infrared  spectrum  using  this 

! 

method.  Fellgett  C19491  and  Jacquinot  C19S41  independently 

t 

t 

showed  the  inherent  advantage  that  the  interferometer  has 

over  grating  and  prism  spectrometers.  This  advantage 

I 

:  results  from  measuring  all  spectral  components 

simultaneously.  This  improvement  was  referred  to  by 


MIRROR 


Figure  1-1.  Layout 
i  nter-f  erometer . 


INTERFEROGRAM 


conventional 


Michelson 


that  the  path  length  AC  is  the  same  as  the  path  length  AB, 
then  the  recombined  signals  are  in  phase  and  thus  add 
constructively.  The  same  constructive  addition  occurs  when 
the  path  length  di-f-ference  BC-AC  is  any  integral  number  of 
wavelengths  of  the  incoming  signal.  On  the  other  hand,  if 
the  path  length  difference  is  not  an  integral  multiple  of 
the  wavelength,  then  the  recombined  signal  will  have  varying 
amounts  of  destructive  interference  depending  upon  the  phase 
difference.  As  the  path  length  AC  is  changed  in  a  uniform 
manner,  by  moving  at  a  constant  rate,  the  electrical 
signal  from  the  detector  is  the  i nterf erogram  of  the 
incoming  optical  signal.  The  Fourier  transformed 


i nterf erogram  yields  the  spectral  content  of  the  incoming 


iO 

The  simultaneous  measurement  o-f  h  i  gh-r  eso  lut  ion 
spectral,  temporal,  and  spatial  characteristics  o-f  the  OH 
airgloM  structure  requires  an  instrument  with  both  a  narrow 
■field  o-f  view  and  high  throughput.  The  standard  Michelson 
interferometer,  when  used  for  high-resolution  measurements, 
has  a  narrow  field  of  view  but  its  low  throughput  would  make 
it  an  order  of  magnitude  less  sensitive. 

The  narrow  field  of  view  limitation  of  a  standard 
Michelson  interferometer  is  illustrated  in  Figure  1-2.  When 
incoming  energy  is  allowed  to  enter  the  interferometer  off- 
axis  <di*0*  )  the  relationship  between  the  displacement  of 
mirror  M^  and  the  actual  path  difference  between  the  two 
mirrors  is  altered.  The  path  difference  or  retardation  is 
no  longer  as  it  is  when  light  is  coming  straight  into  the 
instrument,  but  now  is  2a'cose ,  where  9  is  the  angle  of  the 
incoming  light  with  respect  to  the  entry  normal.  For  an 
instrument  with  a  given  resolving  power  the  maximum  usable 
field  of  view  for  a  standard  Michelson  interferometer  is, 
according  to  Vanasse  C 19771 

n  *2« IR  ,  (1.1) 

max 

where 

A  ^  maximum  field  of  view  in  steradians, 
max 

R  *  resolving  power  of  instrument. 

To  increase  the  throughput  (thus  achieving  a  faster  scan 
rate)  an  optically-compensated  interferometer  was  chosen  for 


use  in  this  study 


There  have  been  many  proposed  methods 
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I 


of  -f  ield-Midenin9  or  opt  leal -componsat  Ion  for  increasing  the 
throughput  of  an  interferometer.  These  techniques  are 
reviewed  by  Baker  CVanasse  19773.  The  method  used  for  the 
instrument  in  this  study  was  first  proposed  by  Connes 
[19S63,  and  uses  opt ical -compensat i on  wedges  or  prisms  in 
each  leg  of  the  interferometer  (specific  details  are 
discussed  in  Chapter  ID.  Optical  compensation  increases 
the  throughput  by  increasing  the  usable  field  of  view  of  the 
instrument.  However,  the  measurement  of  OH  airglow 
structure  requires  a  small  field  of  view.  The  high 
throughput  of  the  compensated  interferometer  was  matched 
(maintaining  temporal  resolution)  to  a  special  optical 
system  which  included  a  large-diameter  telescope  to  obtain 
the  desired  narrow  field  of  view  while  maintaining 
throughput . 

The  interferometer  system  was  now  able  to  simultaneously 
resolve  the  spectral,  temporal,  and  spatial  characteristics 
of  the  OH  airglow  structure.  However,  because  the  measured 
radiation  was  in  the  infrared  it  was  necessary  to  locate  and 
measure  a  structure  occurrence  and  to  characterize  the  total 
structure  into  which  the  interferometer  was  looking.  The 
video  viewing  system  used  is  described  in  the  next  section. 


Isocan  Television  System 


Taylor  [1983-841  participated  in  this  research  by 
providing  and  operating  a  low  light-level  infrared 


television  camera  used  in  conjunction 


lith  the 


interferometer-spectrometer.  This  Southampton  University 
infrared  TU  camera  allowed  any  airglow  structure  to  be 
quickly  and  efficiently  identified  and  permanently  recorded 
as  video  information.  Taylor’s  system  employed  an  image 
intensified  isocon  television  system.  This  isocon 
television  system  is  briefly  described  here  for 
comp  1 eteness . 

The  isocon  tube  is  different  than  other  scanning  image 
tubes  in  that  it  uses  a  different  portion  of  the  scanning 
electron  beam  to  create  the  signal.  The  low-energy  electron 
beam  scans  a  high-resistance  target  as  in  other  tubes; 
however,  the  isocon  signal  comes  from  the  electrons  that  are 
scattered  off  the  target  CSoule  19681  rather  than  the 
reflected  electrons  used  by  conventional  equipment 
(orthicons  use  the  reflected  beam).  The  scattered  signal, 
although  small  in  magnitude,  has  a  high  si gnal - to-noi se 
ratio,  and  is  particularly  well  suited  to  the  viewing  of  the 
low-contrast,  faint  OH  airglow. 

The  TV  system  used  to  image  the  airglow  structures  was 
specially  modified  to  enable  clear  images  of  the  OH  airglow 
patterns  to  be  obtained  in  less  than  a  second  [Taylor  19841. 
The  television  camera  used  was  an  English  Electric  Valve 
Miniature  Isocon,  Type  P1477  fitted  with  a  single-stage 
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image  i  ntensi-f  ier ,  optically  coupled  to  a  55-mm  image  isocon 
tube.  The  camera  has  a  signal -to-noise  ratio  of 
approximately  40  dB  at  10  -ft-candles  (starlight  conditions) 
and  a  dynamic  range  of  about  2000: 1. 

To  further  enhance  the  capability  of  the  camera  to  image 


very  faint  airglow  structures,  it  was  arranged  for  the 
electronic  image  to  be  integrated  on  the  target  of  the  TV 


I 


tube  for  a  period  of  up  to  a  second  (longer  integration 
times  allow  the  image  charge  on  the  target  to  migrate  and 
thus  smear  the  image)  before  being  scanned  and  recorded  onto 
video  tape  CTaylor  19841.  This  technique  is  particularly 
useful  as  it  improves  the  si gnal - to-noi se  ratio  of  the 
airglow  signal  by  nearly  an  order  of  magnitude  17  times  for 


u 


a  1-second  Integration  period)  with  no  significant  loss  of 
temporal  resolution. 

The  camera  has  an  extended  red  spectral  response  (S25) 
with  a  peak  sensitivity  at  about  SOO  nm  and  a  long 
wavelength  cut  off  at  900  nm.  Images  of  the  near  - i nf r ar ed 
OH  structure  were  obtained  by  placing  a  Schott  RG71S  band 


stop  filter  in  front  of  the  camera  lens.  The  combined 
response  of  this  filter  and  the  TV  camera  gave  a  bandwidth 
(half  maximum)  of  715  to  850  nm  and  a  peak  sensitivity 
around  750  nm.  This  bandwidth  is  illustrated  in  Figure  1-3. 
The  location  of  all  the  OH  emission  bands  within  this 
spectral  range  are  indicated.  The  intensity  in  the  zenith 
of  the  OH  emission  within  this  bandpass  is  typically  5  to  15 
kR|  the  principal  emissions  are  the  OH  (9,4)  and  (5,1) 


Th*  Southampton  TV  camaraa  hava  baan  used  -for  many  yaars 


to  "photograph”  tha  naar-in-frarad  OH  airglow  atructura.  An 
axampla  o-f  tha  quality  o-f  tha  vidao  data  gatharad  is  shown 
in  Figura  1-4.  This  vidao  -frama  was  takan  in  August  of  1980 
whila  obsarving  ovar  tha  Swiss  Alps.  The  bright  and  dark 
bands  each  subtend  about  1*  of  arc.  These  data  ware  tha 
basis  upon  which  tha  i  ntarf  aroma'tar  ’  s  field  of  view  was 
designed  in  order  to  resolve  tha  spatial  nature  of  tha  OH 
airglow  atructura. 


Scope  and  Objectives 

Tha  specific  goals  and  objectives  are  outlined  as 
f ol lows: 

1.  Design  and  develop  an  optical  instrumentation  system 
capable  of  quantifying  tha  spatial,  spectral,  and 
temporal  characteristics  of  OH  near- i nf rared  night 
airglow  structure.  Tha  instrumental  field  of  view  must 
be  one  degree  or  less  to  resolve  the  structural 

characteristics  of  the  airglow.  The  system  must  have 
spectral  resolution  of  better  than  3  cm  ^  in  order  to 
provide  spectra  from  which  OH  rotational  temperatures 
can  be  calculated  using  appropriate  algorithms  and 
digital  computer  programs.  The  NESR  must  be  sufficiei.t 
to  resolve  the  OH  near-infrared  airglow  with  scan  times 
of  less  than  1  minute  to  resolve  the  temporal 
fluctuations  of  the  airglow  structure. 


.V.-.'. 
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Us*  th*  i nstrumentat ion  syBt*m  to  moasur*  th*  spectral, 
spatial,  and  temporal  characteristics  o-f  OH  near- 
in-frared  airglow  structure  -from  a  mid-latitude  observing 


Develop  and  apply  signal  processing  procedures  extending 
the  sampling  and  FPT  work  o-f  Ware  [19801  to  extract  both 
radiance  and  rotational  temperature  variations  o-f  OH 
airglow  structure. 

Derive  error  bounds  on  the  measurement  data,  based  upon 
system  spec  i-f  i  cat  ions  such  as  -field  o-f  view,  scan  speed, 
and  spectral  resolution,  as  well  as  instrument 
calibration,  and  signal  processing  techniques. 

Present  th*  observational  results,  correlate  intensity 
variations,  temperature  -fluctuations,  and  structure  with 
simultaneous  near-in-frared  video  images,  and  compare  the 
findings  with  expected  OH  airglow  dynamics  studies  from 
other  investigators. 
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CHAPTER  II 

OPTICAL  INSTRUMENTATION  SYSTEM  DESIGN 

Design  Philosophy 

The  goal  o-f  this  study  was  to  develop  a  technique  -for 
simultaneously  measuring  the  spatial,  spectral,  and  temporal 
characteristics  O'f  OH  near- i nf rared  airglow  structures.  The 
basic  instrumental  approach  chosen  to  provide  the  spectral 
resolving  capability  is  a  Michelson  interferometer- 
spectrometer  which  is  optically-compensated  to  achieve  a 
very  high  throughput.  The  compensation  technique  used  makes 
it  passible  for  obliquely  incident  optical  energy  up  to  5 
degrees  off  axis  to  contribute  to  the  detected  signal 
without  sacrificing  spectral  resolution.  This  resulting 
high  throughput,  within  the  interferometer,  is  matched  at 
the  input  to  a  large  diameter  collecting  telescope  yielding 
a  sub-degree  field  of  view  needed  for  spatial  resolving 
power.  The  entire  optical  system  is  d iagrammat ical ly  shown 
in  Figure  2-1. 

An  optically-compensated  interferometer  has  the  high 
throughput  needed  to  achieve  a  relatively  high  temporal 
resolving  power,  in  other  words,  a  scan  time  of  less  than  a 
minute.  Temporal  variations  are  therefore  identified  while 
maintaining  a  spectral  resolution  of  2  cm  ^  in  the  near- 
infrared.  A  spectral  resolution  nearly  this  high  is 


desirable  in  order  to  unambiguously  compute  OH  rotational 


I 


temperatures  -from  the  measured  spectra.  A  -free  spectral 

range  o-f  0.8  to  1.6  Mm  has  been  obtained  using  a 

cryogenical ly-cooled  intrinsic  germanium  detector.  The 

optical  system  was  designed  to  maintain  the  high  throughput 

2 

capability  (0.285  cm  sr  >  o-f  the  i  nter-f  erometer  while 
operating  at  a  narrow  -field  o-f  «/iew  (<1*>  in  order  to  be 
able  to  resolve  the  spatial  nature  o-f  the  airglow  structure. 

The  design  criteria  and  the  resulting  design  -for  an 
optically-compensated  i  nter-f  erometer  are  given  in  this 
chapter.  Then  an  analysis  is  made  o-f  the  resulting 

i n ter f erometer -spec trometer  system. 

High-Throughput  Inter-f erometer  Design 

The  optical  layout  o-f  a  conventional  Michelson 

interferometer  was  discussed  in  Chapter  I.  Referring  to 

Figure  1-2,  the  retardation  or  path  difference  is  a  function 

of  the  entry  angle  of  the  incoming  energy.  The  relationship 

is,  a=2</cos9,  where  d  is  the  on-axis  (9=0*)  drive  distance 

and  9  is  the  angle  of  the  off-axis  ray.  The  maximum  field 

of  view  for  a  standard  Michelson  interferometer  is, 

A  =2ii (Eq  1.1),  where  A  is  the  maximum  field  of  view 
max  max 

of  the  instrument  and  R  is  the  resolving  power.  When  viewing 
faint  airglow  events  this  limitation  on  throughput  is  a 
severe  one.  The  optical  system  is  compensated  to  increase 
system  throughput  for  temporal  resolution.  This  high- 


throughput  has  been  matched  to  a  collector  system  to 


MIRROR  SURFACE 


Figure  2-2.  Optical  compensation  method  conceived  by  Connes 
C 19563 . 

maintain  the  throughput  while  narrowing  the  -field  o-f  view 
(-for  spatial  resolution). 

The  method  used  in  this  study  to  improve  the  throughput 
o-f  the  interferometer  system  is  the  one  analyzed  by 
Bouchareine  and  Connes  C 19631  and  is  depicted  in  Figure  2-2. 
The  compensated-optics  design  used  follows  that  developed  by 
Despain  et  al .  £19713,  and  the  design  limitations  previously 

derived  will  be  summarized  here  for  completeness.  Figure  2- 
3  shows  the  compensation  analysis  approach  of  Steed  £19781. 
With  the  optical  compensation  prisms  inserted,  the 


retardation  is 
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2c^  CQS  e  *  2tn  cos  ^  -  2t  cos  9  ,  (2.1) 


Mher« 


rstardation  with  compansat ion  prisms  insartad, 
t  -  thickness  o-f  optical  material, 
n  »  index  o-f  ra4raction  o-f  optical  material, 

<f  s  mirror  drive  distance,  AC  -  AB, 

©  =»  angle  o+  0+4 -ax is  ray, 

4  =  re4raction  angle  o4  o44-axis  ray  through  optical 
material . 


The  retardation  equation  shown  above  can  be  expanded  in 
a  Taylor  series  and  like  terms  collected  to  show  the  4ield 
o4  view  dependency  more  directly)  namely. 


2[t<rt  -!)♦</]♦[#  ”  "  -  d]© 

-  »'  .  _£_^  .  ^l,-' 
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(2.2) 


which  is  a  quartic  equation  in  ©. 

Ideally,  the  retardation  should  be  independent  o4  the 
entry  angle  9.  Equation  2.2  shows  that  this  is  not  possible 
by  varying  only  the  drive  distance  d  and  the  index  o4 
re4raction  n.  However,  a  signi+icant  improvement  can  be 
realized  by  designing  such  that. 


^  -  ‘[^] 


(2.3) 


which  eliminates  the  9  term  in  the  retardation  Equation  2.2 

2 

by  4orcing  the  9  coe44icient  to  zero.  This  leaves  only 


Figure  2-3.  Optical  compensation  with  an  optical  section 
inserted  into  one  leg  o-f  a  Michelson  interferometer  CSteed 
19783. 


the  e  term  and  since  6  is  small  this  term  is  relatively 
insignif leant. 

Designing  the  compensation  wedges  to  satisfy  Eq  2.3  and 
then  solving  Equation  2.2  for  an  on-axis  ray  (e^O* ) ,  the 


resulting  retardation  a  ii 

CO 


a  - 

CO  •* 


(2.4) 


Assuming  a  monochromatic  wave  front  and  using  the 
results  of  Eq  2.4,  the  general  retardation  for  the 
compensated  case  is 


a  +  a  . 

CO  CO 


a 


(2.5) 


1. 


Mhmr* 

s  compensatad  ratardation, 

&  s  retardation  tor  an  on-axia  ray. 

CO 

d  a  entry  angle  -for  incoming  ray  in  radians 

n  =  index  o-f  re-fraction  o-f  compensation  prisms. 

The  analysis  -for  the  Connes  [19561  method  shown  in 
Figure  2-2,  in  which  only  one  optical  component  is  driven, 
is  the  same  as  -for  the  system  shown  in  Figure  2.3.  The 
thickness  o-f  the  optical  material  must,  however,  increase 
with  the  drive  distance  </  in  order  to  maintain  optical 
compensation  (see  Equation  2.3).  The  increase  in  optical- 
material  thickness  needed  to  maintain  compensation  is 
obtained  by  driving  one  o-f  the  optical  wedge/mirror 
assemblies  in  synchronism  with  the  drive  motor.  Re-ferring 
to  Figure  2-2,  it  can  be  seen  that  the  re-flective  elements 
or  mirrors  in  this  method  are  created  by  depositing  the 
mirrored  sur-face  to  the  back  side  o-f  each  wedge. 

Steed  [19781  derived  the  limits  on  -field  o-f  view  as  a 
-function  of  resolving  power  based  upon  aberration  limits. 


The  results  are  categorized  into  several  groups.  The  first 
is  chromatic  aberration,  that  is,  differences  in  the 
compensation  because  the  index  of  refraction  of  the  optical 
material  is  a  function  of  wavelength.  The  chromatic  limit 


nin  -  1 )  X 


Sn  t 


(2.6) 


a  -  •  • 


Mh«r* 


>  -fisld  o-f  vi»w  at  wavelength  X  in  steradians, 

*  field  o-f  view  -for  conventional  Micheleon 
M 

i  nter-f  erometer , 

n  =  index  o-f  re-fraction  at  wavelength  X, 

S/7  =  n—n  x/Aere  n  »  index  o-f  re-fraction  at  compensated 

c  c 

wavelength, 

R  s  desired  resolving  power, 

d  =  retardation  -for  on-axis  ray, 

CO 

X  =  wavelength  expressed  in  same  units  as 
The  second  limiting  -factor  on  -field  o-f  view  is  that  o-f 
spherical  aberration.  This  limit  is  also  derived  by  Steed 
C1978].  The  maximum  -field  o-f  view  assuming  only  spherical 
aberrations  is 


n  =  n  4  t2/f) 

S  M 

• 

(2.7) 

The 

third  type 

O-f  aberration  is 

astigmatism. 

In 

the 

Connes 

method  a 

wedge  is  placed 

in  each  leg 

o-f 

the 

i  nter-f  erometer ;  however,  only  one  o-f  the  wedges  is  driven. 
This  simpli-fies  the  mechanical  design  but  because  a  wedge  is 
in  each  optical  path,  astigmatism  aberrations  occur.  The 
severity  o-f  the  aberration  increases  as  the  angle  (k)  o-f  the 
compensation  wedges  increases.  Bouchareine  and  Connes 
C 19633  quanti-fy  this  aberration  as 

"  2n^/tan^0  ,  (2.8> 


wh«r« 


a  >  prism  an9ls  of  the  compensation  wedge. 

The  aberrations  considered  were  each  derived  assuming 
the  net  ef-fect  of  the  distortion  was  a  shift  of  one  fringe 
in  the  interference  pattern.  These  degradations  prove  to  be 
the  limiting  factors  upon  the  maximum  usable  field  of'  view. 
Figure  2-4  is  a  plot  of  Equations  2.6,  2.7,  and  2.8  and 
shows  the  relationship  between  these  various  limits.  At 
least  an  order  of  magnitude  improvement  in  throughput  is 
obtainable  using  the  wedge  prism  compensation  technique. 

The  next  step  in  the  system  design  is  to  ascertain  the 
optical  retardation  for  the  Connes  [19561  method.  Referring 
to  Figure  2-5,  Steed  C19781  showed  the  drive  plane  is 
parallel  with  the  apparent  image  plane  of  the  optical 
wedges.  The  difference  in  optical  path  length  or 
retardation  between  image  points  A1  and  A2  separated  by  a 
drive  distance  X  is 


A  »  2X  5in(p-Y) 


(2.9) 


where 

P*sin^(nsina), 

V  *  0  -  tan  ^  ”  ■" — tana], 

a  a  wedge  angle,  f 

X  »  drive  distance. 


f 


Figur*  2-4.  Plot  o-f  aborration  limita  for  Connao  C 19561 
mathod  o-f  optical  compansation  comparad  with  convantional 
Michalaon  i ntarf aromatar  CStaad  19781. 
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Driv*  System 

Examination  o-f  Equation  2.9  ahowa  that  tha  choica  o-f  a 
large  wedge  angle  reducea  the  actual  drive  diatance,  thereby 
eaaing  the  mechanical  drive  requirements.  However,  the 
choice  must  be  made  in  conjunction  with  the  limits  on 
distortion  presented  in  Figure  2-4.  The  wedges  and  cube 
beamsplitter  used  in  this  instrument  are  made  O'f  a  high 
quality  quartz,  namely,  In-frasil  X,  manu-f actured  by  Amersil 
Corporation.  The  measured  index  o-f  refraction  for  this 
material  is  1.4S.  A  compensation  prism  wedge  angle  of  8* 
was  chosen.  Salving  Equation  2.9  for  this  wedge  angle,  then 
p  =  li.64«,  Y  *  3.S0»,  and  a  »  0.21X.  The  optical 
retardation  for  this  compensation  technique  is  only  about 
20%  of  the  drive  distance,  where  as  in  a  conventional 
Michelson  interferometer  the  retardation  is  equal  to  2d 
where  d  is  the  drive  distance.  Thus,  for  a  given  spectral 
resolution  the  compensated  interferometer  requires  a  drive 
distance  of  2/. 21=9. S  times  longer  than  the  conventional 
Michelson  approach. 

An  optically-compensated  interferometer  with  a  spectral 
resolution  of  2  cm  ^  requires  a  slide  movement  of  3.1  cm. 
This  comparatively  large  drive  distance  is  accomplished 
using  a  gas- lubr icated  platform,  providing  near-zero 
friction,  developed  by  Haycock  C197S].  The  platform 
translates  the  optical  wedge/mirror  assembly  and  as  such 
must  maintain  mechanical  tolerances  as  close  as  passible  to 
the  wavelength  dimensions  of  the  infrared  light  being 


measured.  The  gas  bear  ins  surfaces  are  lapped  to  a 

tolerance  of  about  1  wavelensth  of  5461  A  lisht  and  operate 

-4 

with  a  clearance  of  2.5x10  cm.  At  an  operatins  pressure  of 
6-12  psi  the  bearins  Mill  maintain  an  optical  alisnment  of  1 
arc  second  while  translatins  5  cm  [Haycock  19751.  The 
optical  components  are  then  mounted  on  the  bear i ns~suppor ted 
optical  platform.  Adjustment  for  parallelism  is 

accomplished  on  an  optical  bench  usins  a  laboratory  He-Ne 
laser  as  a  light  source.  The  material  used  in  the  platform 
and  bearing  assembly  is  a  specially  formulated  Invar  alloy, 
chosen  for  its  excellent  temperature  and  stability 
characteristics  and  to  match  the  thermal  expansion 
characteristics  of  the  optical  material.  This  sophisticated 
platform  provides  the  relatively  long  drive  distances 
required  (up  to  5  cm)  as  well  as  the  mechanical  accuracy 
needed  to  translate  the  optical  components  properly. 

The  platform  is  driven  by  a  "voice  coil”  motor  and 
therefore  makes  no  physical  contact  with  the  rest  of  the 
instrument.  The  motor  is  driven  by  a  standard  servo- 
amplifier  with  feedback  from  positional  and  velocity  sensors 
located  within  the  platform  slide.  The  main  signal  driving 
the  slide  is  generated  by  a  digitally-controlled  ramp  which 
can  be  adjusted  for  slide  velocity  and  drive  distance.  The 
slide  controller  can  produce  variable  scan  times  of  5 
seconds  up  to  several  minutes,  and  a  drive  distance  of  up  to 
5  cm.  The  completed  interferometer  uses  a  10-cm  cube 
beamsplitter  and  li.4-cm  diameter  end  wedges  with  a  prism 


detectors  with  the  necessary  detectivity  are  very  di-f-ficult 


to  manu-Facture;  there-fore,  some  compromise  was  necessitated. 
The  required  detector  diameter  and  the  speed  o-f  the 
col  lectins  optics  shown  in  the  example  are  both  somewhat 
impract ical . 

A  detector  system  was  chosen  which  was  readily  available 
and  had  proven  e-f-fective  in  the  past.  The  detector  selected 
was  an  RCA  Ltd.,  solid  state  sermanium  device  with  a  sel-f- 
contained  preamp  1  i-f  ier  and  load  resistor  in  a  liquid- 
nitrogen  dewar.  Figure  2-7  is  a  picture  o-f  the  detector 
dewar .  The  detector  is  S  mm  in  diameter  and  has  a  noise 
equivalent  power  (NEP)  o-f  1.12  x  10  W/iHz  at  1.27  Mm  and 
covers  a  spectral  range  o-f  0.8  to  1.6  Mm. 

It  was  decided  that  a  scan  rate  o-f  30  seconds  would  be 
su-f-ficient  to  resolve  the  temporal  variations  o-f  -first 
interest  in  the  OH  airglow  structure.  This  decision  was 
based  upon  video  data  gathered  by  Taylor  et  al .  C 19801.  A 
30-second  scan,  a  minimum  wavelength  o-f  0.8  Mm,  and  a  drive 
distance  o-f  3.1  cm  <spectral  resolution  o-f  2  cm  set  the 
detector  electrical  bandwidth  at  150  Hz.  The  RCA  detector 
has  a  bandwidth  o-f  600  Hz.  Appendix  A  contains  the  detailed 
spec  i -f  i cat  ions  of  the  detector  system. 

Telescope  Design 

The  spatial  characteristics  of  the  airglow  structure, 
shown  by  Taylor  [1983-841,  dictated  that  the  instrument 
field  of  view  be  less  than  1*  full  angle.  The  goal  of  the 
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optical  dasign  was  to  maintain  the  inherent  high  throughput 


of  the  compensated  interferometer  throughout  the  entire 
system.  All  portions  of  the  optical  path  were  examined  in 
order  to  ascertain  where  the  system  was  throughput - 1 imi ted . 
The  choice  of  a  detector,  and  practical  limitations  on  the 
f-number  of  available  condenser  lenses  proved  to  be  the 
limiting  factors.  An  optical  system  was  designed, 
incorporating  the  interferometer,  using  a  large  diameter 
collector,  which  maintained  the  system  throughput  while 
narrowing  the  field  of  view  to  under  1*.  The  instrument 
system  was  now  capable  of  measuring  the  temporal,  spectral, 
and  spatial  variations  of  the  airglow  layer. 

The  design  of  the  optical  system  began  with  the  detector 

and  choice  of  its  condenser  lens  and  proceeded  back  towards 

the  telescope  collector.  The  entire  optical  system  is 

described  d i agrammat ical 1 y  in  Figure  2-1.  The  diameter  of 

the  condenser  lens  must,  however,  allow  the  converging  beam 

to  pass  through  the  beamsplitter  and  wedges  without 

vignetting.  A  commercially  available  condenser  lens  with  a 

diameter  of  68  mm  and  an  effective  focal  length  of  50  mm  was 

chosen.  The  condenser  lens  and  detector  diameter  of  5  min 

2 

set  the  system  throughput  or  Atl  at  0.28  cm  sr .  Referring  to 

2 

Figure  2-i,  a  throughput  of  0.28  cm  sr  translates  to  a  beam 
diameter  entering  the  interferometer  of  92  mm  and  converging 
in  a  6*  full  field  of  view.  The  physical  dimensions  of  the 
detector  optics  outlined  here  fit  well  within  the  size 
constraints  of  the  basic  interferometer.  The  field  of  view 
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is  also  within  th*  aberration  limits,  shown  in  Figure  2-4, 
-for  an  optically-compensated  interferometer. 

The  calculated  throughput  can  now  be  used  to  design  the 

telescope  needed  to  narrow  the  field  of  view  to  1*  or  less. 

2 

A  throughput,  dictated  by  the  detector  system,  of  0.28  cm  sr 
and  a  1*  field  of  view  design  goal,  set  the  collector 
diameter  at  D  >  AA /2v ( i -cosd ) »44  cm.  A  50.8-cm  <20-inch) 

diameter  system  was  chosen  because  of  availability  and  to 
allow  for  some  error  in  aligning  the  telescope  optics. 
Taking  advantage  of  the  larger  collector,  the  field  of  view 
was  narrowed  to  0.8*  requiring  a  48-cm  diameter  collector, 
within  the  S0.8-cm  mirror  size  and  still  allowing  for  some 
error  in  optical  system  alignment.  For  portability,  an  f/2 
primary  was  specified  to  minimize  the  telescope  length.  An 
overall  system  f  number  of  8  was  selected  to  transfer  the 
telescope  image  to  the  proper  position. 

In  this  application,  the  purpose  of  the  telescope  is  to 
gather  energy  and  transfer  it  to  the  detector,  rather  than 
to  transfer  a  spatial  image.  Therefore,  the  quality  of  the 
optical  image  within  the  field  of  view  is  not  of  as  great  a 
concern  as  it  would  be  in  an  imaging  system.  Many  types  of 
folded  telescope  systems  were  considered.  In  order  to 
simplify  the  optics  Driscoll  and  Vaughan  C 19731  suggests  a 
Dall-Kirkham  type  because  of  the  spherical  secondary,  if  the 
resulting  image  distortion  is  acceptable.  The  distortion 
for  a  Dall-Kirkham  system  of  this  size  was  calculated  using 
the  formula  in  Driscoll  and  Vaughan  119731,  and  found  to  be 


less  than  0.4%  o-f  the  area  (coma  and  astigmatism  distortion 
were  calculated  in  terms  o-f  primary  mirror  area  that  will  be 
degraded)  o-f  the  primary  mirror.  There-fore,  a  Dall-Kirkham 
type  o-f  telescope,  being  more  than  adequate,  was  chosen  for 
its  simplicity  and  relatively  low  cost.  This  type  o-f 
telescope  has  a  spherical  secondary  mirror  and  an  elliptical 
primary  mirror.  The  telescope-equipped  i  nter-f  erometer 
system  with  the  isocon  camera  mounted  on  the  telescope  is 
pictured  in  Figure  2-9. 

The  last  step  in  the  optical  design  was  to  provide  an 
optical  inter-face  between  the  telescope  and  the  detector 
sub-systems.  Two  primary  considerations  were  given  emphasis 
in  this  inter-face  design:  (1.)  imaging  the  detector  on  the 
primary  mirror,  and  (2)  imaging  the  condenser  lens  on  the 
telescope  -focal  plane.  The  detector  may  have  an  uneven 
response  across  its  area,  there-fore  imaging  the  detector  on 
the  primary  mirror  minimizes  the  e-f-fects  o-f  o-ff-axis  rays  by 
illuminating  the  entire  detector  by  light  entering  within 
the  field  of  view  of  the  instrument.  An  image  of  the 
detector  condenser  lens  at  the  telescope  focal  plane  allows 
for  Independent  control  of  the  system  field  of  view. 

An  adjustable  iris  was  then  placed  at  the  focal  plane 
for  adjusting  the  field  of  view  while  still  allowing  the 
entire  detector  to  be  illuminated  by  energy  within  the  field 
of  view.  Careful  choice  of  physical  dimensions  and  optics 
allow  the  transfer  of  each  image  independent  of  the  other. 


The  goal  o-f  placing  a  detector  image  on  the  primary 


mirror  was  accomplished  by  the  use  o-f  two  lenses.  Re-ferring 
to  Figure  2-1,  a  trans-fer  lens,  which  has  a  diameter  o-f  100 
mm  and  a  -focal  length  o-f  250  mm,  is  placed  near  the 
i nter-f er ometer  beamsplitter.  This  lens  produces  an  image  o-f 
the  detector  at  a  -focal  point  near  the  physical  entrance  to 
the  i  nter-f  erometer  housing.  The  magn  i -f  icat  ion  ratio  o-f  this 
trans-fer  lens  -focal  length  to  the  -focal  length  (50  mm)  o-f 
the  condenser  lens  determines  the  detector  image  size  at 
this  point;  there-fore,  the  detector  image  is  25  mm  in 
diameter.  A  -field  stop  is  placed  at  this  point  to  limit  the 
detector  image  size  throughout  the  rest  o-f  the  system. 
Another  trans-fer  lens,  which  has  a  75-mm  diameter  and  a 
-focal  length  o-f  200  mm,  is  placed  at  the  -focal  plane  o-f  the 
telescope.  This  lens  trans-fers  the  detector  image  at  the 
-field  stop  onto  the  primary  mirror.  As  can  be  seen  -from 
Figure  2-1,  two  angled  mirrors  were  required  to  meet  the 
physical  constraints  of  the  design. 

The  second  design  abjective  o-f  the  optical  trans-fer 
system  was  to  place  an  image  of  the  detector  condenser  lens 
on  the  telescope  focal  plane.  A  collimated  image  of  the 
condenser  lens  is  generated  by  the  100-mm  diameter  transfer 
lens  that  is  placed  near  the  interferometer  beamsplitter. 
The  addition  of  a  75-mm  diameter,  200-mm  focal  length  lens 
placed  at  the  field  stop  focuses  this  image  at  the  telescope 
focal  plane.  The  placement  of  this  lens  at  a  detector  image 
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point  does  not  affect  that  image 


An  adjustable  iris  was 
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then  placed  at  the  telescope  -focal  plane  for  control  of  the 
field  of  view.  The  field  stop  and  the  iris  provide  the 


limiting  apertures  for  the  entire  system.  The  optical 
system  described  here  was  verified  by  Harris  C1984]  using  a 
computer-aided  optical  ray-tracing  program. 

The  optical  lenses  used  in  the  design  are  made  of 
commerc ial -grade  optical  glass  and  have  no  coatings.  The 
beamsplitter  and  wedges  are  constructed  of  quartz  and  are 
not  coated.  The  primary  and  secondary  telescope  mirrors  are 
aluminized  reflection  surfaces  with  a  SiO^  coating. 


Instrument  Housing 


The  i nterf erometer-spectrometer  with  its  associated 
optical  and  telescope  systems  was  placed  in  a  30-inch 
diameter,  45-inch  tall,  round  container.  The  interferometer 
placed  in  the  instrument  housing  is  shown  in  Figure  2-8. 
The  telescope  attached  to  the  container  is  shown  in  Figure 
2-9.  The  large  container  size  was  chosen  to  support  the 
telescope  in  a  stable  manner.  The  base  of  the  package  was 
mounted  upon  a  rotating  stand  which  provides  the  ability  to 


move  the  telescope  in  azimuth.  The  rotating  base  has 
attached  three  large  pneumatic  tires  for  system  mobility  and 
three  leveling  screws  to  stabilize  the  interferometer  for 
operation.  The  telescope  is  mounted  to  the  side  of  the 
container  with  a  lO-inch  diameter  ball-bearing  allowing  the 


telescope  to  be  rotated  in  elevation.  The  two  rotating 
joints  provide  the  telescopy  with  complete  pointing  freedom. 
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The  spec i -f  i cations  o-f  the  instrument  with  the  telescope 


are  given  in  Table  2-1. 


TABLE  2-1.  Summary  04  i  nter-f  erometer-spectrometer 

spec i 4 i cat i ons . 


Throughput  . 

Scan  period  (minimum)  .  . 
Collector  diameter  .  .  .  . 

Field  o4  view  (4ull  angle) 

Spectral  range  . 

Spectral  resolution  .  .  . 

Detector  type  . 

Detector  NEP  . 

Dynamic  ...... 

System  sensitivity  .  .  .  . 


.  .  . . .  0.28  cm^sr 

.  1  scan/  30  seconds 

.  . . 30.8  cm 

.  0.8  degrees 

.  .  .  .  .  .  .  .  .  0.8  to  1.6  Mm 

. 2  cm  ^ 

......  Intrinsic  germanium 

1.1  X  10  watts/4hz  at  1.5  Mm 

. . . 80  dB 

.  16  R/cm  ^  at  1.3  Mm 


Figurs  2-9.  Intsr'f •romctvr  system  equipped  with  20-inch 
diameter  telescope.  The  isocon  camera  is  also  mounted  on 
the  telescope. 
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CHAPTER  III 
MEASUREMENT  THEORY 

The  Inter-f erogram 


I ' 


The  electronic  signal  -from  the  i nter-f erometer  detector 


is  a  low--f  requency  time-continuous  signal  called  an 


i  nterf  erogram.  The  i  nter-f  erogram  is  a  scaled  analog  o-f  the 
incoming  light  -frequency.  Using  the  approach  o-f  Loewenstein 
C19711,  and  re-ferring  to  the  layout  o-f  a  Michelson 
i  nter-f  erometer  in  Figure  1-1>  the  i  nter-f  erogram  signal  is  o-f 
the  form 

f  »  i/'ff  ,  (3.1) 

s 

where 

9  s  input  light  wavenumber  in  cm 

{/'  =  rate  of  change  of  effective  path  difference  in 
cm/second , 

f^*  scaled  frequency  output  of  interferometer  (the 
i nterf erogram)  in  Hz, 

The  detector  used  in  this  study  is  sensitive  in  the  region 
from  12500  to  6250  cm  ^  (X  «  0.8  to  1.6  )im)  .  The  scan  rate 

for  the  interferometer  was  set  at  30  secondsi  therefore,  the 
input  light  frequency  was  scaled  to  audio  frequencies  of 


less  than  200  Hz 


Th*  lisht  poMsr  reaching  the  detector  in  an  ideal 
i nter-f erometer  is  CLoewenstein  19713 

P.  .  =  2A^<1+cob  2«dx>  ,  (3.2) 

det 

where 

A  *  amplitude  o-f  on-axis  monochromatic  point  source, 

X  s  e-f-fective  path  di-fference  in  cm. 

In  the  Loewenstein  model,  the  single  line  04  Equation  3.2  is 
replaced  by  an  input  power  density  spectrum  of  the  ‘form 
A^sB(d).  Making  this  substitution  in  Equation  3.2, 

neglecting  the  constant  (dc)  part,  and  integrating  over  ff 
produces 

Z(x)  >  2\Bt9}  coB(2*dx)  efV  ,  (3.3) 

*  • 

where  I(x)  is  defined  as  the  i nterf erogram.  The  desired 
information  B<9)  (the  input  spectrum)  is  obtained  by  taking 
the  inverse  Fourier  transform  of  I(x).  The  signal  produced 
by  the  interferometer  is  referred  to  as  a  "double-sided"  or 
"symmetric"  i nterf erogram.  An  example  of  a  double-sided 
i nterf erogram  is  shown  in  Figure  3-1 (  note  that  the  same 
information  is  available  on  both  sides  of  the  large  center 
which  occurs  at  zero  path  difference.  The  interferometer 
slide  is  moved  an  equal  distance  on  each  side  of  the  optical 
zero  path  difference  to  create  a  symmetrical  signal. 

Since  the  i nterf erogram  I(x)  in  Equation  3.3  is 

symmetric,  the  spectrum  Bid)  may  be  obtained  using  a  simple 
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Figure  3-1  Typicel  “double  sided“  i nteri erogram. 

cosine  trans-form  o-f  the  -form 

m 

6(cr)  »  J  I(x)  cos(2ir<rx)  <fx  .  (3.4) 

The  power  spectral  density  is  computed  with  a  digital 
computer  using  a  Fast  Fourier  Trans-form  (FFT)  .  The  FFT  used 
to  process  the  data  -for  this  study  was  developed  by  Ware 
C19801.  The  digitized  i nter-f erogram  was  produced  using  a  12 
bit  analog  to  digital  converter  on  the  detector  amplifier 
output  and  enough  samples  were  taken  to  perform  a  16,384- 
point  transform.  The  Implementation  of  the  FFT  is  beyond 
the  scope  of  this  paper  but  that  used  in  Fourier  transform 
spectroscopy  was  adapted  by  Forman  [19661  from  the  radar 
signal  processing  work  of  Thomas  Stockham,  then  at  MIT 
Lincoln  Laboratory.  The  FFT  is  explained  extensively  by 
Brigham  [19741. 

The  computation  of  the  FFT  requires  that  the 
i nterf erogram  be  sampled  at  uniform  increments  of  path 
difference.  Ideally,  the  laser  reference  channel  should  be 
a  short  wavelength  laser  beam  which  passes  through  the 
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interferometer  optical  train  as  does  the  signal.  Otherwise, 
the  reference  beam  may  experience  different  motion  stability 
than  the  signal  beam,  due  to  nonidentical  geometry. 
However,  for  simplicity,  a  separate  laser  (X  =  6328  A> 
interferometer  channel  was  used  and  provides  the  uniform 
sampling  function  signal.  The  laser  channel  signal  is 
divided  to  obtain  a  sample  rate  of  650  Hz.  The 
interf erogram  is  digitized  at  a  free  running  rate  of  50  kHz. 
The  values  taken  between  laser  channel  zero  crossings  are 
then  averaged  and  stored  as  the  data  point  for  the  slide 
position  corresponding  halfway  between  adjacent  laser 
crossings.  This  over-sampling  method  will,  according  to 
Ware  C 19803,  minimize  the  system  noise  gain  as  well  as 
minimize  the  effects  of  slide  velocity  variations  in  the 
sampled  data. 


Calibration  Source 

The  determination  of  relative  OH  spectral  line  strengths 
is  sufficient  to  compute  rotational  temperatures.  Since 
this  is  a  major  interest  in  this  study,  a  calibration 
technique  to  produce  a  relative  instrument  response  was 
developed.  The  absolute  calibration  of  an  interferometer- 
spectrometer  is  an  intricate  process!  details  of  the 
approach  are  outlined  by  Wyatt  C 19781.  Ware  119801 
performed  an  interferometer  wavenumber  response  calibration 
as  a  function  of  optical  alignment.  He  found  that  the 


relative  response  must  be  reestablished  each  time  the 
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instrument  is  realigned  at  the  time  of  data  collection 
because  the  instrument  would  only  remain  in  acceptable 
alignment  for  about  2  hours.. 

A  relative  instrument  response  was  obtained  by  causing 
the  interferometer  to  view  a  calibration  source  composed  of 
a  tungsten  bulb  illuminating  a  pair  of  ground  glass 
diffusion  screens  (see  Figure  3-2).  The  intensity  of  the 
illumination  was  controlled  by  placing  a  plate  with  a  small 
3-mm  aperture  between  the  bulb  and  the  diffusion  screens. 
The  second  screen  was  then  imaged  by  a  projector  lens  onto  a 
24-inch  diameter  ground  glass  viewing  screen  placed  ten  feet 
away.  The  net  brightness  of  the  source  was  adjusted  so  the 
calibration  source  intensity  as  viewed  by  the  interferometer 
was  approximately  that  of  the  night  sky.  The  calibration 
source  assembly  was  mounted  inside  of  a  24- inch  diameter, 
10-foot  long  tube.  The  viewing  screen  was  sufficiently 
large  to  completely  fill  the  interferometer  field  of  view, 
thus  providing  a  diffuse  source. 

The  long  length  (3.0  meters)  of  the  calibration  source 
was  chosen  to  ensure  that  the  illumination  on  the  viewing 
screen  was  uniform  across  the  entire  surface  viewed  by  the 
interferometer.  The  quality  of  the  calibration  source  was 
calculated  using  the  geometry  described  above  and  assuming 
the  tungsten  bulb  was  a  blackbody  source  with  correction  for 
the  emissivity  of  tungsten. 


VIEWING  SCREEN 


Figure  3-2.  I nter-f erometer  Calibration  Source. 

Instrument  Response 

A  relative  instrument  response  was  obtained  each  time 
the  instrument  was  realigned  during  the  data  taking  process. 
The  calibration  alignment  process  was  repeated  about  every  2 
hours.  The  i nter-f erograms  obtained  from  the  instrument  when 
observing  the  calibration  source  were  transformed  using  the 
same  FFT  and  apodization  (discussed  later  in  this  chapter) 
routines  as  were  used  for  the  airglow  data  signal  processing 
described  above.  The  calibration  source  spectral  or 
"blackbody"  curves  were  then  averaged  together  (from  S  to  10 
frames)  to  minimize  noise  and  irregularities. 

The  tungsten  bulb  in  the  calibration  source  has  a 
blackbody  equivalent  temperature  of  2370  *K  when  operated  at 
the  specified  current  (750  mA) ,  according  to  Gilway 
Technical  Lamp  C 19821,  its  manufacturer.  The  Plank  equation 
can  be  used  to  calculate  the  spectral  sterance  L^(X)  in 
watts/meter,  for  the  assumed  blackbody  radiation  source  as, 


shown  by  Wyatt  C 19781 


Lb<X) 


he 

X;^T 


) 


> 


(3.5) 


t (exp  - 


TTl 


where 

h  =  6.6262  X  10  (Js)  (Plank’s  constant)  (Wyatt  19783, 

c  =  2.9979  X  10®  (m/s)  (speed  o-f  light)  (Wyatt  19783, 

X  »  wavelength  (m) , 

-23 

k  -  1.3806  X  10  (J/*K)  (Boltzmann’s  constant) 

(Wyatt  19783, 

T  =  absolute  temperature  (*K) . 

The  i nter-f erometer  detector  measures  energy  coming  from 
the  interference  of  two  light  rays  of  the  same  optical 
frequency;  and  the  germanium  detector  operates  in  a  photon 
sensitive  mode.  As  a  consequence,  in  the  analysis  it  is 
desirable  to  manipulate  Plank’s  equation  into  terms  of 
wavenumber  and  quanta  (photons).  The  relationship  between 
sterance  and  photon  sterance  is 

L  =  L„X/(Ac)  .  (3.6) 

p  B 

Noting  the  relationship  between  wavenumber  and  wavelength 
and  that  Equation  3.5  is  a  density  function 

<r  =  IVX  -,  ■  •  •  (3-.  7) 

dff  «  -  1/X^  erx  .  (3.8) 


Using  Equations  3.7 

centimeters.  Plank 

-2  -1 
sr 


and  3.8  and  converting  from  meters  to 
B  equation  can  now  be  shown  in  terms  of 


photons  cm 


sec 


or 
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zero  path  difference  point  is  not  centered  in  the  signal 
time  window. 

The  amplitude  of  the  measured  spectrum  could  be  obtained 
by  the  magnitude  operation  (square  root  of  the  sum  of  the 
squares  of  the  real  and  imaginary  parts).  However,  this 
operation  always  yields  positive  noise  components, 
increasing  the  noise  by  f2.  This  increase  in  noise  is 
especially  harmful  when  several  frames  of  data  are  signal 
averaged  (coadded)  to  help  identify  low  weak  features  in  the 
airglow  emission  spectrum.  The  sampling  phase  shift  error 
may  also  vary  from  frame  to  frame.  Consequently,  a  method 
of  phase  correction  based  upon  the  data  within  each  frame 
must  be  used. 

The  phase  characteristics  of  each  data  frame  could  be 
obtained  by  Fourier  transforming  a  small  data  set  around  the 
''center”  of  the  i nter f erogram.  This  process  would  yield  a 
very  low  resolution  spectrum  from  which  the  slowly  varying 
phase  information  could  be  easily  extracted.  However,  this 
truncation  and  transform  process  is  the  same  as  convolving 
the  original  spectrum  with  the  Fourier  transform  of  the 
truncating  function.  As  described  by  Ware  tl9801,  this 
convolution  is  merely  a  digital  filter  operating  in  the 
frequency /phase  domain.  Hamming  C 19771  showed  that  for  a 
given  filter  width  the  minimum  noise  gain  is  obtained  when 
the  convolving  function  is  rectangular  in  shape.  The  phase 
characteristics  are  not  related  to  the  incoming  signal  in 
any  way  so  the  rectangular  filtering  process  may  be 
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repeated  as  many  times  as  necessary  to  identi-fy  and  remove 
the  phase  anomalies  -from  the  computed  spectrum.  When  the 
■filtering  process  is  complete,  the  remaining  "real"  part  o-f 
the  trans-form  is  used  as  the  measured  spectrum  (power 
spectral  density  -function). 

The  spectral  -frames  were  phase  corrected  with  convolving 
-filter  windows  o-f  3,  5,  9,  17,  and  65  sample  point  widths  in 
a  repetitive  manner.  Figure  3-3  shows  an  instrument 
uncorrected  blackbody  curve  both  be-fore  and  a-fter  phase 
correction,  and  Figure  3-4  shows  the  same  in-formation  -for  a 
spectral  -frame.  Examination  o-f  Figures  3-3  and  3-4  show 
that  the  imaginary  portion  o-f  the  power  spectral  density 
(caused  by  chromatic  variations  within  the  instrument 
optics)  has  been  eliminated  by  the  phase  correction 
algorithm. 

Apodization  and  Interpolation 

The  calculation  o-f  OH  airglow  rotational  temperatures 
was  a  major  goal  o-f  this  study.  The  calculation  o-f 
temperature  requires  the  extraction  o-f  relative  spectral 
line  intensities  -from  the  trans-form  o-f  the  i nter-f erogram. 
The  relative  instrument  response  and  phase  correction  steps 
described  above  are  modifiers  in  signal  processing  the  line 
intensi t ies. 

The  shape  of  an  individual  spectral  line  is  of  interest 
because  the  shape  will  affect  how  the  a  best  estimate  of  the 
line  emission  intensity  is  extracted  from  the  data.  Due  to 
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th«  discrete  nature  of  the  transitions  of  the  excited  OH 
molecule,  the  spectral  lines  entering  the  interferometer  are 
discrete  in  shape  except  for  a  small  finite  width  due  to 
Doppler  and  collisional  broadening  plus  "time  windoMing"  of 
each  photon  source.  The  numerical  computation  of  an  FFT 
requires  that  the  sampling  be  limited  in  length.  The 
truncation  or  multiplication  of  the  i nterf erogram  by  a 
rectangular  ‘‘windoM*  in  the  time  domain  is  the  same  as 
frequency  domain  convolution  of  the  discrete  spectral 
frequencies  with  the  Fourier  transform  of  the  rectangular 
window  function.  The  Fourier  transform  of  the  rectangular 
window  is  a  <sin  x>/x  or  "sine*  function.  This 
characteristic  spectral  line  shape  is  referred  to  as  an 
"instrument  function."  The  high  side  lobe  behavior  of  the 
sine  function  (largest  side  lobe  ->13  dB  down  in  amplitude) 
causes  the  various  lines  in  the  spectrum  to  interact  or  mesh 
together.  Therefore,  it  is  desirable  to  establish  a 
technique  to  suppress  the  sine  function  behavior  and  force 
the  signal  component  into  as  much  of  a  discrete  line  shape 
as  possible. 

The  process  of  suppressing  the  sine  function  side  lobes 
is  called  apodization.  The  apodization  of  the  spectral 
information  can  be  accomplfshed  either  in  the  i nterf erogram 
domain  by  multiplying  the  i nterf erogram  with  another  window 
then  performing  the  FFT,  or  by  convolving  the  transform  of 
the  window  with  the  transformed  i nterf erogram .  The 


convolutian  o-f  the  spectral  data  with  the  transform  of  the 
chosen  window  is  the  method  used  for  this  study. 

Numerous  apod izat ion  functions  have  been  tried  for  use 
with  Fourier  transform  spectroscopy.  In  choosing  an 
apodizing  function  the  trade-off  is  among  instrument 
function  width  (resolution),  side  lobe  attenuation,  and 
computational  efficiency.  Norton  and  Beer  [19763  computed 
over  IlOO  different  apodization  windows  and  plotted  each  as 
a  function  of  central  peak  width  versus  height  of  maximum 
side  lobe.  Their  study  indicated  an  optimal  boundary 
existed  between  the  two  plotted  parameters.  Figure  3-5 
shows  the  Norton  and  Beers  limit  with  some  specific 
functions  also  shown.  V/agin  C 19801  independently  computed 
over  3000  functions  and  also  showed  that  this  boundary 
exists.  However,  he  went  on  to  analytically  show  that  for  a 
chosen  characteristic  (either  side  lobe  suppression  or 
resolution)  an  optimal  apodization  function  can  be  computed. 
Harris  [19783  presents  an  excellent  analysis  of  many  of 
these  apodization  functions  and  shows  the  trade-offs 
associated  with  each.  Nuttall  [19813  shows  some  corrections 
to  Harris’s  work  and  presents  additional  functions  for 
consideration. 

The  present  interferometer’s  full-width  half-maximum 
(FWHM)  instrument  function  is  1.8  cm  The  OH  spectrum  in 
the  l-Mm  wavelength  region  has  line  separations  on  the  order 
of  lO  cm  This  occurrence  renders  the  prime  consideration 
in  the  choice  of  an  apodization  function  as  one  of  side  lobe 


attenuation  and  ease  o-f  computation, 


rather  than  maximum 


resolution . 

As  can  be  seen  -from  Figure  3-5  the  Hamming  window 
(sometimes  referred  to  as  the  “minimum  2-point"  window)  is 
located  on  the  optimal  boundary.  This  -function  is 

comparatively  simple  to  compute  because  it  contains  only  two 

terms  and  provides  side  lobe  attenuation  of  -43  dB  (a  20  dB 
improvement  over  that  provided  by  the  sine  function).  The 
additional  side  lobe  attenuation  is  sufficient  because  the 
data  collected  for  this  study  have  s i gnal - to-no i se  ratios  of 
about  lOO,  thus  placing  the  side  lobe  behavior  below  the 
noise  level.  The  Hamming  window  yields  good  spectral 

resolution  by  providing  a  FWHM  central  lobe  of  2.7  cm  ^ 
which  is  sufficient  to  identify  the  OH  spectral 

characteristics.  This  window  also  has  an  asymptotic  side 
lobe  roll-off  of  6  dB  per  octave  [Harris  1978J. 

The  Hamming  convolution  function  takes  the  form 

k*-«-6  r 

I,  =  S  I,  0.53836  sincCwc (k6d« ) 3 
k=-6 

*  0.46164  rsinctwCc(k6(y>  )-«-13>-*-sinc<TCc(k6ff>  )-13>1  | 

2  ■* 


<3. 10) 
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k  *  counter  for  integer  sample  numbers  centered  at 
transform  point  nearest  9t , 


»  wavcnumbvr  at  intagar  aampla  in  cm  , 

>  wavanumbar  at  point  whara  intanaity  ia  baing 
computad  in  cm 

>  intanaity  o-f  trans-formed  data  at  intagar  trana-form 
wavanumber  <T , 

»  intensity  o-f  apodizad  data  at  wavanumnbar  dt  > 

»  number  o4  sample  points  par  wavanumber. 


)(<-TRIANGLE 
<-KAISER  BESSEL 


HAHHING-> 


NORTON  AND  BEER'S  LIMIT-?^ 


H  10-5 


NOHmiLlZtO  HHW  WIDTH  CW/Hj] 

Pigura  3-S.  Optimal  apodization  -functions, 
resolution  vs.  side  lobe  attenuation  CEspy  198<l]. 


showi ng 


The  discrete  spectral  data  points  computad  by  the  FFT 
occurred  at  l.S-cm  ^  intervalsi  there-fore,  as  can  be  seen 
-from  Equation  3.10,  the  convolution  computation  was  summed 
over  ±  4  sample  numbers.  Examination  o-f  Equation  3.10  shows 
that  this  equation,  in  addition  to  suppressing  side  lobes, 
can  be  used  to  interpolate  spectral  values  between  the 


each  i nterf erogram  was  also  stored  on  digital  tape  -for  later 


processing. 

The  information  of  interest  within  each  spectral  frame, 
needed  for  rotational  temperature  calculations,  was  the 
intensity  of  a  few  specific  OH  spectral  lines  <the 
temperature  model  is  explained  in  Chapter  IV) .  In  order  to 
identify  the  spectral  features  of  interest,  each  of  the 
transformed  frames  of  data  was  plotted  on  paper  using  an  x>y 
recorder  and  a  simple  computer  routine  which  performed  a 
straight  line  connect  between  each  of  the  discrete  points 
computed  by  the  FFT.  Each  of  the  plotted  frames  was 
examined  manually,  and  the  lines  of  interest  were 
identified  by  FFT  sample  number.  The  lines  selected  were 
bright  lines  within  each  OH  band  which  were  far  enough  apart 
as  to  not  be  contaminated  by  other  features  and  situated  at 
wavelengths  of  minimal  atmospheric  absorption  by  *^2^  and 
CO^.  Examination  of  the  plots  showed  that  the  position  of 
each  line  never  varied  from  frame  to  frame  more  than  ±1  FFT 
sample  number.  The  selected  FFT  sample  point  was  then 
recorded  as  the  initial  position  of  each  of  the  spectral 
lines  of  interest. 

The  extraction  of  spectral  line  intensities  was 
accomplished  using  the  Hamming  apod i ration  function  to 
interpolate  between  the  discrete  FFT  data  points.  Allowing 
for  additional  error  in  the  chosen  position  of  each  line, 
the  interpolation  routine  was  operated  in  a  sample  number 


window  13  data  points  around  that  manually  chosen  line 


position.  The  line  intensity  was  then  computed  at  0.01 
sample  number  increments  (O.OIS  cm  ^  increments)  over  the 
data  point  window  usin9  the  Hammin9  -function.  The  routine 
saved  the  maximum  value  -found  within  the  sample  window  as 
the  line  amplitude  -for  that  spectral  feature. 

The  instrument  calibration  -function  and  the  correct 
spectral  line  amplitudes  could  now  be  used  in  the 
temperature  model  for  extraction  of  OH  rotational 


temperatures  as  a  function  of  each  Interferometer  scan. 


CHAPTER  IV 


OH  ROTATIONAL  TEMPERATURE  MODELING 

Introduct ion 

Hydroxyl  is  a  minor  atmospheric  constituent,  residins 

in  a  thin  layer  (about  7  km  thick)  at  an  altitude  near  87  km 

IBaker  et  al .  19831.  Although  minor  in  concentration 

4  6  3 

(between  10  and  10  molecules/cm  at  night,  Baker  119783) 
it  is  the  major  atmospheric  near- i nf rared  airglow  radiator 
at  night.  The  radiation  generated  by  OH  occurs  in  spectral 
bands  known  as  the  Meinel  bands,  since  their  discovery  and 
identification  by  Meinel  119503.  After  their  discovery,  the 
measurement  of  these  bands  has  been  of  great  interest  to  the 
atmospheric  science  community.  OH  airglow  emission  band 
measurements  contain  information  about  mesospheric 
populations  and  temperatures,  and  as  such  provide  insight 
into  the  dynamics  of  the  entire  middle  atmospheric  region. 
A  rotat ional -temperature  model  is  developed  in  this  chapter 
similar  to  the  technique  used  by  Hill  et  al .  C19793.  The 

field  measurements  are  then  fit  to  the  model  in  a  least- 
squares  sense.  The  quality  of  the  fit  will  also  determine 
error  bounds  on  both  relative  band  intensity  and  absolute 


rotational  temperature 


The  Excited  OH  Radical 


The  hydroxyl  radical  is  a  diatomic  heteronuc 1  ear 
molecule  and  has  a  nonzero  electric  dipole  moment.  The 
molecule,  therefore,  has  the  ability  to  be  easily  excited 


and  radiate  electromagnetic  radiation.  The  OH  Meinel 
radiation  bands  in  the  infrared  are  caused  by  vibration- 
rotation  transitions  within  the  ground  electronic  state  of 
the  OH  molecule.  (The  electronic  absorption  and  emission 
spectra  of  OH  occur  in  the  ultraviolet.) 

Hydration  of  ozone  and  perhydroxyl  reduction  are  the 
primary  processes  for  the  creation  of  v ibrat i onal ly-exc i ted 
OH  in  the  earth’s  upper  atmosphere  [Baker  19781.  The 


vibrational ly-excited  states  of  OH  are  quantized  and  are 
populated  according  to  certain  dipole  selection  rules 
described  by  Hertzberg  [19713.  The  quantum  numbers 
associated  with  the  v ibrat i onal 1 y-exc i ted  states  take  on 


values  V  »  O, 1 , 2 , 3 , . . . , 9  and  the  allowed  transitional 
changes  during  emission  are 


Av  —  1,2,3,  ...,9 


(4.1) 


where  v’  is  the  upper-state  energy  level  and  v*  ’  is  the 
lower-state  energy  level.  The  Meinel  radiation  bands  are 
known  by  the  value  of  6v  associated  with  the  transition. 
The  band  sequences  measured  in  this  study  are  Avs2  (4,2  and 
3,1)  and  Av»3  (8,3  and  7,4)  transitions  since  they  have 
energy  differences  in  the  near- i nf rared  region.  A  typical 
measured  spectrum  is  shown  in  Figure  4-1. 


transitions  in  molecular  rotational  ener9>'  in  addition  to 


the  changes  in  vibrational  levels.  These  angular  momentum 
changes  modulate  the  vibrational  transitions  causing  the 
intricate  spectral  structure.  Appendix  A  includes 
additional  in-formation  on  OH  radical  transition  theory. 


Rotational  Temperature  Model 


In  extracting  temperatures  from  OH  spectral  data  it  is 
assumed  that  the  OH  molecular  population  is  in  rotational - 


state  thermal  equilibrium.  The  population  distribution  over 


the  di-f-ferent  quantum  numbers  can  then  be  described  by  the 


N  .  =  number  o-f  molecules  in  vibrational  level  v*  and 

v'  J* 

rotational  level  J*  at  temperature  T, 

»  2J’ -H  or  the  statistical  weight  o-f  state  J’  with 

its  (2J' ■•’1 ) --fold  degeneracy,. 

C  .  >  a  constant,  equal  to  N  . /Q  where  Q  is  the 

v»  ’  ^  v*  r  r 

partition  function  for  rotational  level  v, 

E^,j,  =•  energy  of  the  rotational  state,  and  equal  to 

F  (J* )Mc,  where  F  . (J’ )  is  the  term  value  for 
v’  v’ 

the  upper-rotational  state, 
ft  «  Boltzmann’s  constant. 

The  volume  emission  rate  in  photons  sec  ^  cm  ^  of  a 
central  line  is  given  by  Mies  C19743  as 


^v»J»,v»»J»»  “  *^v*J»  ^v’J’,v'’J''  ’  (4.3) 

where 

^  I  T«  t f -n  f  “  intensity  of  emission  from  upper-state 

N/  w  f  V  J 

v’J*  to  lower-state  v’’J*’, 

^vtji  *  Einstein  transition  probability  of 

spontaneous  emission  from  state  v'J*  to 
state  v’ ’ J»  » . 

Inserting  Equation  4.2  into  Equation  4.3  gives  the  volume 
emission  rate  of  each  rotational  line  arising  from  a 


transition  from  an  upper-state  v'J’  to  a  lower-state  v''J'' 


where 


»  the  value  o-f  the  model  at  wavenumber  i, 

A  »  relative  total  integrated  band  intensity, 

0 

3 

K  A(2J'-t’l)  »  <7  S,  which  is  the  Einstein  coe-f  t  i  c  i  en  t 

times  the  J*  state  degeneracy  or  the  wavenumber 

cubed  times  the  line  strength, 

B  =  hctkJy  where  T  is  rotational  temperature  in  ‘K, 

*  term  value  -for  upper-state  transition, 

=  relative  instrument  response  at  wavenumber  i, 

C  >  summation  over  all  lines  within  the  band,  is  a 

scaling  factor  to  ensure  the  entire  band  intensity 

is  contained  within  the  A  term. 

The  model  shown  above  is  similar  to  the  intensity 

calculation  shown  in  Equation  4.5.  The  difference  is  that 

the  model  pertains  to  an  entire  band  of  transitions  and  is 

modified  by  to  appear  like  the  measured  data.  All 

quantities  within  the  model  are  now  known  except  the  band 

intensity  A  ,  and  the  temperature  term  B  .  Within  each  OH 
0  0 

band,  the  measured  data  points  can  be  subtracted  from  the 
model  at  the  same  wavenumber  and  any  differences  will  be 

contained  within  small  changes  of  A  and  B  times  the 

•  0 

derivative  of  the  model  at  points  A  and  B  . 

0  0 


r  3Mj<A^,B^)  ^ 

X  A  B  1 

3M. (A  ,B  )  1 

1  0  0 

L  3  A  J 

▼  o  o  1 

0  ^ 

3B  J 

9 

•  0 


where  is  the  measured  data  point.  Equation  4.6  can  be 
shown  for  (J )  data  points  within  a  band  in  matrix  form  as 


’  r 


where  the  le+t  matrix  is  the  di  +  -ference  between  each 


measured  line  intensity  within  a  band  and  the  line  intensity 
as  predicted  by  the  model*  The  right  hand  side  of  Equation 
4.7  must  be  equal  to  the  preceding  difference. 

The  computation  of  Equation  4.7  requires  the  solution 
for  each  of  the  derivatives 


and , 


dB 


Rj  Ij  exp(-F^  ) 
3A~  “  N 


(4.8) 


exp(-F  B  ) 

t  9 


N 


[  t  = 

j 


(  I  F  exp< -F  B  > )  3  - 
j  j  J  • 


N  F 


J 


(4.9) 


where  N  »  £  t  exp(-FjB  )  or  the  normalization  factor.  All 
J  j  • 

the  terms  in  Equation  4.7  are  computed  and  submitted  to  a 
least-squares  fitting  routine,  developed  by  Lawson  and 

Hanson  C 19741  which  returns  the  values  of  AA  and  AB  .  The 

0  • 


least-squares  routine  is  used  to  sol*  e  problems  of  the  form 


where 

A  a  the  derivative  matrix  o-f  Equation  4.7, 

X  s  the  unknown  or  4  matrix, 

B  a  the  di-f-Ference  matrix. 

To  besin  the  process  a  starting  temperature  was  chosen 
and  all  terms  computed.  The  result  of  this  first  -fit  is  to 
determine  an  appropriate  initial  value  for  A  (A  is  merely 

•  t 

a  relative  intensity  term).  Once  the  initial  values  for  B 

• 

and  A^  were  ascertained,  the  least-square  fitting  routine 
was  operated  in  a  loop,  the  derivative  matrix  vectors  were 
each  scaled  to  unity  value  before  submission  to  insure 
convergence  of  the  fit,  and  each  time  the  routine  would 

return  values  for  4A  and  aB  .  The  4  values  were  added  to 

0  • 

the  previous  values  for  A  and  B  ,  the  model  values 

•  • 

recomputed  and  resubmitted  to  the  routine.  The  value  of  4B 

I 

was  tested  after  each  iteration  until  the  temperature  change 
(temperature  is  part  of  )  from  the  previous  fit  was  less 
than  0.9*  K.  When  the  least-squares  routine  converged  to 
the  0.9  *K  limit  the  values  for  A  (relative  band  intensity) 
and  T  (the  temperature  part  of  B^ )  were  stored  along  with 
the  time  and  date  of  the  frame. 

The  constant  terms  used  in  the  model  were  computed  by 
Espy  C 19841  using  the  molecular  constants  measured  bv  Coxon 
and  Foster  [19821  and  OH  dipole  moments  derived  by  Werner  et 
al.  C 19831.  The  values  of  these  terms  are  given  in  Tables 
4-1  through  4-4.  The  tables  show  the  wavenumber,  a  typical 
relative  Instrument  response,  the  term  value  for  the  upper- 


TABLE  4-1.  Molecular  data  tor  OH  Avs2,  <4, 2)  band. 


Line 

Wavenumber 
<cm-> ) 

Relative 

Response 

Term  Value 
(cm-‘ ) 

Line  Strength 

Pi 

(4)-N- 

6159.634 

O.  147 

13421.921 

0.23415100E4-13 

Pt 

(A)* 

6174.994 

O.  166 

13514. 125 

O.  19663726E-I-13 

Pi 

(3)* 

6200.353 

0.  199 

13320. 758 

0.  16226162E-I-13 

Pt 

(3>* 

6219.070 

0.224 

13428.966 

0. 13226441E+13 

Pi 

(2)« 

6238.114 

0.247 

13248.923 

0.88080260E-I-12 

Pt 

(2)« 

6260.766 

0.275 

13377. 615 

0. 70291000E-t-12 

Qi 

<3) 

6301.517 

0.326 

13421 . 922 

O. 25067365E+12 

Qi 

<2) 

6309.948 

0.336 

13320.758 

0. 71495030E+12 

Qi 

(1) 

6315.828 

0.344 

13248.923 

0.  11681950E-f  13 

Ri 

( 1 ) 

6368.381 

0.411 

13428.966 

0.  60586570E-I-12 

Ri 

(1) 

6387.661 

0.440 

13320. 758 

0.74062960E-H12 

Ri 

<2) 

6411. 112 

0.473 

13421.921 

0.  12600000E-I-13 

<* 

lines  used  in  tit) 

• 

TABLE  4 

-2.  Molecular 

data  tor 

OH  Av»2,  (3, 

1)  band. 

Line 

Wavenumber 

Relative 

Term  Value 

Line  Strength 

<cm-> ) 

Response 

(cm-* ) 

Pi  <4)* 

6480.230 

0.572 

10352.446 

0.27218330E-*-13 

Pi  (4)* 

6495.579 

0.593 

10443.293 

0.22819520E4-13 

Pi (3)* 

6522. 332 

0.626 

10247.071 

0. 18861154E+13 

Pt  (3)* 

6541.239 

0.650 

10354.206 

0. 15350O37E+13 

Pi  <2)* 

6561.409 

0.679 

10172.300 

0. 10240617E+13 

Pt  <2)* 

6584.562 

0.712 

10300.410 

0.81566820E-M2 

Qi  (3) 

6627.706 

0.770 

10352.446 

0.57940270E+12 

Qi  (2) 

6636. 171 

0.782 

10247.071 

0.82784900E-t-12 

Qi  (1) 

6642.081 

0.  789 

10172.300 

O.  13556844E4-13 

Rt  (1) 

6697.063 

0.867 

10354.206 

0. 35196955E+12 

Ri  (1) 

6716.853 

O.  892 

10247.071 

0.86348040E+12 

Ri  (2) 

6741.555 

0.931 

10352.446 

0. 14703511E+13 

(4'  lines  used  in  tit) 


TABLE  4 

-3.  Molecular 

data  -for 

OH  Av»3,  (8, 

5)  band. 

Line 

Wavenumber 

Relative 

Term  Value 

Line  Strength 

( cm-* ) 

Response 

<cm-‘ ) 

Pi  (5) 

75SS. 135 

1.000 

24163.096 

0.34607000E+13 

Pi  <S) 

7569.641 

0.997 

24248.912 

0.47194710E-H3 

Pi  (4) 

7600.862 

0.992 

24055.028 

0.41913920E-I-13 

Pf  (4) 

7617.468 

0.990 

24152. 797 

0.35340550E-i-13 

Pi (3)* 

7642.343 

0.984 

23971.237 

0.29201550E-I-13 

Pi  <3)* 

7661.252 

0.985 

24083.736 

0.23944700E-H3 

Pi  (2)* 

7679.631 

0.978 

23911.582 

0.  15952644E4-13 

Pi  <2)* 

7700.808 

0.972 

24042. 113 

0.  12846050E-M3 

Qi  (3) 

7726. 136 

0.960 

24055.028 

O.  95265240E-M2 

Qt  (2) 

7739.286 

0.952 

23971 . 237 

0.  13444194E-M3 

Qt  ( 1 ) 

7748.311 

0.945 

23911.582 

O.  21738O30E-*-13 

Ri  (1> 

7791.408 

0.900 

24803. 736 

O.  11463823E-I-13 

Ri  (1) 

7808. 166 

0.874 

23971.237 

0.  13900488E-M3 

Ri  (2) 

7823.076 

0.837 

24055.028 

0. 23899920E-H3 

('X'  lines  used  in  -fit) 

TABLE  4 

-4.  Molecular 

'  data  -for 

Av=3,  (7,4) 

band . 

Line 

Wavenumber 

Relative 

Term  Value 

Line  Strength 

<cm-* ) 

Response 

(cm-* ) 

Pi 

(5)* 

8049.708 

0.786 

21763.099 

0.66072740E-I-13 

Pi 

(5)-» 

8064.112 

0.779 

21847.340 

0.36979800E-H3 

Pi 

(4)* 

8096.397 

0.763 

21649. 131 

0.40669700E-I-13 

Pi 

(4)» 

8113. 112 

0.734 

21743.319 

O.426633O0E-M3 

Pi 

(3)* 

8138.903 

0.740 

21560.824 

0.35302450E-M3 

Pi 

(3>* 

8158.204 

0.723 

21672.328 

0.28898620E-t-13 

Pi 

(2)* 

8177.237 

0.711 

21497. 994 

0.  19277462E-H3 

Pi 

(2)« 

8199.208 

0.696 

21628.204 

0.  15491604E-M3 

Qt 

(3) 

8227.210 

0.680 

21649. 131 

O.  11407934E-H3 

Qt 

(2)X- 

8240.067 

0.672 

21560.824 

0.  16141494E-M3 

Qi 

( 1 )« 

8249.071 

0.668 

21497.994 

O.  26171900E-I-13 

Ri 

( 1 

8294.714 

0.646 

21672. 328 

0.  13787860E-M3 

Ri 

( 1 )-» 

8311.901 

0.642 

21560.824 

0.  16761780E-*-13 

Ri 

(2)* 

8328.374 

0.635 

21649. 131 

0.  28788900E-I-13 

(* 

lines  used  in 

-fit) 

state  transition,  and  the  theoretical  line  strensth  -for  each 
line.  The  wavenumber  shown  is  the  average  o-f  the  two  A- 
doubled  lines;  the  term  value  is  also  the  average  o-f  the  two 
lines,  and  the  line  strength  given  is  the  sum  of  the  two 
line  strengths.  The  sum  of  the  line  strengths  was  used 
because  the  measured  line  is  actually  the  sum  of  the  A- 
doubled  pair.  Data  for  all  the  strong  lines  are  listed  for 
use  in  the  normalization  factor  of  the  model  but  only  those 
with  were  used  in  the  fit.  The  reason  all  were  not  used 
in  the  fit  was  that  many  were  determined  to  have  been 
contaminated  by  water  vapor  absorption  or  by  other  emission 
lines  CRoychourdhury  19831. 


Error  Analysis  and  Testing  of  Model 


The  model  used  here  is  over-determined  with  only  two 
unknowns  and  at  least  four  equations  (in  (8,5)  band  fit, 
more  in  other  bands)  and  as  such,  the  additional  information 
can  be  used  to  estimate  the  accuracy  of  the  least-squares 
fit.  Examination  of  Equations  4.7  and  4.10  show  that  when 
the  fitting  routine  is  complete,  the  derivative  matrix  A  is 
left  as  a  2X2  upper-triangulated  matrix.  The  rest  of  the 
terms  in  the  A  matrix  are  of  little  significance  because  the 
least-squares  routine  has  manipulated  the  data  into  the 
upper-triangulated  portion  in  order  to  solve  the  equations. 
The  upper-tr iangulated  matrix  now  takes  on  the  form 


(4.11) 


wher* 


R  =»  upp*r-tr iangulated  0064^  1  c i ents  returned  in  the 
derivative  matrix  A. 

Upon  completion  o4  the  4itting  routine,  the  values  le4t 
in  the  matrix  x  must  be  equal  to  zero,  except  for  noise  in 
the  fit,  and  the  upper  two  values  in  the  vector  B  must  also 
be  equal  to  zero,  also  except  for  noise,  because  the  system 
of  equations  has  been  solved.  But  because  noise  does  exist 
in  the  fit  then  some  residual  remains.  Therefore,  the 
equation  system  is  left  in  the  form 


R  X  »  n 


(4.12) 


where  n  is  the  noise  vector. 

In  order  to  calculate  the  accuracy  of  the  models’  fit  to 
the  two  system  unknowns  (temperature  and  band  intensity)  it 
is  necessary  to  obtain  the  covariance  matrix  which  contains 
the  variance  of  each  parameter  in  the  least-squares  fit.  It 
can  now  be  shown  by  manipulation  of  Equation  4.12  that  the 
covariance  matrix  is  CLawson  and  Hanson  19741 


<X  x^>  »  R  <n  n^>  R^ 


(4. 13) 
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wher* 

<x  x^>  =  covariance  matrix, 

R  »  inverse  o-f  R, 

T  indicates  the  transpose  operation, 

<n  n^>  =  the  variance  o-f  white  uncorrelated  noise. 

The  noise  in  the  system  -fit  is  assumed  to  be  white  and 

uncorrelated  because  all  the  mathematical  manipulations 

per-formed  were  linear  in  nature. 

A  measure  o-f  the  noise  in  the  system  is  returned  by  the 

least-squares  -fitting  routine.  The  returned  parameter  is 

the  Euclidean  norm  o-f  the  residual  vector  and  is  called 

R  .  The  variance  o-f  the  system  -fit  can  be  computed  from 
norm 

R  as 

norm 


=  (R  ^)  /  (M-K)  , 

norm 


<4. 14) 


where 

2 

0  »  variance  of  fit, 

W  s  number  of  equations  in  the  fit, 

K  >B  number  of  unknowns  in  the  fit. 

The  covariance  matrix  can  be  solved  in  terms  of  the  variance 
of  the  model’s  fit 


T  2  T 

<x  X  >  »  0  I  (R  R) 


(4.15) 


where  1  is  the  identity  matrix. 

The  standard  deviation  for  each  of  the  fit  variables  can 

be  computed  from  Equation  4.15  in  terms  of  the  original 

2 

variables  of  the  residual  matrix  R  and  the  variance 
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<r,  2 

2,2 


r  >2 

1,1  2,2 


(4.  16) 


where  9.  is  the  standard  deviation  O'f  the  intensity 
A 

function . 

The  standard  deviation  on  the  variable  B  (temperature  is 
part  o-f  B)  is 


=  — L  2 

B  2 

L  •'2,2 


(4. 17) 


Re-ferrln9  to  Equation  4.5  for  the  relationship  between 
temperature  T  and  the  variable  B,  the  standard  deviation  of 
the  temperature  can  be  calculated 


<r_  “  i  <B  ) 

T  B  0 


(4. 18) 


where  9^  is  the  standard  deviation  of  the  temperature  as 
calculated  by  the  model. 

The  model  presented  herein  provides  a  method  for 
calculating  relative  band  intensity,  rotational  temperature, 
and  the  standard  deviation  on  both  parameters.  Appendix  D 
presents  the  computer  programs  which  implement  the  data 
processing  algorithms  developed. 

As  verification  of  the  quality  of  the  model,  the  process 
outlined  was  used  to  compute  the  intensity  and  temperature 
of  synthetic  spectru  generated  by  Espy  [19841.  The  signal- 
to-noise  ratio  of  the  synthetic  spectra  was  varied  from  a 
low  of  2  to  a  high  of  20  and  then  each  was  submitted  to  the 
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much  higher  than  the  signal -to-noise  ratio  o+  the  test 
■frames}  there-fore,  the  uncertainty  should  continue  to 
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CHAPTER  \> 


RESULTS 

Intraduct  Ian 

Th*  h  i  9h-thrau9hput ,  narraw  -field  a-f  view 

inter-ferometer-spectrameter  was  taken  ta  several  air9law 
observation  sites  in  the  western  United  States  in  the  sprin9 
a-f  1983.  Over  several  months  of  observation,  one  of  the 
bri9htest  air9low  structure  events  was  seen  and  recorded  on 
Tune  IS,  1983,  thus  fulfillin9  the  primary  9oaI  of  this 
study.  The  observation  was  made  from  Sacramento  Peak,  New 
Mexico.  The  site  is  located  at  105*48*16”  west  lon9itude, 
32*47’57"  north  latitude  at  an  elevation  of  9370  feet.  The 
results  presented  in  this  chapter  were  taken  from  this  site 
between  June  13  and  June  15,  1983.  A  complete  catalo9  of 

the  results  computed  from  the  interferometer  data  is 
presented  in  Appendix  C. 


Back9round 

The  interferometer  with  an  i ma9e- i ntens i f ied  infrared 
isocon  camera  mounted  on  and  coali9ned  with  the 
interferometer  telescope  (see  Fi9ure  2-9)  was  used  to 
collect  data  at  the  observation  sites.  A  second  Isocon 
camera  was  mounted  on  its  own  tripod  and  was  used  to  search 
the  skies  for  indication  of  structure  prior  to  movin9  the 
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mor*  bulky  i  nter-f  •rotn*t*r .  Th*  in'frared  cameras  were 
supplied  by  the  University  o-f  Southampton,  Ensland.  Taylor 
£1983-841  operated  the  camera  equipment  and  provided 
experience  and  expertise  in  the  search  -for  airglow  structure 
events.  The  cameras  were  the  “eyes”  -for  the  i nter-f erometer 
in  determining  i-f  any  OH  airglow  structure  was  present. 
Once  an  area  o-f  OH  airglow  structure  was  located  in  the  sky 
the  earner a- i nterf erometer  system  was  positioned  to  view  that 
area.  The  data  -from  the  two  coaligned  systems  were  used  to 
correlate  the  viewed  and  calculated  intensities  and 
temperatures. 

An  in-frared  radiometer  CHuppi  19761  was  also  used  to 
gather  trend  and  absolute  intensity  data  o-f  the  OH  activity. 
The  radiometer  looked  in  the  zenith  and  has  a  field  of  view 
of  9*.  The  radiometer  spectral  bandpass  was  centered  at 
1.53  Mm  wavelength. 

Infra.'ed  Isocon  Camera  Results 

The  Isocon  camera  systems'  response  is  in  the  near- 
infrared  from  about  700  to  850  nm.  The  response  curve  and 
the  OH  transitions  which  occur  within  this  bandpass  are 
shown  in  Figure  1-3. 

Figure  5-1  shows  a  video  frame  taken  with  the  large, 
independently-mounted  isocon  camera  at  8: 15  hrs.  UT  on  June 
15,  1983.  The  field  of  view  for  the  photo  is  28*  vertical 
and  37*  horizontal  at  an  azimuth  of  about  320*  and  the 
bottom  of  the  picture  beginning  at  10*  elevation. 


phot 


Examination  o-f  this  -frama  shoMs  7  distinct  bri9ht  OH 


•mission  bands.  Th»  bands  ar*  movin9  northward  in  a 
direction  normal  to  the  bands  <movin9  towards  lower  ri9ht 
hand  corner  of  -frame)  .  These  structural  bands  extended 
across  the  entire  northern  hal-f  o-f  the  sky.  The  structure 
depicted  in  Figure  5-1  was  observed  and  recorded  be9innin9 
at  7:30  hrs.  UT  cantinuin9  until  10! IS  hrs.  UT  on  this  date. 

Usin9  the  video  record,  Taylor  C 1983-841  calculated  the 
apparent  wavelength  and  velocity  o-f  the  band  structure,  with 
techniques  described  by  Hapgood  and  Taylor  119821.  Using 
the  location  and  altitude  o-f  the  observation  site,  radius  o-f 
the  earth,  viewing  angles  of  the  camera,  and  assuming  that 
the  band  structure  resided  on  a  spherical  shell  at  an 
constant  altitude  of  87  km,  the  apparent  temporal  wavelength 
derived  was  2411  km,  and  the  apparent  period  was  1411 
minutes.  Based  on  these  two  values  the  apparent  velocity  is 
2812  meters/second  which  agrees  closely  with  other  similar 
observations  CTaylor  et  al .  19801. 

The  small  isocon  camera  which  was  mounted  on  and 
coaligned  with  the  interferometer,  was  operated  during  the 
same  time  frame  as  the  large  camera.  The  field  of  view  of 
the  small  camera  is  nearly  square  being  13*  vertical  by  IS* 
hor i zontal . 

Figures  5-2  through  5-4  show  a  sequence  of  video  frames 
beginning  at  7:32  hrs.  UT  and  taken  at  8  minute  intervals 
with  the  small  camera.  Each  frame  again  shows,  in  more 
detail,  the  OH  emission  structure.  The  *X"  in  each  frame 
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Fi9ure  5-3.  Isacan  IR  camerai  photo,  day  166,  7:40  hra.  UT, 
field  of  view  »  13*  vertical  IS*  horizontal,  interferometer 
looking  at  ”X"  on  a  dark  band. 


ng  at 


marks  ths  canter  of  the  i  nter 'f  erometer  ’  s  -field  o-f  view 
within  the  video  -frame.  The  interferometer’s  field  of  view 
is  0.8*  full  field  and,  therefore,  is  about  0.4  inches  in 
diameter  on  each  frame.  This  small  field  of  view  permits 
the  interferometer  to  resolve  the  bright  and  dark  portions 
of  the  structure.  The  interferometer  is  viewing  at  an 
elevation  angle  of  17*  and  an  azimuth  angle  of  328*  in  the 
sequence.  The  three  photos  show  one  of  the  periods  where 
the  interferometer  was  able  to  view  a  bright,  a  dark,  and  a 
bright  b  and  in  sequence  as  the  structure  moved  through  the 
field  of  view. 

Figures  S-S  and  5-6  show  another  similar  sequence  of 
video  frames  observed  later  in  the  night  beginning  at  8:31 
hrs.  UT  at  an  interval  of  11  minutes.  This  later  series  of 
frames  again  shows  the  interferometer  viewing  a  dark  then  a 
bright  emission  band  as  the  structure  moves.  The 
interferometer  was  viewing  at  an  elevation  angle  of  15.5* 
and  an  azimuth  angle  of  340*  in  these  last  two  figures. 

Taylor  C 1983-841  compared  the  bright  and  dark  bands 
recorded  by  the  small  isocon  camera  with  the  OH  (3,1)  band 
intensity  plots  computed  using  the  interferometer  data 
(shown  later  in  Chapter  V  and  in  Appendix  C>.  In  the  time 
period  from  7:30  to  10:15  hrs.  UT  the  two  independent 
records  show  that  for  each  bright  or  dark  band  depicted  in 
the  video  data  a  corresponding  increase  or  decrease  in 
intensity  is  also  shown  in  the  interferometer  data.  This 


Figure  5-6.  Isocan  IR  camera  photo,  day  166,  8:42  hra.  UT , 
field  of  view  =  13*  vertical  15*  horizontal,  interferometer 
looking  at  "X"  on  a  bright  band. 
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In-frared  Radiometer  Data 

The  radiometer  was  used  to  show  trend  and  absolute 
radiance  values  for  the  OH  in-frared  activity  as  viewed 
within  a  9*  field  of  view  in  the  zenith.  Radiometer  data 
are  included  to  show  for  the  night  before  (June  14,  1983) 
and  the  night  of  (June  15,  1983)  the  recorded  structure,  OH 
airglow  emission  intensity  exhibited  some  unusual  trends  and 
modulations.  Figure  5-7  shows  the  radiometer  data  for  three 
specific  days  and  is  calibrated  in  k i 1 o-Ray 1 e i ghs  (kR)  of 
intensity  of  OH  (3,1)  band  emission.  The  scale  has  been 
shifted  to  show  the  three  days  all  on  one  chart. 

Day  162  is  a  typical  curve  of  near - i nf rared  OH  activity 
and  was  taken  from  White  Sands  Missile  Range  (4300  feet 
elevation)-  and  located  about  30  miles  from  the  Sacramento 
Peak  site.  As  can  be  seen  from  the  Day  162  curve,  the 
intensity  of  the  1.53  Mm  radiation  steadily  decreases  after 
local  sunset  (sunset  occurred  about  3:30  hrs.  UT )  and 
throughout  the  night. 

The  Day  165  radiometer  curve,  however,  shows  a  markedly 
different  trend  in  the  OH  activity.  The  post-sunset 
decrease  starts,  but  at  about  two  hours  before  local 
midnight  (or  5:00  hrs.  UT)  the  OH  activity  dramatically 
Increased.  After  the  peak,  the  OH  activity  remains  high  but 
also  shows  some  modulation  that  could  be  Interpreted  as 


r" 


-'.V. 

v.v. 


I 


ij 

*  -L-  •. 


f 


{ 


.  53m  CHRNNEL 


90 


distinct  large-scale  wave  structure  passing  overhead.  The 
period  o-f  these  intensity  waves  ranges  -from  30  to  60  minutes 
with  a  modulation  (or  contrast  ratio)  o-f  -from  5  to  10 
percent.  The  waves  were  not  distinguishable  in  the  zenith 

with  the  in-frared  camera  equipment  because  the  zenith 

emission  intensity  is  too  -faint  and  their  large  size 

exceeded  the  camera’s  -field  o-f  view.  A  camera  was  lowered 

in  elevation  in  an  attempt  to  view  the  waves  (lower 
elevations  view  the  structure  obliquely  thus  the  waves 
appear  compressed  together  and  brighter  due  to  the  van  Rhijn 
e-f-fect)  but  no  structures  were  distinguishable  by  the  camera 
throughout  the  entire  night. 

The  radiometer  curve  in  Figure  5-7  for  day  166  also 
shows  enhanced  OH  activity.  Although  not  as  intense  as  the 
previous  night,  a  similar  pre-midnight  maximum  is  observed.. 
After  the  pre-midnight  maximum,  the  decrease  is  slower  than 
normal  and  again  some  distinct  wave  structure  is  observed. 
The  structure  has  periods  of  about  25  minutes  and 

modulations  in  intensity  of  about  5  percent.  Shortly  after 
local  midnight  (7:00  hrs.  UT) ,  the  moon  had  set  far  enough 

to  lower  the  sensitive  camera  to  the  horizon  to  search  for 

observable  structure.  The  measurements  gathered  throughout 
the  remainder  of  the  nignt  with  the  cameras  and 

interferometer  also  exhibit  high-contrast  OH  emission 


structure 


Direct  Comparison  o-f  Interferometer 
and  Radiometer  Intensities 

Prior  to  moonset  on  UT  days  165  and  166  the 
interferometer  Mas  pointed  to  the  zenith  and  was  viewing  the 
sky  coincident  with  the  radiometer.  Therefore  a  direct 
comparison  is  made  between  the  two  records.  Figure  5-8 
shows  the  OH  (3,1)  band  intensity  from  3:30  to  7:15  hrs.  UT 
on  day  165.  The  radiometer  data  were  used  to  provide  the 
absolute  calibration  for  the  interferometer,  but  as  the 
figure  shows  the  trends  in  the  interferometer  intensity 
curve  track  the  radiometer  curve  closely.  The  intensity 
peak  appears  at  about  5:40  hrs.  UT  on  this  day.  The 
standard  deviation  of  the  interferometer  intensity 
calculation  is  about  2%  (the  standard  deviation  curve  is 
included  in  Appendix  C).  The  interferometer  was  realigned 
at  6:30  hrs.  UT.  Just  prior  to  realignment,  the  standard 
deviation  of  the  intensity  increased.  After  alignment  the 
intensity  appeared  to  have  increased.  These  observations 
make  the  rapid  intensity  decline  at  6:00  hrs.  UT  suspect,  as 
an  Instrument  alignment  drift  problem. 

Figure  5-9  presents  the  interferometer  OH  (3,1)  band 
intensity  for  UT  day  166.  Again  it  compares  favorably  with 
the  radiometer  curve  for  the  same  day  in  Figure  5-7.  The 
intensity  peak  occurs  on  this  day  at  5:30  hrs.  UT.  The  same 
reservations  about  alignment  drift  occurred  on  this  day  at 
6: 00  UT  as  happened  on  the  previous  day. 
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Figur*  5-8.  OH  (3,1)  band  radianc*  and  rotational 
tamparature,  viawing  angla  ■>  zanith,  day  165,  3:30-7:30  hrs. 
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Th»  calculated  interferometer  OH  radiance  and 
temperature  curves  for  the  other  OH  Meinel  bands  measured  in 
this  study  are  included  in  Appendix  C.  The  absolute  scales 
of  radiance  for  the  other  OH  bands  Mere  based  on  the  (3,1) 
band  (the  radiometer  was  calibrated  in  terms  of  the  (3,1) 
band)  and  usin^  the  band  intensity  ratios  developed  by 
Turnbull  and  Lowe  119831  ( I (4, 2) (3, 1 ) xl . 03, 
I(8,S)*I (3, l)xO. 14,  I (7,4)«I (3, 1 ) xO. 09) .  Only  when  the 
interferometer  is  vieMing  in  the  zenith  are  the  intensity 
curves  plotted  in  absolute  terms.  The  geometry  of  looking 
at  loM  elevation  angles  for  the  other  data  precluded 
inclusion  of  absolute  scales.  Also  shown  in  the  Appendix 
are  the  standard  deviations  for  all  the  calculations. 

Rotational  Temperatures  for 
Interferometer  Zenith  Data 

Included  in  Figures  S-8  and  5-9  are  the  calculated 
rotational  temperatures  for  days  165  and  166  when  the 
interferometer  is  viewing  in  the  zenith.  The  temperature 
for  day  165  indicates  a  rise  in  temperature  and  intensity, 
beginning  at  about  4:00  hrs.  UT,  with  a  mean  temperature  of 
about  175  *K.  The  total  rise  in  temperature  is  15  *K  with 
the  peak  occurring  at  5:20  hrs.  UT.  The  temperature  curve 
for  day  166  shows  the  same  increasing  trend  as  on  the 
previous  day.  The  mean  is  about  170  *K  with  a  total  rise  in 
temperature  of  15  *K.  The  temperature  peaks  on  day  166  at 
5:10  hrs.  UT.  On  both  days  as  the  interferometer  views  the 
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z»nith,  th«  p^oik  of  th»  temporatur*  curva  pracada*  tha  paak 
of  tha  intanslty  curva  by  about  20  minutas  (discusaad  mora 
i n  Chaptar  VI ) . 


Rotational  Tamparatura 
Smoothing  Algorithm 

Tha  curva  of  calculatad  rotational  tamparaturaa  appaarad 
quita  noisy  at  timas.  Consaquant 1 y ,  smoothing  algorithms 
wars  amployad.  Tha  smoothing  tachniqua  usad  is  a  3  frama 
wida  sliding  window  whara  the  averaging  is  a  weighted  one. 
Only  3  data  points  ware  usad  so  as  to  not  degrade  tha 
temporal  resolution  of  tha  data  (this  is  about  a  2  minute 
time  window)  and  still  provide  sufficient  smoothing.  Tha 
weighting  is  accomplished  using  tha  reciprocal  of  tha 
standard  deviation  squared  (or  reciprocal  of  variance)  for 
each  data  point  as  its  weighting  factor  in  tha  running  sum; 

therefore,  if  a  particular  data  point  has  a  large 

uncertainty  it  is  weighted  lass  in  tha  average.  This 

algorithm  was  chosen  because  in  data  sets  with  large 
differences  from  data  point  to  data  point,  with  significant 
differences  in  tha  variance  of  each  point,  this  weighting 
tachniqua  bast  identifies  tha  mean  curva  through  all  data 
points  CBevington  19691. 

Tha  specific  algorithm  usad  for  tha  rotational 


temperature  smoothing  is 


T. 

1 


(S.  1) 


Z  <T  /<ff  >  ) /  s  (l/{ff  )  ) 
Li*!  *  'i  i«l 


wh*r» 

s  weighted  averag*  temparatura  at  tima  i, 

»  tamparatura  at  tima  i, 

standard  davlation  o-f  tamparatura  . 

Tha  standard  davlation  on  tha  smoothad  curva  Is  also 
changad  dua  to  tha  walghtlng  and  is  racalculated  using 

F-  -  f  E  (l/tff  >^>1  ^  ,  (5.2) 

i  Li.i  i  •• 

whara 

9^^  •  avaragad  standard  davlation  at  tima  i. 

Figure  5-10  shows  tamparatura  and  standard  deviation 
curves  ba-fora  smoothing  and  Figure  5-11  shows  tha  same 
curves  a-ftar  smoothing.  All  tha  tamparatura  curves  in 
Appendix  C  show  a  sat  o-f  tamparatura  curves  -for  both  ba-fore 
and  a-fter  smoothing. 

Intarf aromatar  Recorded  Structure 

This  research  project  was  undertaken  in  an  attempt  to 
quantify  the  intensity  modulations  and  anticipated  OH 
rotational  temperature  modulations  associated  with  airglow 
structure  events  typified  in  Figure  1-4.  The  presentation 
in  this  section  is  tha  result  of  tha  processed 
i ntarf aromatar  OH  (3,1)  band  records,  which  ware  observed 
during  tha  airglow  structure  event  of  June  15,  1983.  Tha 
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OH  (3,1)  band  data  hava  the  hi9hest  signal -ta-noise  ratio 
and,  there-fore,  are  -felt  to  be  the  most  reliable.  The 
presentation  is  given  in  three  segments  corresponding  to 
when  the  instrument  was  moved  in  viewing  elevation  or 
azimuth.  A  complete  catalog  o-f  all  the  observed  OH  band 
intensities  and  rotational  temperatures  is  contained  in 
Appendix  C. 

The  i nter-f erometer- i socon  camera  system  was  lowered  to  a 
viewing  elevation  angle  o-f  17*  and  an  azimuth  angle  o-f  328* 
at  about  7:30  hrs.  UT  (a-fter  moonset)  on  this  date.  Figure 
S-12  show  the  intensity  modulations  (and  standard  deviation) 
as  seen  by  the  i nter-f erometer  until  about  8:30  hrs.  UT.  The 
bright  band,  dark  band,  bright  band  sequence  show  a 
modulation  in  intensity  o-f  about  20X.  This  same  sequence  is 
also  shown  in  the  isocon  video  -frames  in  Figures  5-2  through 
5-4.  Following  the  last  bright  band,  the  intensity  falls 
off  by  40%,  indicating  a  relatively  dark  band.  The 
simultaneous  plot  of  the  OH  rotational  temperature  shown  in 
Figure  5-13  has  a  mean  temperature  of  about  165  *K.  The 
modulations  seen  are  in  phase  and  correspond  to  the 
increases  and  decreases  in  intensity  in  Figure  5-12.  The 
magnitude  of  the  temperature  modulations  are  from  5  to  8  *K. 
The  standard  deviation  on  the  temperature  calculation  is 
about  ±3  *K. 

The  interferometer  was  moved  to  a  different  viewing 
location  at  about  8:30  hrs.  UT  and  remained  there  until 
about  9:15  hrs.  UT.  The  new  viewing  angles  were  15.5* 
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elevation  and  340*  azimuth 


The  intensity  modulations  at 


this  viewing  position  are  shown  in  Figure  S-14.  The  'first 
dark  band,  bright  band,  dark  band  sequence  in  this  figure 
show  a  modulation  of  40X.  The  intensity  modulations  shown 
after  the  first  bright  band  are  much  smaller,  being  on  the 
order  of  lOX.  The  video  frames  in  Figures  5-5  and  5-6 
correspond  to  the  first  dark  band  and  occur  just  after  the 
peak  of  the  very  bright  band  of  Figure  5-14.  The  rotational 
temperature  plot  for  this  same  time  period  is  shown  in 
Figure  5-15.  The  mean  temperature  is  about  165  'K.  with  a 
maximum  modulation  of  lO  •K.  The  uncertainty  on  the 
temperature  calculations  is  about  ±4  *K  during  this  time 
frame.  The  temperature  and  intensity  modulations,  at  the 
beginning  of  this  measurement  period,  are  again  in  phase 
with  each  other.  The  lOX  intensity  modulations  shown  in 
Figure  5-14  occurring  after  8:50  hrs.  UT,  however,  do  not 
show  any  discernible  modulations  in  the  temperature. 

The  last  observations  of  the  night  began  at  about  9:15 
and  lasted  until  lo: 15  hrs.  UT.  The  viewing  position  of  the 
i nterf erometer  was  again  changed  to  15.5*  elevation  and  309* 
azimuth.  Figure  5-16  presents  the  intensity  modulation 
record  for  this  period.  The  curve  shows  the  first  sequence 
of  structure  exhibiting  a  modulat  on  of  40X  and  three  more 
bright  and  associated  dark  bands  with  modulations  of  about 
20X.  The  corresponding  rotational  temperature  plot.  Figure 
5-17,  once  again  has  a  mean  temperature  of  165  *K.  The 


temperature  modulations  are  from  5  to  8  *K  with 


each  rise 
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in  temperature  having  a  corresponding  rise  in  intensity  o-f 
at  least  20%.  The  uncertainty  on  the  temperature 
calculations  is  again  about  14  *K. 

The  presentation  in  this  chapter  o-f  the  interferometer- 
recorded  structure  for  June  15,  1983  is  focused  on  the 
information  extracted  from  the  OH  (3,1)  band.  Examination 
of  the  Appendix  C  records  for  the  other  OH  bands,  presents 
an  additional  observation.  The  two  4v*2  band  temperatures 
track  each  other  within  the  standard  deviation  of  the 
calculations.  The  two  Av«3  band  temperatures  track  each 
other  within  the  standard  deviation  of  the  calculations. 
The  4va3  band  temperatures,  however,  are  consistently  from 
15  to  32  *K  hotter  than  the  Av*2  band  temperatures 
Idiscussed  more  in  Chapter  V/Z).  Figure  5-18  is  an  example 
with  additional  data  available  in  Appendix  C.  Table  5-1 
outlines  a  summary  of  the  results  presented  in  this  chapter. 

TABLE  5-1.  Summary  of  OH  airglow  structure  measurement 
results  for  June  15,  1983. 


1.  Apparent  structure  period  ........  1411  minutes 

2.  Apparent  structure  temporal  wavelength  .  2411  km 

3.  Apparent  structure  phase  velocity  .  .  .  2812  meters/sec 

4.  In-tensity  modulations  measured . lO  to  40  % 

5.  Rotational  temperature 

modulations  measured . 5  to  10  *K 


6.  Phase  relationship  of  recorded  rotational 

temperature  and  intensity  measurements  ....  In  phase 


7 


Mean  zenith  OH  (3,1)  band  intensity 
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rotational  papulations  ara  in  trua  tharmal  aqui librium  and 
are  thus  strictly  Boltzmann  distributad.  A  slight  daviation 
from  this  assumption  would  causa  a  systamatic  arror. 
Anothar  passible  source  of  error  results  from  the  assumption 
that  all  OH  airglow  radiation  is  emitted  from  a  thin  uniform 
layer)  whereas  in  reality,  the  layer  is  about  7  km  in 
thickness.  In  addition,  at  low  viewing  elevation  angles 
%10*  the  layer  geometry  is  much  more  complicated.  At  these 
low  elevation  angles  atmospheric  extinction,  van  Rhijn 
effect,  and  curved  spherical  geometry  potentially  have  a 
significant  impact  on  the  interpretation  of  the 
measurements . 

The  interferometer  spectral  response  calibration  is  very 
sensitive  to  alignment.  The  instrument  typically  remained 
in  alignment  for  about  2  hours.  As  can  be  seen  from  the 
increase  in  the  standard  deviation  as  a  function  of  time 
(7; 30  to  lO: IS  hrs.  UT,  day  166),  the  alignment  changed 
significantly  and  this  change  could  account  for  a  portion  of 
the  uncertainty.  Therefore,  the  temperature  modulations 
obtained  from  the  spectral  data  are  felt  to  be  more 
accurately  defined  than  is  suggested  by  the  standard 
dev i at i on . 

The  assumption  that  the  low  observation  elevation  angles 
associated  with  the  structure  measurements  did  not  unduly 
impact  the  computed  temperature  is  supported  by  the  fact 
that  the  temperature  at  6:4S  hrs.  UT  on  day  166  was  163  *K 
(see  Figure  5-9)  and  the  temperature  at  7:30  hrs.  on  the 


same  day  was  166  *K  <see  Figure  5-13).  These  temperatures 
represent  the  values  computed  just  be-fore  and  just  after  the 
interferometer  "look  direction"  was  changed  from  the  zenith 
to  near  the  horizon. 

Rotational  Temperature  and 
Intensity  Modulations 

The  ranges  of  temperature  and  intensity  modulations 
observed  in  the  OH  Melnel  airglow  structures  are  given  in 
Chapter  V).  The  "adiabatic  oscillation"  and  "IGW"  modeling 
of  the  OH  Meinel  airglow  variations  mentioned  in  Chapter  I 
utilize  a  parameter  which  is  readily  calculated  from  the 
intensity  and  temperature  modulations.  This  parameter  is 
the  ratio  of  the  change  in  emission  intensity  normalized  by 
the  mean  emission  intensity,  divided  by  the  change  in 
temperature  normalized  by  the  mean  temperature  and  is  deemed 
useful  in  studies  of  the  OH  airglow  structure  phenomena. 
The  parameter  is  usually  represented  by  the  Greek  letter  eta 
(t))  and  is  defined  as  follows: 

»  C  (6l/i)  /  (aT/T )  3  ,  (6.1) 

where 

Al  a  change  or  modulation  in  emission  intensity, 

1  *  mean  value  of  the  emission  intensity, 

AT  >  change  or  modulation  in  rotational  temperature, 

T  •  mean  value  of  the  rotational  temperature. 


The  value  o-f  is  potentially  use-ful  in  d  ist  i  nguish  i  n9 
between  chemical  processes  which  give  rise  to  the  OH  airflow 
emission  and  temperature  structure.  The  physics  o-f  this 
parameter  is  discussed  by  Krassovsky  C 19723  and  Weinstock 
C197S1.  Pendleton  C 198S3  has  summarized  the  essential 
-features  of  this  parameter  in  Figure  6-1.  In  this  figure, 
the  11  value  is  plotted  versus  the  ratio  (H/H  ),  where  H  is 

X 

the  appropriate  atmospheric  scale  height,  and  is  the 
scale  height  (near  8S  km)  of  minor  species  "x".  Here  the 
letter  "x"  represents  either  oxygen  (O)  or  hydrogen  (H) . 
The  simple  adiabatic-oscillation  model  of  Krassovsky  C1972I 
yields  i)  values  which  are  independent  of  <H/H^),  whereas  the 
gravity-wave  model  of  Weinstock  C19783  yields  (H/H^) 
dependent  values.  The  range  of  i)  values  expected  on  the 
basis  of  values  of  (H/H  )  inferred  from  several  measured 

X 

atomic  oxygen  profiles  is  also  shown  in  the  figure.  The 
information  in  Figure  6-1  indicates  that  values  of  i)  in  the 
range  from  3  to  6  might  be  expected  based  on  current 
gravity-wave  modeling  and  the  ozone  hydration  process. 

Using  the  numbers  for  the  intensity  and  temperature 
modulations  and  means,  presented  in  Chapter  V,  the  range  of 
calculated  values  for  i)  are  f''om  8  to  12.  These  values  are 
about  a  factor  of  2  greater  than  those  shown  in  Figure  6-1. 
In  view  of  the  relatively  large  standard  deviations  on  the 
temperature  determinations,  the  nominal  factor  of  two 
disparity  between  predicted  and  measured  i)  values  is  not 
deemed  significant. 
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Figur*  6-1.  Expactvd  valuvs  for  i)  (vortical  axis)  basad  on 
gravity-wava  modals  and  oxygan  maasuramants  CPandlaton 
198S1 . 


The  1)  values  calculated  for  the  time  periods  when  the 
interferometer  was  viewins  in  the  zenith  were  obtained  by 
assumins  the  pre-midnight  increase  in  both  temperature  and 
intensity  reflected  a  wave-like  disturbance.  The  zenith¬ 
viewing  T)  values  associated  with  the  major  pre-midnight 
(I,T)  fluctuations  on  UT  days  165  and  166  were  found  to  be 
consistent  with  values  deduced  from  the  low-elevation- 
angle  data.  The  consistency  of  these  two  sets  of 
calculations  lends  credence  to  the  idea,  once  again,  that 
viewing  near  the  horizon  had  little  impact  on  the  i) 
determinations  although  the  modulation  in  both  intensity  and 
temperature  may  have  been  impacted  by  the  geometry  of  the 
measurements . 

The  field  of  view  <FOV)  of  the  interferometer  is  0.8* 
full  field.  Consideration  of  the  OH  structures  in  the 
nominally  13*  by  IS*  video  frames  suggests  that  the  0.8*  FOV 
of  the  i nterf erometer  results  in  a  horizontal  spatial 
integration  over  about  one  half  cycle  in  the  quasi-period 
structures.  This  integration  will  degrade  the  horizontal 
spatial  resolution  of  the  interferometer  measurements.  If 
the  wave  is  assumed  to  be  sinusoidal  in  nature  and 
restricted  to  a  very  thin  spherical  shell,  a  simple 
integration  over  1/2  cycle  indicates  that  the  rotational 
temperature  modulations  could  be  degraded  by  a  factor  of 
about  1.5.  The  horizontal  intensity  structure  is  expected 
to  be  more  complex  than  the  temperature  structure  CWeinstock 
19781.  However,  if  a  similar  "degradation  factor"  were 


applicable  to  the  intensity  measurements,  then  the 
sat  is-f  actory  agreement  o-f  the  values  deduced  from  the 
zen i t h - V i ew i ng  and  near -hor i zon-v iew i ng  measurements  could 
be  explained. 


Temperature  Differences  Observed 
Between  Bands 

Rotational  temperature  differences  between  the  high-v' 
and  low-v'  levels  characterized  both  the  zenith-viewing  and 
near-hor i zon-v iewi ng  measurements.  The  high-v'  rotational 
temperatures  were  consistently  higher  than  those  obtained 
from  low-v'  bands.  The  bands  used  for  detailed  comparison 
are  the  high-v'  <7, 4)  band  and  the  low-v'  (3,1)  band.  These 
bands  were  chosen  because  of  favorable  instrument  response 
and  alignment  stability  in  the  spectral  regions  of 
occurrence.  The  temperature  differences  observed  ranged 
between  IS  and  32  'K.  The  smaller  difference  applied  when 
looliing  in  the  zenith,  and  the  difference  gradually 
increased  throughout  the  night  as  the  telescope  was 
physically  moved,  affecting  the  alignment.  The  increase  in 
the  temperature  difference  is  largely  attributed  to 
instrument  alignment  drift  because  the  standard  deviations 
on  the  temperature  calculations  (whici  also  reflect 
misalignment)  increase  simultaneously  with  (and  at  about  the 
same  rate)  the  increase  in  temperature  difference. 

Explanations  were  sought  for  this  difference.  It  was 


discovered  that  an  error  in  the  line  strength  constant  for 


the  OH  (7,4),  Pt  (4)  line  had  been  entered  into  the 
processing  system.  The  constant  was  about  20%  low  in  value. 
In  order  to  assess  the  impact  o-f  this  error,  a  synthetic  OH 
spectrum  was  generated  assuming  a  typical  temperature  with 
the  associated  error  entered,  and  a  Boltzmann  plot  was  made. 
A  line  was  -fit  to  the  points  in  a  least-squares  sense  and  an 
associated  temperature  extracted.  The  OH  (7,4)  band  model 
used  12  lines  in  the  -fit:  there-fore,  because  the  model  -fit 
is  least-squares  in  nature  the  error  from  this  incorrect 
constant  was  found  to  be  less  than  1%. 

Another  possible  explanation  of  the  high-v',  low-v' 
temperature  difference  is  slightly  different  (%  1-3  km) 
emitting  altitudes.  The  ratio  of  aT  to  the  mean  temperature 
f,  if  different  for  the  high-v’  low-v*  measurements,  would 
support  such  a  difference  in  mean  emitting  heights 
CPendleton  19853.  The  AT/T  ratio  within  each  observational 
time  frame  was  calculated  and  the  difference  in  the  ratio, 
between  the  high  and  low  rotational  levels,  was  found  to  be 
less  than  10%,  with  less  than  30%  difference  among  all 
frames.  These  differences  in  the  AT/T  ratio  between  high-v' 
and  low-v*  bands  when  viewed  in  terms  of  the  calculated 
standard  deviations  does  not  provide  evidence  for 
differences  in  mean  emitting  altitud.‘. 

Within  the  standard  deviation  of  the  calculations,  the 
rotational  temperature  differences  between  the  high-v*  and 
low-v*  levels  appears  to  be  real.  The  magnitude  for  this 


d  i-F-f  erence,  as  was  mentioned  in  Chapter  I 


is  within  the 


range  reported  by  Krassovsky  and  Shagaev  C 19773. 

Temperature  and  Intensity 
Phase  Relationship 

The  modeling  o-f  IGU’s  mentioned  earlier  predicts  that 
the  changes  in  the  OH  Meinel  rotational  temperature  should 
be  in  phase  with  the  IGW  CHines  19601.  The  change  in 
intensity,  however,  should  be  related  to  the  IGW  temporal 
structure  in  a  potentially  more  complicated  manner.  This 
situation  arises  partly  because  o-f  the  -finite  chemical  time 
constant  associated  with  mesospheric  ozone.  The  appropriate 
time  constant  for  the  cool  <T  *  160  *K)  summer -mesopause 
conditions  is  about  25  minutes  CPendleton  19853.  For  IGW 
periods  somewhat  in  excess  of  this  value,  chemical 
conversion  of  is  expected  to  be  a  significant  factor  in 
the  phenomenology,  whereas  for  much  smaller  IGW  periods 
chemical  conversion  should  be  unimportant.  Examination  of 
Figures  5-8  and  5-9  shows  that  the  temperature  maxima  lead 
the  intensity  maxima  by  about  20  minutes.  This  apparent 
phase  difference  may  relate  to  the  aforementioned  O  time 
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constant,  but  it  would  be  premature  to  draw  this  conclusion. 
It  is  suggested  that  additi.jnal  attention  be  given  to  this 
interesting  possibility. 

The  small-scale  structures  observed  near  the  horizon  on 
UT  day  166  exhibit  an  in-phase  relationship  between 
Intensity  and  temperature.  The  only  exception  to  this  is 
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when  the  measured  intensity  modulation  was  less  than  about 


20%.  Under  these  conditions,  no  direct  correlation  between 
intensity  and  temperature  could  be  drawn.  Since  the 
measured  period  ot  the  small-scale  OH  Meinel  structures  was 


CHAPTER  VI 1 


CONCLUSIONS  AND  RECOMMENDATIONS 


Overview 


The  goal  o-F  this  study  was  to  design,  develop,  and 
operate  an  instrument  system  capable  o-f  performing 
simultaneous  spatial,  spectral,  and  temporal  high-resolution 
OH  airglow  measurements.  The  design  herein,  and  the 
resulting  data  demonstrate  the  effectiveness  of  the 
technique.  An  airglow  structure  event  which  occurred  on 
June  IS,  1983  was  measured  with  the  interferometer- 
spectrometer  system.  OH  Meinel  intensities  and  rotational 
temperatures  were  obtained  for  the  peaks  and  troughs  of  this 
wavelike  structure. 


Conclusions 


The  following  are  the  specific  accomplishments  of  this 
study.  The  areas  addressed  pertain  both  to  the  instrument 
designed  for  airglow  structure  measurements  and  to  the  data 
processing  techniques  used. 

1.  An  optically-compensated  interferometer  for  high 

2 

throughput  (Aa>0.28S  cm  sr ) ,  was  matched  to  a  large 
area  collector  (SO-cm  diameter)  to  narrow  th'e  field  of 
view  (0.8* >.  A  noise  equivalent  spectral  radiance 
(NESR)  (sensitivity)  of  16  R/cm  ^  at  l.S  Mm  was 
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achieved.  In  comparison,  a  conventional  Michelson 
i  nter-f  erometer-spectrometer  with  the  same  detector, 
collector  area,  resolution,  and  scan  time  Mould  have  a 
NESR  o-f  128  R/cm  a  factor  of  8  less  sensitive  than 
the  one  developed  for  this  study.  When  compared  with  a 
conventional  Ebert  spectrometer,  using  the  same  detector 
and  operated  at  the  same  resolution,  a  sensitivity  of 
208  R/cm  ^  could  be  achieved,  a  factor  of  13  less 
sens! t i ve. 

2.  A  spectral  resolution  of  2  cm  ^  was  sufficient  to 
resolve  the  OH  emission  line  structure  for  the 
extraction  of  OH  rotational  temperatures.  Based  upon  a 
rotational  line  separation  of  10  cm  ^  and  the  Hamming 
apodization  function  used,  the  spectral  resolution  of 
the  instrument  could  be  lowered  to  no  more  than  4  cm  ^ . 

3.  Based  upon  the  video  records,  the  bright  or  dark  bands 

of  OH  structure  subtend  about  1*  of  arc  at  these  low 
elevation  angles  IS* > .  The  apparent  temporal 

wavelength  was  24i 1  km,  with  a  period  of  14±1  minutes, 
and  an  apparent  phase  velocity  of  28±2  meters/second. 

4.  The  interferometer  system  field  of  view  was  measured  at 
0.8*.  The  Interferometer  FOV  is  sufficiently  narrow  to 
independently  view  a  "bright"  or  a  "dark"  structure 
band.  Based  on  this  limited  data  set  and  simple  wave 
geometry,  the  field  of  view  could  be  as  large  as  7*  if 
these  structures  were  viewed  in  the  zenith. 


Th»  measur*d  intansity  modulations  (contrast  ratios)  for 
the  OH  airglovi  structures  ranged  between  20  and  40%  with 
interferometer  recorded  periods  of  14  minutes  The 
calculated  standard  deviation  was  typically  3%. 

The  mean  calculated  OH  Meinel  rotational  temperature  for 
the  aforementioned  event  was  165  *K.  The  measured 
modulation's  in  rotational  temperatures  associated  with 
the  changes  in  intensity  ranged  from  5-10  *K  and  are  in 
phase  with  the  intensity  modulations.  Typical  standard 
deviations  on  the  rotational  temperature  calculations 
ranged  from  2-7  *K.  The  mean  temperature  and  magnitude 
of  the  temperature  fluctuations  are  consistent  with  both 
IGU  theory  and  previous  mesospheric  temperature 
measurements. 

The  least-squares  model  used  to  extract  band  intensity 
and  rotational  temperature  provides  a  computational 
efficient  way  (convergence  to  final  values  occurred 
within  4  iterations  of  the  fitting  routine)  to 
simultaneously  derive  these  values.  The  model  also 
provides  a  measure  (standard  deviation)  of  how  well  the 
data  fit  a  Boltzmann  distribution. 

The  insight  needed  to  identify  what  is  being  observed 
from  the  air  glow  layer  can  only  be  provided  by  the 
camera  (or  similar)  video  system.  The  measurement  of  OH 
airglow  structure  events  with  the  interferometer  system 
would  not  be  possible  without  the  simultaneous  use  of 
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the  isocon  camera,  bacausa  Mith  tha  intarf aromatar 
alona,  exactly  what  was  bain9  viewed  would  be  unknown. 

Recommendations  ‘for  Future  Research 

The  interferometer  worked  as  designed  and  proved  to  be 
an  excellent  tool  for  this  type  of  study.  The  model 
developed  for  the  extraction  of  rotational  temperatures  and 
intensities  is  an  accurate  technique  and  provides 
computational  flexibility.  However,  several  suggestions  are 
made  for  consideration  for  future  work. 

1.  The  optical  path  within  the  interferometer  is  very 

complex.  There  are  20  optical  surfaces  through  which 
the  incoming  energy  must  pass  before  reaching  the 
detector.  Assuming  a  typical  loss  of  4%  per  surface, 
80%  of  the  incoming  signal  is  lost  before  reaching  the 
detector.  This  complex  optical  path  should  be 

redesigned  to  minimize  the  number  of  optical  elements 
and  optically  coat  the  remaining  elements  to  minimize 
reflection  loss. 

2.  The  physical  size  of  the  instrument  should  be  reduced  to 
facilitate  portability  to  remote  sites. 

3.  The  instrument  is  very  sensitive  to  optical  alignment. 
The  optical  components  need  to  be  mounted  in  a  more 
stable  manner  to  hold  their  location  better. 

4.  The  alignment  is  very  sensitive  to  temperature.  It  is 
necessary  therefore  to  temperature  control  the 
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changes  as  the  outside  temperature  -fluctuates. 

The  telescope  pointing  system  needs  to  be  automated. 
The  sensitivity  of  alignment  to  physical  motion  as  well 
as  the  need  to  point  the  instrument  to  an  area  of 
interest  demand  that  the  positioning  of  the  system  be 
automated . 

An  infrared  camera  system,  like  the  one  provided  by  the 
University  of  Southampton,  needs  to  be  permanently 
incorporated  as  part  of  the  interferometer  system.  The 
interferometer  cannot  effectively  gather  data  on  airglow 
structure  if  the  location  of  the  structure  is  not  known 
and  the  camera  provides  this  input. 

The  instrument  must  be  provided  a  better  means  of 
calibration.  The  blackbody  sources  used  in  this  study 
give  a  reasonable  indication  of  alignment  and  instrument 
response  but  as  the  system  alignment  drifts  the 
calibration  is  less  meaningful.  Perhaps  a  technique 
utilizing  OH  spectral  line  pair  ratios  which  are 
independent  of  rotational  temperature  but  sensitive  to 
alignment  could  be  used  as  a  dynamic  measure  of 
instrument  alignment. 

A  recommendation  is  made  to  investigate  other  detectors, 
in  order  to  extend  the  ability  of  the  interferometer, 
with  a  wider  spectral  bandwidth  and  hi-aher  sensitivity. 
The  RCA  detector  used  in  this  study  is  an  excellent 


detector  where  it  is  sensitive  but  is  somewhat  limited 
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L2.  The  extraction  o-f  line  amplitudes  'from  the 
i ntcr-f erometer  could  be  improved  in  two  areas.  First, 
the  approximate  locations  o-f  the  spectral  lines  of 
interest  are  found  in  a  manual  manner.  A  template  usin9 
a  synthetic  spectrum  could  be  designed  for  each  OH  band 
and  a  correlation  routine  could  be  used  to  automatically 
search  the  raM  data  for  the  location  of  the  lines. 
Secondly,  the  apodization  routine  used  to  extract  the 
actual  line  amplitude  from  the  data  should  be  modified 
to  calculate  a  line  area  rather  than  amplitude.  The 
area  routine  would  provide  for  the  averaging  out  of 
noise  whereas  the  amplitude  routine  always  searches  for 
the  most  positive  peak. 

13.  The  model  should  be  modified  to  to  use  the  *Q*  branches 

of  the  OH  bands.  To  do  this  the  molecular  constants  for 

both  Q.  and  would  need  to  be  averaged  as  one  because 
12 

the  interferometer  does  not  resolve  the  two  groups. 
This  addition  should  add  more  accuracy  to  the  model 
because  the  Q  branches  are  the  largest  lines  within  each 
band . 

14.  The  model  could  be  modified  to  include  a  third  variable, 
water  column  content,  based  on  the  several  lines  within 
the  OH  bands  which  are  severely  affected  by  water 
absorption.  The  model  then  could  provide  additional 


information  about  the  atmosphere 
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Appendix  A 
OH  Transitions 

This  appendix  brie-fly  describes  why  the  radiation  from 
the  OH  molecule  is  so  complex.  The  many  spectral  lines 
generated  by  the  excited  radical  are  depicted  in  Figure  4-1. 
The  molecule  rotates  and  vibrates  s Imu 1 taneous 1 y «  and  each 
of  the  motions  are  quantified  and  interactive. 

The  total  molecular  angular  momentum  without  electronic 
spin  i<  is  also  quantized  and  is  identified  by  the  quantum 
number  K.  The  quantity  K  is  comprised  of  two  parts 

^ 

K  *  N  +  A  .  (A. 1 > 

The  vector  H  is  the  nuclear  angular  momentum  and  the  vector 
A  is  the  angular  momentum  of  the  orbiting  electron  cloud 
projected  onto  the  internuclear  axis.  The  quantum  number  A 
associated  with  the  electronic  orbital  momentum  can  take  on 
a  value  of  -t-l  or  -1  depending  upon  which  way  the  electron 
cloud  is  orbiting  with  respect  to  the  nuclear  rotation.  The 
double  degeneracy  of  A  leads  to  the  so-called  A  splitting  of 
each  state;  however,  the  split  is  less  than  1  cm  ^  at  low 
rotational  speeds  CBaker  19783,  which  is  less  than  the 
instrument  resolution  used  for  this  study,  therefore  the  A- 
sp 1 i t  lines  will  be  considered  as  one. 

The  quantum  number  K  can  take  on  values  K  =  1,2,3,... 

The  selection  rule;  however,  is 
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6K  =  O,  ±1 


(A.  2) 


The  collection  o-f  lines  within  each  band,  grouped  according 
to  their  respective  AK,  are  called  branches.  The  branch 
with  AK  =  0  is  called  the  Q  branch,  that  -for  AK  =  +1  is  the 
P  branch,  and  that  -for  AK  =  -1  is  called  the  P  branch. 

The  OH  molecule  has  an  odd  number  o-f  electrons.  This 
imbalance  results  in  a  net  electronic  spin  angular  momentum 
S  and  is  represented  by  quantum  number  S.  The  odd  electron 
gives  rise  to  an  even  multiplicity  2S^1.  Since  the  total 
number  o-f  electrons  is  odd,  S  is  hal-f  integral  (S  =±1/2), 
each  transition  state  is  a  doublet.  It  is  sometimes 
convenient  to  consider  the  total  electron  angular  momentum  A 
as  a  separate  entity.  The  total  electron  angular  momentum 


n  =  A  +  S 


(A.  3) 


There-fore,  each  vibration-rotation  transition  will  split 
into  two  separate  spectral  lines  according  to  whether  A=3/2 
or  A=l/2. 

The  OH  molecule  is  very  light  and  as  a  consequence  the 
odd  electron  spin  3  is  only  weakly  coupled  to  the 
i  nternuc  lean  axis.  Thj>  molecule  is  there-fore,  modeled  as 
Hand’s  case  (b)  CHertzberg  19713.  The  total  molecular 
angular  momentum  5  can  now  be  formed 


•t 

5  =  K  -►  S 


(A. 4) 


I 
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As  can  be  seen  -from  Equation  A. 4,  -for  each  value  o-f  K  there 
S’^e  two  values  tor  J.  Each  branch  o-f  the  OH  spectra  must 
also  tat<e  on  two  values.  It  A=3/2,  then 

J  =  K  +  1/2  =  1.3,  2.5,  3.3,  ...  .  <A.3) 

These  values  -for  J  lead  to  a  set  o-f  spectral  lines  known  as 

P^,  Qj,  and  branches.  I-f  11=1/2,  then 

J  =  K  -  1/2  *  0.5,  1.3,  2.5,  ...  .  (A. 6) 

These  values  for  J  give  rise  to  a  set  of  spectral  lines 

known  as  *3^,  and  R^  branches.  Additional  information  on 

the  physics  of  OH  molecule  is  readily  available  in  the 
literature,  among  them  are  Baker  119783,  Hertzberg  119713, 


and  Mies  119743 
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Appendix  B 

RCA  Limited  Germanium  Detector  Spec  i -f  i  cat  ions 


Senes  TS-71 
TM  71 

GERMANIUM  PHOTODIODE  PREAMPLIFIER  SYSTEMS 
For  Detection  of  Radiation  at  Wavelengths  from  0.8  —  1.7  Micrometers 
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-  NEP  (1.42.  SOO.  1)  =  10  -  •  WHl  '  ' 
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-  Cooitd  fo  77  K 

-  Raspoftsivity  10  •lO'VW  *' 


B 

In  order  to  achieve  the  good  noise  performance,  it  is  necessary  ro  ccoi  oolh  ine  oetector  anc  ;re  pre¬ 
amplifier.  Dewar  A  is  for  liquid  nitrogen  coolant;  the  r»oia  time  -s  acoui  i0  nours.  tne  we  cni  is  wCprox- 
imateiy  3  pounds  and  the  overall  size  is  rougnly  9  cm  diameter  ov  16  cm  lonq.  Dewar  3  is  fc'  i  cuid 
nitrogen,  freon.  COj.  or  any  other  convenient  coolant:  tne  ncid  t-me  for  Mcuid  nitrccen  is  aocut  T 
hours,  weight  is  about  7  pounds,  and  the  overall  sizeis  rougnly  1 1  cm  diameter  by  Cl  cm  long  C  tner 
Dewar  allows  thruput  from  tne  siOe  or  bottom  «f  specified. 


Sceoai  optics  Th«  stanoarp  window  <s  ouaru.  Oifferem  wmoow  material,  scec  a  •  or  co'^ce'^s  '--  jdt  cs  cj''  -e  I'tied 
provided  no  suostantiai  mechanical  redesign  is  necessary.  Tne  norrr^ai  fteid  o»  v.e«v  icr  tr*  oetector  is  c  c  .p  to  90- 


If  the  customer  desires  fo  operate  the  rw-Zi  system  at  very  hign  frpcuenc  es  t  s  coss  o'e  fo  I'.ioe  for  "'eduoncv  res¬ 
ponse  Other  special  features  may  oe  aooed  to  these  oeveiocmentai  svsiems  T-e  r^ecessarv  cr'.mqes  can  se  m.i.ie  for  a 
s.maii  additional  cost. 
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PCA  Lim-ted  —  Sesparcn  Ljboroiories 
S:e-Annp.de-Be  evue.  Due. 
reiepnone  tSt^i  463*3f000 


The  resDonsitiv>lv  vanes  accoffli'^g  to  the  scectrai  curve  Sorre 
variation  n  ins  curve  'S  aossice  accorc-ng  to  me  customer 
reauirements. 
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Figure  B-1.  Typical  detector  spectral  response 


TABLE  B-1.  Germanium  detector  technical  data. 


Operating  Temperature 

Noise  equivalent  power 
(guaranieeo)  WHa  -  * 

Best  noise  equivalent  power 
(achieved  m  the  past)  WHa  - 

Noise  level  (at  output) 

Impedance  level  (at  output) 

Responsivity  v  w 

Linear  range  tor  power  w 

Useable  limit  (power)  w 

Frequency  characteristic 

Detector  area  (circular)  cm- 


TESena^  71 
T7  K 

NEP  (1.42.10.1)^1110 

NEP(1.42.23.1)s1h10 
'1-V 
-  5000 

-  5  I  10 
10  *  -  10 

-  6i  10  • 

I 

0.2 


TM  Series  71 
77  K 

NEP(V42.S00.1)*  IslO* 

NEP(1. 42. 4000.1)  s  3al0 

^  to  ,V 
-  SOC 
-  5  ■  10 

10  ‘  -  to 
^  2  1  10  • 
flat.  3CB  at  SOO  Hz 
0.2 


Operating  volts  (typical)  v 
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Appendix  C 


I nter -f er ometer  Data  Catalog 


The  i  nter-f  erometer  data  -for  observations  taken  on  June 
13-15,  1983  are  contained  within  this  appendix  in  their 
entirety.  All  data  presented  were  recorded  at  Sacramento 
Peak,  New  Mexico.  The  figures  are  all  organized 
chronologically  beginning  on  day  165  at  3:30  hrs.  UT  and 
proceeding  through  day  166,  10:15  hrs.  UT . 

Both  of  the  observation  days’  records  begin  with  the 
interferometer  viewing  in  the  zenith.  At  moonset  on  day  166 
the  interferometer  and  camera  systems  were  lowered  to  view 
near  the  horizon.  The  low  elevation  viewing  period  during 
day  166  is  divided  into  three  time  frames,  corresponding  to 
when  the  interferometer  was  adjusted  in  viewing  location. 
The  time  period  and  viewing  position  are  identified  in  each 
figure  caption. 

Within  each  time  segment  of  the  data  presentation,  the 
figures  are  organized  according  to  the  OH  Meinel  band  from 
which  the  data  were  calculated.  First  is  the  OH  (4,2)  band, 
second  the  OH  (3,1),  third  the  OH  (8,5)  band,  and  fourth  the 
OH  (7,4)  band.  Each  of  the  band  groups  shows  curves  for 
first  the  intensity,  second  the  rotational  temperature,  and 
third  the  smoothed  rotational  temperature  calculation  with 
the  associated  standard  deviation  for  each. 
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Figure  C-3.  OH  (4,2)  band  smoothed 
and  standard  deviation,  viewing 
3:30-7:30  hrs.  UT. 
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Figure  C-6.  OH  (3,1)  band  smoothed  rotational  temperature 
and  standard  deviation,  viewing  angle  =  zenith,  day  165, 
3:30-7:30  hrs.  UT. 
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Figure  C-7.  OH  (8,5)  band  radiance  and  standard  deviation 
viewing  angle  =  zenith,  day  165,  3:30-7:30  hrs.  UT. 
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Figure  C-12.  OH  (7,4)  band  smoothed  rotational  temperature 
and  standard  deviation,  viewing  angle  ■  zenith,  day  165, 
3:30-7:30  hrs.  UT. 
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Figure  C-29.  OH  (3,1)  band  rotational  temperature 
standard  deviation,  viewing  angle  =  17*  El.  328*  Az 
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Figur*  C-30.  OH  (3,1)  band  smoothed  rotational  temperature 
and  standard  deviation,  viewing  angle  s  17#  328*  Az., 
day  166,  7:30-8:30  hrs.  UT. 
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Figure  C-32.  OH  (8,5)  band  rotational  temperature 
standard  deviation,  viewing  angle  =  17"  El.  328*  Az . , 
166,  7:30-8:30  hrs.  UT. 
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Figure  C-33.  OH  '3,5)  band  smaothed  rotational  temperature 
and  standard  deviation,  viewing  angle  =  17*  El.  328*  Az., 
day  166,  7:30-8:30  hrs.  UT. 
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Flgur*  C-38.  OH  (4,2)  b«n4  rotational  tamperatur*  a 
standard  daviation,  viaMing  angla  «  15.5*  El.  340*  Az.,  d 
166,  8:30-9:15  hrs.  UT. 
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Figur*  C-39.  OH  (4,2)  band  smoothad  rotational  tamparatura 
and  atandard  daviation^  viawing  angla  *  15.S*  El.  340*  Az., 
day  166,  8:30-9:13  hra.  UT. 
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Fi9ur*  C-46.  OH  (7,4)  band  ralativ*  intansity  and  standard 
daviation,  viawing  angla  «  15.5*  El.  340*  Az. ,  day  166, 
8:30-9:13  hrs.  UT. 
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Figure  C-49.  OH  (4,2)  band  relative  intensity 
deviation,  viewing  angle  »  15.5*  El.  309* 
9:15-10:15  hrs.  UT. 
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Flgur*  C-33.  OH  (3,1)  b»n<l  rotational  tamparatur* 
atandard  daviation,  viawing  angla  «  IS. 5*  El.  309*  Az., 
166,  9:l5-10:i5  hrs.  UT. 


a  a 


Figure  C-56.  OH  (8,S)  band  rotational  temperature  and 
standard  deviation,  viewing  angle  =  IS. 5*  El.  309*  Az.,  day 
166,  9:iS-10:i5  hrs.  UT. 
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Figure  C-S8.  OH  (7, <47  band  relative  intensity  and  standard 
deviation,  viewing  angle  «  15.5*  El.  309*  Az . ,  day  166, 
9: 15-10; 15  hra.  UT. 
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Figure  C-60.  OH  (7,4)  bend  smoothed  rotational  temperature 
and  standard  deviation,  viewing  angle  «  15. S*  El.  309*  Az., 
day  166,  9:15-10:15  hrs.  UT. 
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F0PMAT< 16) 


P£AD  IN  nOL£CULPP  CONSTANTS  AND  R£LATIV£  RESP0S£ 


DO  110  I‘1,ITL 

R£AD( 5. J20)NAVN( I > ,RLRS( I ) ,F(I> .C(I) 
FORMAT <B2 5.6. £1 5.6, BlS.o.DJS.S) 
CONTINUE 


^EAD  IN  TOTAL  NUMBER  OF  FRAMES  FROM  LINE  AMP  FILE 


READ(6.  IOOUFMS 
URITEi? ,Z05)1FMS 
FORMAT (14) 


BE61H  TO  PROCESS  EACH  FRAME  OF  DATA  IN  TURN 
DO  9000  IF *2. IF  MS 
READ  IN  HEADER  INFO 


DO  200  JF*2.20 
READ(6,200)IHDTA(JF ) 
CONTINUE 


READ  IN  LINE  AMPLITUDES 


BO  220  JF-2,ITL 
F.EAD<a.220)SN.XLSTR<JF> 
F.)RM0T<F20.2  .£25.1) 
CONTINUE 


HR1TE<4.60'>IHDTA(1)  .IHDTACZ) 

FORMAT ('  FRAME  ’ . I S , ’ - ’ . 13 . '  HAS  BEEN  READ') 


SCALE  LINE  AMPLITUDE  DATA  FOR  FFT  «  POST  AMP  6AINS 


IFFT‘IHDTA(S) 

IGM*IHDTA( 20) 

CALL  GNSCLilFFT.IGN.ITL.XLSTR) 


SET  UP  INITIAL  VALUES  FOR  PROCESS 


MDA'24 

MDB^24 


c 


FOPTSiiH  yj  pp.,j6P.f)» 


I 


PIT,  4  >-0.0 

1  '1' 

flT(7)=-1.0 

liS 

f:T<.s>^i.o 

2S4 

SOTO  5000 

:ss 

C 

/56 

C 

7-4  SAND 

2S7 

C 

JS8 

4000 

:SND*5 

1S9 

Bf^<45) 

160 

1S6N‘39 

161 

I£ND^S2 

162 

»0-12 

1 6Z 

fJ7<l'>‘1.0 

264 

fIT<2>^1.0 

i  u  5 

TIT<.3>  =  1,0 

166 

PIT^  4 )=1.0 

167 

PIT< 10)^1 .0 

16S 

F17( 11>‘1.0 

169 

tlT( !3>‘1.0 

170 

FI7(14}=1.0 

171 

6070  5000 

172 

c 

173 

9000 

CONTINUE 

174 

60  70  10000 

175 

c**** 

176 

c 

TENAENATUAE  CALCULATION 

177 

c 

17S 

5000 

70*190.0 

17  9 

B0*NCI</T0 

ISO 

N*I£N0-IBSN*2 

ISl 

N*2 

1S2 

HB*1 

i  ^'4 

ANN*0,0 

1S4 

A0*0 . 0 

1S5 

KAANK*0 

1S6 

C 

1S7 

c*** 

FILL  AAAA'/S  NITH  BANO  BATA 

ISS 

c 

1$9 

KJ*1 

190 

DO  300  J0*IB6N,I£ND 

191 

EB<NJ)*FiJO>-DF 

192 

CB(N3)*C<ja> 

193 

RI</<J)*ALAS(JO) 

194 

B1<KJ)*XLSTA<JQ) 

195 

KJ*KJ*1 

196 

300 

CONTINUE 

197 

c 

1  9$ 

CALCULATE  ARRAYS  ’A'  AND  '1 

/O', 

c*** 

INITIAL  BAND  INTENSITY  ’AO 

200 

r 

C 


.% , 


,• 

r 


FORTROH  IK!  PROaRkH 


call  l1J3Lk»<tB,Cc!.M0A.M.AHH.B0  > 

CALL  SLCBi  DI  ,I1DA.R.A0.AMH,£B.CB,R1  ,BO,B> 
CALL  Hfn(B.MBA,».fir,2> 

CALL  BLDA<A,ltDA,R.AO,RI  ,£S,CB.AHH.BO> 
CALL  HF1T<A,IIDA.H.FIT,2) 

CALL  ASCAL£<A,MDA,M,N,C.CL> 


CALL  HF71(A,HDA,H,N,B,»nB,NB,7AU ,KRAHK,RH0RR.H,B,1P1 
AO=S<J)/SCL<J> 

BO-HCK/TO 

H*2 


B£G1H  2T£RAT:0N  TC  SOLUC  FOR  T£MP£RATLIP.£  ANB  INTEHSn'/ 


BO  3 SO  ILOP^l.SO 

CALL  MBLNN ( £B ,C3 , MBA . ANR .60  > 

CALL  BLBB<B1 ,RBA,H.A0.AhR,£3,CB.RI,B0.S) 
CALL  NFlT<B,MOA,tt,FIT,J> 

CALL  BLBA<  A,  MBA.  ft,  AO,  RI  ,£B.CB.AHI1. BO  > 
CALL  HFIT<A, MBA, N, FIT, 2> 


CALL  ASCAL£<A,MBA,M,H,SCL) 

CALL  HFTI<A.MBA,M,N,B,MDS,NB, TAU .LRANK ,RN0RM ,H ,G , I P > 


TEST  FOR  C0HV£R6£NC£  OF  TEMPERATURE 

Al*B(l)/SCL(2>i-A0 

B1^B<:2)/SCL(2>*B0 

TEMPO’HCK/BO 

T£MP]’HCK/B2 

A0=Al 

B0‘B1 

IF<ABS(  TEMPO-TEMPI > .L£. 0.5  >  GO  TO  jrO 

CONTINUE 

CONTINUE 

IF  (FRANK  .EO.MOt  GOTO  <f93 

CALL  NOISE(A. MBA ,M,N, RMORM . KRAMK . MO , SCL , HCK . 50 , CBA .SBT > 


MRITE  OUT  BATA  TO  FILE  UNIT 


MRITE(7  ,^kS)lLOP. AO, SBA. TEMPI  .STIT 
F0RMAT(IF.4E:5.<i> 

IF  (IBNB.ECi.2)G0T0  2000 
IF(IBNB.E0.3)G0TO  3000 
IF < IBNB .£0.4) GO  TO  4000 
IFdBNB.EO.S'GOTO  ‘='000 


FORTRON  :u  SVBROUTIHE 


FORTRAN  :<J  ■■■USROUTINE 


FILE  NRNE  '  SNSCL 


URITTEH  BV  RRRRIS  NEAL 

UTAH  STATE  UNIOERSITV  JUNE  26,19S4 

SUBROUTINE  SCALES  BATA  FOR  FFT  BAIN  AND  POST  ANR  BAIN 


23  100 

24  C 

25  C 

26  C 


SUBROUTINE  SNSCL (IF  FT . I  BN .ITL ,KLSTR) 
aiHENSIGH  XLSTR(c,0> 

SCALE  FOR  FFT  GAIN 
Si.  *0.0 

IF( IFFT.Eu.9>  SL-1.0 
IFdFFT.tO.O)  SL'O.S 
IFdFFT.EO.T)  SL~0.2S 
IFdFFT.£Ci.I0)SL=2.0 
DO  JOO  1*1  dTL 
XLSTRCI  )*>iiSTRd  )*SL 
CONTINUE 

SCALE  FOR  POST  AMP  BAIN 
SL*0.0 

IFdBN.EO.S)  SL*4.0 
IFd&N.E0.10)SL*2.0 
lFdSH.E0.20)SL*1.0 
IFdSN.EU.S0)SL*0.4 
IFdGN.EO.lOO>SL*Od 
DO  200  1*1  dTL 
XLSTRd  )*SLSTRd  >»SL 
CONTINUE 


RE  TURN 
END 


‘V. 


i 

\ 


i 


I 

I 

I 


I 


i 


'^J  I^J»J  f.  •'i  ».'WJ  I'i  PJ 


C  FOUrP.ftH  m  ZUSROUTINS 


2  C  F>-:RTRRH  IV  SUSROUTINi 

-  c 

J  C  FIL£  NRM£  •  HDLNM 

■t  C 

5  C  URITTEN  BY  PRRRJS  N£RL 

a  C  UTRH  STATE  UNIVERSITY  JUNE  2o.J9S4 

7  C 

$  C  SUBROUTINE  CALCULATES  THE  H0RHALIZ1N6  FACTOR  FOR  THE  SOLTZNAN 

->  C  TEMPERATURE  MODEL 

JO  C 

11  SUBROUTINE  MDLNMC £B ,CB .MDA .M , ANM ,B0) 

12  C 

13  DOUBLE  PRECISION  EB< MDA ) .CS( MDA ) .ANM . D£ XP 

14  ANM^^O.O 

15  DO  100  I‘1,M 

li,  ANM=ANM*(CS(  I  )*D£XP<-EB<  1  )»B0)  > 

17  100  continue 

IS  RETURN 

19  END 

20 


tl'. 


X 


210 


i 


i 


I 


r 

n 

»  I 


3 

9 

10 

11 

12 

13 

U 

15 


i  <:• 

19 

20 
21 


F0P.TRaH  IV  I-LIFROLiriliE 

FOPTPaH  IV  SUBPOUTIUE 

fill  NftNE  =  BLUB 

UPITTiN  BY  PAfi/tlS  NEf)L 
UTRH  STATE  UNIVERSITY 


100 


JUNE  £6,7v«3 


SUBROUTINE  BUILDS  THE  ARRAY  "B"  FOR  SOLTZHAN  MODEL 
LEAST  SQUARE  fITTINS  ROUTINE 

SUBROUTINE  BLDBiDI . MBA ,N ,A0 . ANN ,EB .CB , R I ,B0 , B > 

DIMENSION  RKMDA)  ,DHKDA  >  .S<MDA) 

DOUBLE  PRECISIOtt  EBi  NBA)  .CB'  MDA1  .DEUP 

DO  100  I-1,M 

,il10D=<A0*RI  (  I)  *CB(  I  )»<BE>iP<  -EB<  I)»B0>  >  )  /ANN 

B<  I>=DI(  I  >~>iM0D 

CONTINUE 

RE  TURN 

END 


'/yf 

••■I-.'. 

L. 


U 


r  t 


^  2m. .V  »*«  .\-‘r  .'li 


mtmrnmsmsm 


0  rx  «  o-  o  'N  >  'o  0  (v  0-.  o  »♦>  >  «•*»  o  lx 

•N  *^1  ''I  »•«  •N  »-5  W  ’••  4  Cl  ? I  Cl  C*!  Cl  C| 


ORTRaN  I{l  Si;?ROLITIN£ 

FORTRAN  I'J  SUBROUTIHE 

FILE  NRME  =  ASCRLE 

MRITTEH  BY  PARRIS  NEAL 

UTAH  STATE  UNIUERSITY  JUNE  2c,,19$4 

SUBROUTINE  NORMALIZES  THE  VECTORS  IN  THE  ARRAY 
THE  SCALING  FACTOR  IS  RETURNED  IN  ARRAY  SCH I > 

SUBROUTINE  ASCALE'.  A  .  NBA  .N,N.  SCL  > 

DIMENSION  AiMDA.tn  ,SCL(N> 

DOUBLE  PRECISION  SM.DSORT 

DO  100  I=1.N 
SN^O.O 

DO  2iO  J-i,M 
SM^SM*<AiJ,I)*A(J,.I)) 

CONTINUE 

SCL<I}^DSaRT(SN> 

CONTINUE 

SCALE  ARRAY  A  BY  THE  SCALE  FACTORS  SCL<U> 

DO  200  I^I.N 
DO  210  J^J.N 

IF<SCL< I KLT. 2 .OE-iO>  SO  TO  20S 
0(J,I)-AuT,I)/SCL(I> 

SO  TO  220 
A(J.I>^0.0 
CONTINUE 
CONTINUE 
RE  TURN 


213 


fc.T  •_”*  V-  •r 


fORTf^N  IV  PROGRfiM 


ki 


A 

C 

FORTSH  10  PR06P.4-H 

3 

L 

c 

Flit  UPM£  =  HFTI 

■f 

c 

5 

■3LIBP.0UT2N£  HFTI  <  P .  HOP  ,H  ,S  .  MOB  ,HB ,  TPU  ,KP.PHK  ,PMOPH  ,H  .5 

<y 

c 

7* 

c 

C.L.LPHSON  PHD  P.J.HPHSOH,  J£T  PP.OPULSIOH  LPBOP.PTO'P.V ,  1 

3 

c 

TO  PPP£PR  IN  ’■SOLOINS  L£PST  SOUPRF^S  PP.OBLBHS',  PREHTJCE 

9 

c 

SOLVE  LEPST  ■SaUPRES  PROBLEM  USIHG  PL60RITHM,  HFTI 

10 

c 

11 

DIMENSION  P<MDP,N> ,S(MDB,NB> ,H<N) .G<N) , RNORMi MB ) 

12 

INTEGER  IP(N1 

1 S 

DOUBLE  PRECISION  ■SM ,DZERO ,BBLE 

14 

SZER0=0. 

1 5 

DI£RO=O.DO 

1  o 

FPCT0R=0.001 

17 

c 

IS 

H^O 

19 

LDIPG^MINOi M,H1 

20 

IF  (LDIPG.LE.O)  go  to  270 

■21 

DO  SO  0=1,LDIPG 

22 

IF  (J.EO.l)  GO  TO  20 

2S 

c 

24 

c 

UPDPTE  SQUPRED  COLUMN  LENGTHS  PND  FIND  LMPK 

■25 

c 

26 

LMPX*J 

<77 

DO  10  L‘J,N 

2$ 

H<L  >*H<.L)-P<J- 1  ,L>*»2 

29 

IF  <H(L>,GT.H(LMPX>>  LMPX-L 

30 

10 

CONTINUE 

31 

IF (DIFF (HMPXfFPC TOR*H <  L MPX  > , HMPX > )  20 .20, 50 

32 

c 

S3 

c 

COMPUTE  SQUPRED  COLUMN  LENGTHS  PHD  FIND  LMPX 

34 

c 

3  5 

20 

LMPX=J 

Jo 

DO  40  L^J.N 

37 

HiL)-0. 

30 

DO  30  I=O.M 

39 

30 

H(L  >-H{L  )*P(  I  .L  >-»-»2 

40 

IF  (H<L> .6T.H<LMPK>)  LMPX^L 

41 

40 

CONTINUE 

42 

HMPK'HILMPX  > 

43 

c 

44 

c 

LMPX  HPS  BEEN  DETERMINED 

45 

c 

46 

c 

DO  COLUMN  INTERCHPNGES  IF  NEEDED 

47 

c 

40 

50 

CONTINUE 

V* 

IP<J1=LMPX 

50 

IF  <  :p<J).£O.J>  go  to  70 

JUN  li 
IL,  J9T 


LJ 


.’ %*  s 


W  ■i 


c 


F OPT AH  10  P POOR AH 


:oi 

SO  2a0  JS-J.HS 

202 

c 

JOS 

c 

SOLO£  THE  K  BV  K  TRIAHGUlAP  SVSTEH 

104 

c 

JOS 

BO  2 JO  L  =  J,K 

J06 

■SH=DZ£RO 

JOT 

1=KPJ-L 

JOS 

IF  (l.EO.KO  SO  TO  200 

J09 

1PJ’1*J 

JJO 

DO  J90  J^IPJ,K 

JJJ 

J90 

SH*SH*A< I .J  >»DBL£i > 

JJ2 

200 

SMJ^SH 

US 

2J0 

S< 1 .JB)^<B< I .J8>-SHJ )/A(l U  > 

1J4 

c 

US 

c 

COMPLETE  COMPUTATION  OF  SOLUTION  VECTOR 

2  J  o 

c 

JJT 

IF  (H.ECl.N)  60  TO  240 

US 

BO  220  J^KRJ.N 

JJ9 

220 

B(J ,JB)-SZ£RO 

J20 

DO  230  1=J,N 

J2J 

230 

CALL  HJ2  (2.1 ,KRJ ,N ,A(l U> .NBA .6t 1) ,B( J .JB) , J ,MDB,J> 

J22 

c 

J23 

c 

RE-ORDER  THE  SOLUTION  VECTOR  TO  COMPENSATE  FOR  THE 

J24 

c 

column  INTERCHANGES. 

J2S 

c 

J2o 

240 

DO  250  JJ*J.LD1A6 

J27 

J’LD1A64J-JS 

J2S 

IF  (  IP<JKEO.J>  60  TO  2S0 

J29 

L=lP(jy 

ISO 

TmP*B(L,JB> 

JSJ 

B{L,J8)=B(S.SB> 

1S2 

SiJ.JBT^TMP 

JSS 

2  SO 

CONTINUE 

JS4 

260 

continue 

JSS 

c 

1 3*^ 

c 

THE  SOLUTION  VECTORS.  X .  ARE  NON 

137 

c 

IN  THE  FIRST  N  RONS  OF  THE  ARRAY  5< , > . 

2  j  V 

c 

23 

2T0 

NRANK^K 

2  40 

RE  TURN 

4.  Uj 
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FOP.TP.RH  1<J  PR03KPH 


C 
C 
C 

5  C 

6  C 

7  C 
i  C 
'>  C 

20  C 
22  C 
22  C 
2  S  C 
2  4  C 
2  5  C 

4  i2>  I* 

27  C 
2S  C 
2'>  C 
20  C 
22  C 
22  C 
2S  C 

24  C 

25  C 

26  C 

27  C 
.0  C 

29  C 

30  C 
32  C 

32  C 

33 

34 

35 

JO 

37  C 

31 

39 

40 

42  C 

42 

43  20 

44 

45  20 

4/3 

47 

42  30 

49  C 

50 


FOP.TRPH  20  PPOOPPM 
FILE  HPHE  '  HI2 

SUBROUTINE  H22  i NODE . LPIUOT .L2 .M .0 , lUE .UP . C . ICE , I CJ . NCO i 
C.L.LRNSOH  PHD  P..J.hRNSOH.  JET  PROPULSION  lRBORRTORY  .  .  *ri  JON  2. 
TO  APPEAR  IN  ’SOLUIH6  LEAST  SQUARES  PROSLENS' .  PREHT ICE- UAL.. ,  2  ■‘I 

CONSTRUCTION  AND/OR  APPLICATION  OF  A  SIN6ELE 
HOUSEHOLDER  TRANSFORMATION. .  <J»  I  *■  U«Ul**TWB 

MODE  OF.  2  TO  SELECT  ALSORITHM  H2  OR  H2 

iPIOOT  IS  THE  INDEX  OF  THE  PIVOT  ELEMENT. 

L2,  M  IF  L2  .l£.  N  the  TRANSFORMAT  I  ON  HILL  RE  C0NST7VCTED  T '/ 

ZERO  ELEMENTS  INDEXED  FROM  LI  THFOUGH  M.  IF  LI  .OT.  N 
THE  subroutine  DOES  AN  INDENTITY  TRANSFORMATION. 

U( > .lUE .UP  ON  ENTRY  TO  H2  U< >  CONTAINS  THE  »IO0T  VECTOR. 

lUE  IS  THE  STORAGE  INCREMENT  SETNEEN  ElEMEN'^S. 

ON  EXIT  FROM  HI  U<  >  AND  UP 
CONTAIN  QUANTITIES  DEFINING  THE  VECTOR  U  OF  THE 
householder  transformation.  on  entry  to  H2  Uf  ■ 

AND  UP  SHOULD  CONTAIN  QUANTITIES  PREVIOUSLY  COMPUTE 
BY  HI.  THESE  MILL  NOT  BE  MODFIED  BY  M2. 

CO  OH  ENTRY  TO  HI  OR  H2  CO  CONTAINS  A  MATRIX  MMICH  MILL  BE 

REGARDED  AS  A  SET  Of  VECTORS  TO  MMICH  THE  HOUSEHOLDER 
TRANSFORMATION  IS  TO  BE  APPLIED.  ON  EXIT  CO  CONTAINS  TE 
SET  OF  TRANSFORMED  VECTORS. 

ICE  STORAGE  INCREMENT  BETMEEN  ELEMENTS  OF  VECTORS  IN  C(> 

ICV  STORAGE  INCREMENT  BETMEEN  VECTORS  IN  C-.  >. 

NCV  NUMBER  OF  VECTORS  IN  CO  TO  BE  TRANSFORMED.  IF  NCV  .LE. 
NO  OPERATIONS  MILL  BE  DONE  ON  CO. 

SUBROUTINE  H12  <  MODE . iPIVOT ,L 2 . M .U , I UE ,UP .C . ICE . ICV . NC V > 

DIMENSION  U< lUE .M) ,C<  2  > 

DOUBLE  PRECISION  SM.D.DBLE 
ONE^i. 


IF  <O.GE .LPIVOT .OR.LPIVOT.GE .LI .0R.L2 .GT .M>  RETURN 
CL^ABS  <  U (  I  ,LPIVO T  >  .> 

IF  {MODE. £0.2)  GO  TO  <»0 

******  CONSTRUCT  THE  TRANSFORMATION  ****** 

DO  10  J*LI.M 

CL*AMAXl(ABSiU( 1 .J>  > .CL  > 

IF  (CL)  130. ISO. 20 
CL1NV*QHE/CL 

SM*<DBL£(U(2  ,lPIVOTO*CLIHV>**2 
DO  30  J‘L2,M 

SM*SM* (DBLE(U( 2 .J > )*CL INV >**2 

CONVERT  DBLE.  PREC.  SM  TO  SNGL .  PREC.  CM 2 


/ 


i 


fi'iRTP.fth  7U  PP.06PAM 


CL=CL*^-‘3PTi  SPl  ) 

:f  <m  2  .lp:i.iot>  >  zo.so.40 
CL*-Cl 

UP’Ui  2  .LPI</OT>-CL 
LH2.LP:OOT'fCL 
00  TO  ro 

PPPLT  7H£  TP»P2F0PltPri0H  I*U*iU«*T1/B  TO 

IF  <CL>  130.230.70 
IF  (NCW.if.C)  PETUPH 
BABBLE  <  UP  )»U<  1  .LPIum> 

B  PUST  BE  POhPOSTTlUE  HEPE .  IF  5*0, 


IF  (B>  OO. 230. 130 
S^OHE/S 

I2‘2-ICU*ICE*t LPIUOT- 2 > 
IHCP^ICE»(.L1-LPIU0T> 

00  120  J-l.NCU 
I2*I2*ICU 
I3*I2*INCP 
14*13 

SP*C(I2)»BBLE<UP) 

DO  90  l*Lt.H 

SP*SP4C (13^ •DBL £<U<1,I>> 

I3*I3*ICE 

CONTINUE 

IF  (SN)  100.120,100 

SN*SN»n 

C(12)•C(I2)^SN»0BLE^UPy 
00  110  I*L1.M 

C< 14  i*C< I4>*SN*DSL£<U(  1 ,1 i ) 
14*14*IC£ 

CONTINUE 

CONTINUE 

PETUPN 

FNn 


C 


PE TUPN 


J  I?  0»  U| 
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F-iRTRUH  IV  FUNCTION 

FORTRf)N  IV  FUNCTION 
FILE  NNNE  =  OIFF 
FUNCTION  dIFFOi,'/) 

C.L.  LANSON  AND  fl.J.  HANSON,  JET  PROPULSION  LABRATORY,  7  JUN  297i 
TO  APPEAR  IN  ’SOLVINB  LEAST  SQUARES  PROBLEMS',  PRENTICE-HALL  JRT'i 

DlFF*X-y 

RETURN 

END 


f  -  o  C  l 


220 


fORTRUH  lU  ZUBPOLITINE 
4  C  fCRTP.AN  10  SUSP.OUTItte 

r 

3  C  FILE  NAME  -  NOISE 

4  C 

5  C  NPITTEN  BY  PPRP.1S  HERE 

6  C  UTRH  STATE  UNIOERSITY  JUNE  29,  19S4 

7"  C 

S  C  SUBROUTINE  TAKES  RESIVUE  FROM  LEAST  SQUARE  FITTINS  ROUTINE 

9  C  -HFTl’  AND  CALCULATES  THE  STANDARD  DEUIOATION  OF  THE 

JO  C  BAND  INTENSITY  "A"  AND  THE  TEMPERATURE  "T” 

JJ  C 

12  SUBROUTINE  NO ISE< A .HDA .N .N , RNORN .KRANK ,MQ , SCL ,HCK ,B0 .SDR . SBT ) 

IS  C 

14  DIMENSION  A< MDA .N > .SCL < N > 

15  C 

lo  OARN^ ( RN0PM**2 > / <  NO- KRANK ) 

17  0NA‘( < { A<2 •2)**2)*<A( 1 ,2>»*2) > / <  (A< 1 , 1 )*A(2 .2> )**2>  >*UARN 

10  SBA^(SaRT<ONA>>/SCL(l> 

SDB‘(S0RT<<JARN/<A<2,2)*»2>)>/SCL(S) 

SDT^< SDB»HCK )/< B0«*2 ) 

RE  TURN 
END 


i 


1 

E 

1 

c 

FOP.TROH  IV  PPOSKPH 

r 

2 

C 

FOP.TP.fiH  IV  PPOePPH 

S 

c 

C 

FIL£  NAH£  -  LIHPHP 

1 

4 

c 

5 

c 

MPITT£N  BV  PPPRIS  N£OL 

6 

c 

UTAH  STPT£  UNIV£PSITy  HP'/  23,1904 

/ 

c 

THIS  PP.UGPPH  P£finS  FFT  FIL£S  FROM  TPP£ 

y 

9 

c 

DPTP  IS  CONVEP.TEV  TO  FORTRPH  IV  FOP.HRT 

►  j 

70 

c 

HEPDEP.  IHFO  IS  PERU 

7  7 

c 

TAPE  HUHBEP. 

1 

72 

c 

DPTP  FP.PHE  HUHSER 

1 

7J 

c 

FFT  EXPONENT 

24 

c 

P/O  SPIN 

75 

c 

DPTE  <Y£PR.DPY.H0UR,H1H,S£C> 

76 

c 

FFT  IS  P.EPD.  PHPSE  CORRECTED,  PHD  LINE  PHPLITUDES  FOUND 

77 

c 

THE  SPNPLE  HUMSEP.  NHEP.E  THE  LINE  IS  EXPECTED  TO  BE 

75 

c 

IS  REPD  FRON  P  DISH  PILE  PND  P  SEPRCN  IS  USED  TO  FIND 

7^- 

c 

THE  PCTUPL  PEPK  VPLUE  USING  P  NPNMIN6  PPODIJPTION  ROUTIN 

1 

20 

c 

1 

22 

c 

SUBROUTINES  USED  PP.E 

22 

c 

TPPE.TPPEC  -  TPPE  DRIVE  UTILITIES 

» 

23 

c 

SFCOHV  -  CONVERT  12  SIT  TO  FOP.TRPN  IV  36  BIT 

24 

FFTP.H,P.KOSC  -  DISK  DRIVE  UTILITIES 

25 

c 

PHPSEC  -  PHPSE  CORRECTION  ROUTINE 

i’ 

26 

c 

PPODZ  -  HPHHING  PPODI ZPTION/ INTERPOLPTION 

.•>7 

c 

DTINE  -  CONVERTS  HEPDEP.  INFO  TO  DPTE/TINE 

1 

2  V 

c 

UNIT  -  INPUTS  FRPNE  HUHBERS  TO  BE  PROCESSED 

1 

2 

c 

SEPP.CH  -  SEPP.CHES  FOR  PEPK  PHPLITUDE  OF  LINES 

JO 

c 

REOLPT  -  LOPDS  MEMORY  NITH  P.EPL  PPRT  OF  FFT 

J7 

c 

32 

COMMON/ P.KOSC /BUFO  , BUFF  EP.<  2047  i  ,HBUFO  ,HBUFF  <  2047  )  ,XRING<  1024)  , 

33 

1 

YRING<I024) ,IDPTP 

■. 

34 

DIMENSION  IPPSS<  5) ,INFRPM<30> , IHDTP(IO) , IGPINdSi , IPOS(GO> 

35 

DOUBLE  PRECISION  NPVNFGO.i 

M 

36 

c 

1 

37 

c**»* 

REPD  IN  INITIPL  VPLUES 

r 

3S 

c 

39 

ITFMS^O 

40 

NGRPS‘0 

42 

CPLL  L  INIKNGRrS  .  INFRPM ,  IGPIN ,  I TFMS ) 

42 

c 

. 

43 

c 

NRITE  TO  OUTPUT  FILE  TOTPL  *  OF  FP.PMES  PROCESSED 

L 

44 

c 

1 

45 

MP.ITEio,t:,06>  ITFMS 

1 

4t^ 

c 

1 

c**** 

MINDON  UIDTHS  FOR  <*  SPMPLEOi  PHPSE  CORRECTION  ROUTINE 

\ 

w 

IPPSSy I) -3 

E 

50 

!PPSS<2>'S 
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PROGPUM 

IPPXSi.  j  )•»“ 
IPaC.S(  4  >^27 
IPPSSi  5  >-»o5 
NSZ=22 
NPASS-5 
inorP‘1 


5S  C»»»* 

59  C 


6/ 

6S 

69  C 

70  C**** 
72  C 

72 

73 

74  262 


7 3  607 

79  270 


C 

c*»** 

c 


INniPLlZE  FPT  INPUT  TAPE  OPIUE 

NP1T£<4, 260> 

FOPNATC  ’> 

MPl  TEC  4,200') 

FOPNATC"  INPUT  TAPF.  SP.IUE  NUHSEP  ".»> 

P.EAD(4.220>:2P.U 

FOPNPTC 14) 

2F  slUP.U.EO.l-inP.U-O 
10PP=4 

call  TAP£<IT)PU.2,SUFO.O) 

PEAD2N&  IN  LINE  POSITIONS 

NP1T£(4,260) 

MP2T£(4,262.) 

FOPNATC’  PEADINB  IN  LINE  POSITIONS’) 

PEAVC  5,606)ITL 

00  270  1^2,ITL 

P.£AD<S ,  607  >  IPOS  CD  .  NAUN  ( I  '■ 

.  FOPNATC  16,825.6.) 

CONTINUE 

SPACE  ANEAO  TO  FIRST  TAPE  NEAOER 
CALL  TAPE<10RU,4.SUF0,2) 

SPACE  TO  RECORD  SPECIFIED  BY  INFRANiD 


90 

92  220 

92  C 

93  C**** 

94  C 


IC0UNT=1NFRAH< 1 )- 1 

IF  CICOUNT .EO.O)  60  TO  220 

IC0UNT~IC0UNT*2 

CALL  TAPE<IDRU,4,BUF0,IC0UNT) 

CONTINUE 

START  READING  DATA 
JX^2 

DO  2000  IJ=2.N6RPS 

DO  2000  IFM^INFRANCJZ) .INFRANCOI*!) 


f^h**.*******.****^******  ******* 


*  •  V  V-  •  .  "  •  '>^J.  --i.  • -V-  • 


Il 

- 

f0P.7Rf>H  10  PPriGRUM 

1 

.'0/ 

c 

PSPDING  HERDER  INFO 

:o2 

c 

p 

20  J 

r 

?iRIN6<30)  =  TRPE  # 

:04 

c 

>l/ilN6C24>^DA7f)  fP.RME  6 

1 

205 

c 

KRIN6(14>=FET  EXPONENT 

/0<S 

c 

XRIN6<5)  =R/D  GRIN 

107 

c 

XR1NG<6-  ^•DATE/TIME 

:os 

c 

i 

109 

1 

110 

c 

I 

111 

CALL  TAP£( IDR0.7,XRING.2S6) 

112 

CALL.  SPC0N0(XRING.XRING,256> 

11  j 

IHDTA(  1  .>-XR!NG(30y 

114 

lHDTA<2)-XRlNG(24y 

115 

1HDTA<0)=XRING< 14> 

12*^ 

IHDTA<9)=XRING( 5> 

- 

117 

1HDJA< lOl’IGAtNL IJ > 

1 

22$ 

CALL  DT1M£(XR1NG(6> .YEAR, DAY ,H0UR,N1N ,S£C) 

I 

119 

IHDTAi 3}‘Y£AR 

J 

ISO 

IHDTA(4)=DAY 

ISl 

1HDTA<  5)=H0UR 

( 

ISS 

IHDTA<6)‘N1H 

■ 

123 

1HDTA<7)’S£C 

124 

NRIT.£<3,525)IHDTA<  !.■>  .IHBTACl) 

■ 

12S 

525 

FORHATC  ACTUALLY  READ  TAPE  «  ’,16,’  FRANE  0  ’,16) 

1 

126 

C 

1 

127 

c**»* 

READING  IN  FFT  RECORDS  5,6  UHICN  CONTAIN  ALL  LINES  FOR 

", 

120 

c 

4-2, 5-1 , $-5,7-4  ON  RANDS 

129 

c 

ISO 

MRITE<4.160) 

232 

URITE(4.520>IFN 

V 

132 

520 

FORNATC  READING  IN  TRANSFORH  “.II,"  FROM  TAPE") 

‘1 

233 

C 

1 

134 

SPACE  TO  RECORD  5 

1 

135 

c 

236 

CALL  TAPE< IDR0.5.BUF0.4) 

2  3T 

c 

13S 

c**** 

READ  RECORD  5 

139 

c 

140 

CALL  TAPE(  IDRO  ,7  ,RUF0,204$-> 

141 

CALL  SPC0NO(BUFO,BUFO,2O49> 

142 

ICYLD‘0 

1 

143 

IERR*0 

7 

144 

CALL  FFTRNL IDATA,ICYLD,1 .BUFO.IERR) 

145 

c 

146 

J#### 

READ  RECORD  6 

" 

147 

c 

* 

140 

CALL  TAPE(IDR0.7.BUF0.204$) 

i4'^> 

c 

150 

CALL  SPC0H0<.  BUF0,BUF0.204$  ) 

224 


FORTRf>H  lU  PROGRRH 


:5i 

:CVLI!-i 

jsi: 

CALL  FP  TRK  (  lOFirR  .  ICVLO  ,  2  ,  BUFO .  lERR  > 

JSS 

C 

254 

c**** 

REROy  TO  PHASE  CORRECT  DATA 

2  55 

c 

2  5^ 

HR2TE<4..260) 

257 

HR2TE(4.S40)1FH 

25$ 

540 

FORMAT<"  PHASE  CORRECTING  FFT  FRAME  4  ".I4> 

2  5$ 

CALL  PHASECiNSZ ,HPASS.IPASS> 

260 

NR2T£( 4.550) 

262 

550 

FORNRTf"  PHASE  CORRECTION  COMPLETE") 

262 

C 

263 

c**** 

URITE  OUT  HERDER  DATA  TO  LINE  AMPLITUDE  FILE 

264 

c 

2  65 

BO  605  10^2.20 

26*3 

UR1TE<6.60'3>IHDTA*  10.) 

267 

606 

F0RMAT<I6) 

26$ 

605 

CONTINUE 

269 

C 

270 

C*»** 

READ  PHASE  CORRECTED  DATA  IN  FOR  LINE  AMPLITUDE  SEARCH 

272 

c 

272 

CALL  REALPT 

273 

c 

274 

c»**» 

READ  IN  HAMMING  TABLE  VALUES  INTO  COMMON  AREA  ’BUFFER' 

27  5 

c 

276 

RE  AD (7) BUFO . BUFFER 

277 

REMIND  7 

27$ 

c 

279 

c**** 

START  TO  FIND  LINE  POSITIONS  AND  AMPLITUDES 

2$0 

c 

2S2 

MRITE(4,$50)IFM 

2S2 

$50 

FORMATC  SEARCHING  FRAME  <*  ’.14) 

2S3 

DO  900  IT=2.ITL 

2S4 

IPSTcfIPOSL  IT) 

2  $5 

XOFS^O.O 

2  $6 

VLFS-0.0 

2  $7 

CAL  L  SEARCH  ( I  PST.  AOFS  ..ULFS  > 

2$S 

C 

2S9 

c»*** 

MRITE  SAMPLE  NUMBER  AND  AMPLITUDE  TO  DISK 

2  90 

c 

292 

MRI TEL  6. 7$S)X0FS,VLFS 

2  92 

7$5 

FORMATLF  20.2,E25.6) 

293 

900 

CONTINUE 

294 

C 

295 

c»*** 

SPACE  TO  NEXT  HEADER 

296 

c 

297 

CALL  TAPELIDRO .4,BUF0,2) 

2  9$ 

c 

1  Cl  Cl 

2000 

CONTINUE 

200 

C 

FORTPdN  IV  PROePAM 


JO  I 

c**** 

SPACE  TO  NEXT  GROUP 

J02 

c 

JOS 

IF( 1J.£0.HGPPS>  GO  TO  2000 

204 

JZ=JZ*2 

205 

ICOUNT=< < lNfPAM(JZ>-IHFPaH<J. 

206 

CALL  TAPE <IDPV . 4 . BUfO . ICOUN T 

207 

c 

200 

2000 

CONTINUE 

209 

c 

210 

c»*** 

RENINO  TAPE 

21 1 

c 

212 

CALL  TAPE( IDRV .J.BUPO.O^ 

213 

c 

r»  rjc-j  r>r>nr>c><^c>rir^r>r»rjrjr»r>r>r>fjr»c^r>c^r>r>r>t*»r>r>r>r-> 
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r-ORTP.f>N  /(/  PHRSi  CORRECTION  CUBP.OUTIHE 

FORTRAN  lil  PHASE  CORRECTION  SUBROUTINE 


RON  BROUN 
EE  DATA  FACILITY 
L06AN,  UTAH  S4S2S 

THIS  SUBROUTINE  TAKES  A  FFT  FILE  FRON  A  BISK  DRIOE  AND 
PERFORMS  A  PHASE  CORRECTION  BY  MEANS  OF  SUCESSIUE 
AOERASINS  AND  PHASOR  MULTIPLICATION. 

A  MAXIMUM  OF  FINE  PASSES  ARE  ALLOUED  FOR  THE  AUERASING. 

EACH  PASS  AOERASES  A  SIUEN  NUMBER  OF  DATA  POINTS.  IT 
IS  IMPOTRANT  THAT  THE  TOTAL  NUMBER  OF  POINTS  A0EP.A6ED 
<0N  ONE  SIDE  OF  THE  POINT  OF  INTEREST)  OF  ALL  PASSES 
DOES  NOT  EXCEED  IOCS.  IF  SO.  THE  SUBROUTINE  HILL  RETURN 
UITHOUT  PERFORMING  THE  CORRECTION. 

CALLING  LIST:  PHASEC ( N . NPASS . IPASS ) 

UHEREiN  IS  THE  LOG  BASE  S'  OF  THE  HUMBER  OF  DATA  POINTS 

NPASS  IS  THE  NUMBER  OF  PASSES  AND  CANNOT  EXCEED  FI'JE 
IPASS  IS  AN  ARRAY  OF  THE  NUMBER  OF  DATA  POINTS  TO  BE 
AVERAGED  ON  THE  NTH  PASS.  THE  SUMATION  FROM  1  TO 
NPASS  OF  S»IPASS<N)*I  CANNOT  EXCEED  !0i4.  IF  SO,  THE 
ROUTINE  HILL  RETURN  UITHOUT  OPPERATION. 

NOT  NOTED  IN  THE  CALLING  LIST.  BUT  EQUALLY  IMPORTANT.  IS  A  COMMOF 
SLOCK  NAMED  RKOSC  UITH  AT  LEAST  6144  FREE  LOCATIONS  UITH  THE 
6I4STH  LOCATION  INDICATING  THE  DISK  DRIVE  TO  BE  USED. 

REVISION  HISTORY! 

ORIGINAL  VERSION  FON  BROUN  i  JUL  SC 

SUBROUTINE  PHASEC ( N . NPASS . I  PAS ) 

COMMON / RKOSC /BUFO .BUF  F  ERi 2040 > .HBUF( 1047) .XRING( 1024 > . 

I  YRING(  1024'-  .IDATA 
DIMENSION  IPASU) 

DOUBLE  PRECISION  XDOUB I .  XD0UB2  .  XDOUBC .  .^B0UB4 .  .iDOUBS 
I  YDOUBS  .  YD0UB2  .  YD0UB3  .  '/B0US4 .  /D0UB5 

SET  NUMBER  OF  RECORDS  IN  FFT 
NREC’N-10 

IF (NREC.LE.OIRE TURN 
NREC*2»*NREC 

NPAS^HPASS 
IF ( NPAS.L T.J) NPAS*I 
IF  <NPAS.GT .S)NPAS*5 
JEND^NREC/2 

JENV^i 2*JEND*2-HREC)»2040 
2  F0RMAT(X,SFB.2) 

JEND  POINTS  TO  THE  LAST  POINT  IF  ON  THE  LAST  RECORD. 


Rf>H  PURSE  COP.P.ECriOM  SUBROUTINE 


PSS16N  SCO  TO  NEXTl 
aSSlGN  SCO  TO  NEKTS 
RESIGN  SCO  TO  NEXT3 
RSSIGN  SCO  TO  NEXT -I 
IDir^O 
IS0=1 

GOTO  <5O.-tO.S0,2O.i0)  ,NPRS 

N1TTS=IPRS( S) 

IDIE^ZDIF^NITTS 

ISS^ISO 

15-ISO 

IP5-1S0*N1TT5*2 

ISC=lP5-*-NZTT5 

XDOUBS^O.ODO 

vnouBs^o.czic 

ASSIGN  500  TO  NEXT 4 

NZTT4=IPAS( 4) 

II>IF-‘IT}IF*NITT4 

IS4-ISQ 

14*130 

ZP4*IS0*NITT4*2 

ZS0*IP4*NITT4 

XD0UB4*0.0D0 

YD0UB4*0.0D0 

ASSIGN  400  TO  NEXTS 

NITT3*IPAS<S) 

IVIF  =  II1IF*NITT3 

IS3*IS0 

13*130 

IP3*IS0*N1TT3*I 

IS0‘IP3*N:TT3 

Xi}0UB3*0.0V0 

yD0UB3*0.0D0 

ASSIGN  300  TO  NEXTS 

NITT2*IPAS(2) 

IVIF*IBIF*NITT2 

IS2*IS0 

I2*IS0 

IP2*IS0*NITT2*i 

IS0*IP2*NITT2 

XDOVB2*O.ODO 

YI>0UB2*0,0D0 

ASSIGN  200  TO  NEXTl 

NITT2*IPAS< 1) 

IDIF*IDIF*NITT1 

ISI*ZSO 

II*ISO 

IPt*ISO*NITT 1*1 
IS0=1P1*N1TT1 


c 


fORTReiN  I'J  PHASE  CORRECTION  SUBROUTINE 


iOl 

XI10UB1=O.ODO 

102 

VDOUBl^O.ODO 

1  Oi 

IF ( ISO. GT .1024) RETURN 

104 

C 

lOS 

C 

I-Sl  THRU  ISS  ARE  POINTERS  TO  THE  FIRST  LOCATION  OF  THEIR 

106, 

C 

RESPECT lOE  RING  BUFFERS.  ISO  IS  USED  AS  A  LIMIT  ONL'/. 

107 

C 

11  THRU  15  are  pointers  THAT  NRAP  AROUND  IN  THEIR  BUFFERS 

108 

C 

A  PATH  HAS  BEEN  ESTABLISHED  USING  ASSIGNED  GOTO  STATEMENTS, 

109 

c 

THEREBY  RETAINING  GENERALITY  HHILE  INCREASING  SPEED. 

110 

c 

111 

IBUF=0 

112 

c 

112 

c 

IBUF  IS  THE  SUBSCRIPT  OF  BUFFER  SUCH  THAT  HBUF(IBUF) 

114 

c 

IS  THE  FIRST  LOCATION  TO  BE  LOADED  TO  OR  STORED  FROM  THE  BUFFER. 

US 

r 

ll-i 

IPTR=2048 

117 

JPTR=IPTR-2*IDIF 

US 

c 

119 

c 

DATA  ALNAYS  COMES  OUT  OF  THE  BUFFER  AT  IPTR  AND  IS  STORED  AT  JPTf 

120 

c 

UHICH  IS  ALNAYS  OFFSET  FROM  IPTR  BY  THE  DIFFERENCE  IDIF , 

121 

c 

122 

KHT^O 

122 

DO  60  I^l.ISO-1 

124 

XRING(  I  )-*0. 

125 

VRING<I)‘0. 

126 

60 

CONTINUE 

127 

DO  70  1-1024,2047 

128 

BUFFER(1>*0. 

129 

70 

CONTINUE 

120 

C 

121 

IERR=0 

122 

SO 

CALL  FFTRfUIDATA.HNT .0,HBUF< IBUF  > ,I£RR) 

122 

DO  85  1=1BUF*205C ,1BUF*4094 , 4 

1 2  4 

BUFFERL I )--BUFF£R(  I > 

125 

BUFFEP.<1*1)^-BUFFEP.(1*1  > 

12<y 

8  5 

CONTINUE 

1 27 

90 

KNT^KNT*! 

128 

IBUF 2048- IBUF 

129 

100 

K*SUFFER( IPTR) 

140 

Y^BUFFEPA  IPTR*!) 

141 

XDOUBl^XDOUBUK-XRINGKIl) 

142 

YD0UB1’YD0UB1*Y- YRINGLII) 

142 

XR1NG( 11 )’X 

144 

YRIN6<11)‘Y 

145 

X^XDOUBl 

146 

/•YDOUBl 

147 

XMA6‘rSQRT ( X»X*Y*Y  > 

148 

JF(XMAG.N£.0.0>G0T0  150 

149 

X‘1.0 

150 

XMAG=X 

v  .i*  w  ’Jij  Vo  'j  -j  -§  -f  ■>*  --j  •■  -•!  4  -I  V  w  V.  V  V  v-  w'  it'  V-  %/■  VI  VI  VI  s.n  vn  v»  vn  v*i  vt 

C*i  *«*  o  C  ■»  '^1  G  C.h  <  ►O'  O  C'  CH  O  0  Oi  “I  c  U  4,  <»J  K''  Ci  "O  <«.  N  Cv  *.•»  4*  <*i  K  ^ 
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FORTRAN  10  PHASE  CORRECTION  SUBROUTINE 


ISO 


200 


250 


jOO 


3  SO 

tP4 
'95 
’.9tS 
■97 
t9S 
‘  «.*/  «.■/ 

I'OO 


V-X/XMAfi 
V '/=■// KNAG 
XO^KRINGdPi  ) 
V0=VRIHG<1P2> 
)i=F)\*FO*V'/»YO 

'/•»xx#yo-yy»xo 

iPi‘iPi*i 

IF(IP1.S£.IS0>IP1»ISJ 

Il^Il-^l 

IFdl.GE.ISOIIl’ISl 
GOTO  NEXT! 

XI}0UB-2*XD0UB-i*X-XRING<I-2y 
YB0UB2  =  VT}0UB-2*V-YRING(I2) 
XRINGC  I  -2>  =  X 
VRINGf  1-2)  =  '/ 

X^XDOUB-2 

y^'/D0US2 

XNAG^SQRT (  X*X*'/»V  ^ 

IF(XNA3.N£.0.0 IGOTO  250 

X^l.O 

XNAG‘X 

XX^X/XHAG 

VY^Y/XHAG 

X0=XRIN6(1P2> 

Y0*YRIH6< IPS) 
X*XX»XO*YY»YO 
Y^XX»YO-YY*XO 
IP2*IP2*1 

IF(IP2.G£.ISI) IP2‘IS2 
12-^12*/ 

IF<12.G£.1S1) I2=IS2 
GOTO  NEXT 2 

XD0UB3-XD0UB3*X-XRING< 13) 
VB0UB3=YV0UB3*Y- YRING<I3) 
XRING< 13 )-X 
YRING< I3)^Y 

x-xnouB3 

Y=VTJ0UB3 

XNAG^SOP.T  (  X*X*Y*V  •> 

1F<XHAG.NE.0.0)G0T0  350 

X^J.O 

XHAG^X 

XX*X/XNAG 

YY*Y/XHAG 

XO^XRING< 1P3) 

Y0*YRING(IP3) 

X*XX«XO*YV*yO 

Y^XX*YO- YY»XO 

IP3=IP3*1 

IF  dP3  .G£.IS2)IP3=‘IS3 


's.  r*i  >  >0  K  <0  »•  O  'N  lo  >  «  'O  Is  <0  Os  O  «•  f«j  >  «  >0  K  <0  O-  O  C<  lo  >  »•’>  ■?)  Is  <0  Os  o  m  >  «  >0  Iv  «&  o 

O  O  O  O  O  O  O  O  O  •'I  ’N  C-i  c#  c<  Ci  »o  »*i  »•>  »n  rn  »o  »*»  »n  >  >  V  >  >  >  >  >  "T  *»• 

'M  C|  C4  Cl  N  04  M  Cl  Cl  C|  C|  C|  C-l  C^  C|  C|  C|  Cs|  .Ci  Cl  C|  C|  Cl  Cl  C|  C|  C^  C|  Cl  C|  C|  Cl  C|  C|  C|  C|  Cl  C|  C|  C|  C|  C|  C|  C| 


FORTKflH  JIJ  PHASE  CORRECTION  SUBROUTINE 


I3‘IS*i 

IF( lS.aE.IS2>lS-ISS 
SOTO  HEXTS 

400  XBOUS4^XnOUS4*X-XRIH6< 14> 

yV0UB4‘yV0UB4*Y-yRIN6< I4> 

XRIHG( I4>*X 

yRiHG<  n>‘y 
X^XB0UB4 
y=yj}0UB4 

XHR6^SQRT<X»X*y*y > 
lf<XNR6.NE.0.0)60TO  450 
X^l.O 
XHRS^X 

450  XX^X/XMRS 

yy^y/xnas 

X0=XR1NS<1P4) 
yo-yRiNS< IP4) 
x^xx*xo*yy*vo 
y*xx*yo-yy»xo 

IP4^IP4*i 

IF<lP4.SE.tS3)lP4’lS4 
14=14*i 

IF(14.S£.IS3U4^1S4 
SOTO  NEXT 4 

500  xn0US5^XD0UB5*X-XRlNS<JS> 

yi)ouB5*yBouBs*y-  ypiHSdS) 

XRIN6(I5)*X 
yRiN6(is)*y 
X*XI}0US5 
y=rynous5 

XNOS’SaRT(X*X*y*Y.> 

IF<XNR6.NL.O.O)GOTO  550 
X*1.0 
XHRG‘X 

550  XX^X/XMRG 

yy=y/xnRG 
X0=XP.ING<  IPS) 
yO’ypiNG< IPS) 

x^xx»xo*-yy*yo 
Y-xx*yo-yy*xo 

IP5*lP5-fI 

IF< IP5.SE.XS4)IP5'IS5 
15*15*1 

1F<15.GE.1S4)15*IS5 
C 

C  AT  THIS  POINT.  THE  SUMS  OF  ALL  THE  PASSES  APE  IN  X  AND  Y 

C  AND  THESE  POINTS  COPPISPOND  TO  NHAT  IS  IN  THE  BUFFER 

C  AT  JPTP. 

C 

C50  $00 


BUFF£P<JPTP)*X 


CH  C*> 


FoetTR^N  HI  PHRSE  CORRECTION  SUBROUTINE 


BUFFER':  JPTR*J)^T 
lPTR=IPTR*-2 

IF< IPTR.GT. 4095)IPTR-0 
JPTR=JPrR-i-2 

IFiJPTR.BT. 409S)JPTR*0 

IFdPTR.NE.O.RND.  IPTR.NE  .2040  .RNB .  <  JPTR  .  HE .  JEND  .  OR 
KNT.LE.NREC) >6070100 
IF(KNT.LE.l>6OT0  $20 

CRLL  FF TRK < lORTR. UN T-2,1, HBUF ( ISUF > , I ERR > 
IF<KNT.LT.NR£C>60T0  SO 
IF(KNT.GT.NR£C>RE TURN 
DO  S40  I-1PTR.IPTR*I02S 
SUFF£R<  2 )*0. 

C0NT1NU£ 


fORTf>/)N  lu  time  cove  decode  SUSROUTIliE 

EOPTPfiN  IV  TIME  CODE  DECODE  SUDPOUTINE 
FILE  NOME  IS  'DTIME' 

SENE  NONE 

UTOH  STOTE  UNIVERSITY 
LOSON  UTOH  04322 
JULY  2,  1979 

DTIME  DECODES  TIME  FROM  4  NORDS  REOD  FROM 
TIME  CODE.  ORROY  N  CONTOIHS  THIS  DOTO  IN 
ORDER  OF  REOD. 

SUBROUTINE  DTIM£(N, YEOR .DOY .HOUR .MIN .SEC  > 


"N"  IS  yhE  ORROY  OF  NORDS  TO  BE 

DECODED.  N  S.NOULD  BE  DIMENSIONED 
OT  LEOST  4.  Nil>  SHOULD  BE  THE  FIRST 
VOLUE  REOD  FROM  THE  TIME  CODE. 

"YE OR"  IS  THE  YEOR 

"DRY"  IS  THE  NUMBER  OF  DOYS 

“HOUR"  IS  THE  NUMBER  OF  HOURS 

"MIN"  IS  THE  HUMBER  OF  MINUTES 

"SEC"  IS  THE  NUMBER  OF  SECONDS.  SEC  IS  HOT 

TO  BE  ON  INTEBER. 


REVISION  HISTORY 


2  JUL  1979 


SENE  0.  MORE 


IHITIOL  VERSION. 


SUBROUTINE  DTIME < N . YEOR .DOY .HOUR ,MIN .SEC > 
DIMENSION  N<4> 

IHITIOLISE 

N1‘N<1) 

N2’N(2> 

N3‘N( 3) 

N4‘N<4) 

IF<HI .LT .0)N1’N1*4090 
IF  (N2.LT  .0)  N2^N2*409ii 
IF<N3.LT.0)N3^H3-^4O9S 
IF(H4.LT.0)N4mN4*4096, 

SECONDS  COLCULOTION 

Il'Nl/16 
12=11/ IS 
I3=N2/1S 
14=13/0 

SEC  =  1 0*  <  1 3~  0*  1 4  \  N2~  l<i*I3  >  4 .  l  2 C"" .  0 1  <  1 1  ~  1 S* 1 2  )  ^ 


FORT Ran  10  PROBRRM 


51 

$00 

52 

ISO 

5S 

54 

140 

55 

56 

57 

150 

58 

59 

160 

60 

oJ 

170 

62 

oS 

ISO 

64 

o5 

oo 

190 

67 

o$ 

o9 

200 

70 

71 

C 

72 

z 

7S 

c 

74 

700 

7  5 

76 

77 

710 

78 

C 

Tf> 

c 

SO 

c 

$1 

7  SO 

*?'•/ 
tr*  A 

$S 

$4 

SIO 

$S 

300 

S<i 

0? 


UR1T£( 4 . ISO) 

FORHRK’  ’> 

URl T£ ( 4 , 1 40 >  UPON < 1 > 

FORMRTC  UaO£UURS£R  ' 

UR2T£(4.J30) 

MRIT£(4,1S0) 

F-ORRRTC  T£RM  ORLU£  FOR  THIS  HOV£MUHB£R  * 
R£Ae(4,J60>F<l) 

F0RHf)T(B12.a) 

HR1T£<4,170) 

fORHRT<’  LIH£  STRENSTH  FOR  THIS  Hf)0£NUHB£R 
R£RD(4,JS0)C(I) 

FORHRTiSJS.S) 

URIT£(  4..JS0> 

UR1T£<  4,1 PO  ' 

format <'  IS  THIS  CORRICT  '  (O’H.l-Y)  ’■!> 

R£aD<4,100)10HS 

IFdRHS.tQ.O)  GO  TO  $00 

COHTINU£ 

GO  TO  7S0 

6£T  T£RH  ORLU£S  FROM  OLD  DATA  F1L£ 

R£AD(6,i00)lTQ 
DO  710  1*1, ITQ 
R£A0<6,S10)X,y,F(l >,€<:> 

COHTINU£ 

HRIT£  OUT  COMRL£T£  F1L£ 

URIT£<7 .lOOlITL 
DO  SOO  :=1,ITL 

HRIT£<7 ,S10)MA0N<I > ,RLRSP( I > .f(I),C(I> 
F0RMAT(Bl5.a,£15.a,Bj2.'i.D25.S> 

C0NT1NU£ 

£MD 


rORTRftN  H>  FROSRAH 


51 

IPRGSi  5  >-»<35 

c 

5S 

c 

IHPUT  DRIVE  UUMSER  RHD  INITIALIZE  TAPE  DRIVE 

5V 

c 

53 

MRI  TE<4. 500 ,) 

5d 

500 

FORNATCl  INPUT  TAPE  DRIVE  HUNSER  ".*) 

57 

REA0<4,510)IDRV 

S$ 

510 

E0RHAT<I4) 

59 

IE<IBRV.£Cl.l)ISRV^2 

60 

INITIO 

67 

IREC‘1 

o2 

IOPR‘4 

6J 

CALL  FFTFIH IN  I T . IDRV . IREC . lOPR , XRIN6< 1 ) .DUEC) 

c 

o5 

c 

DETERMINE  MHAT  RECORD  IS  TO  BE  READ 

66 

c 

o7 

NRIT£(4.520) 

6<? 

520 

EORMAT("  INPUT  PFT  E PANE  HUMBER  TO  BE  PROCESSED  ",*) 

69 

READ{4,5101IFILE 

70 

c 

77 

c 

READ  IN  HEADER  FOR  FRAME  HUMBER  “IFILE" 

72 

c 

7  J 

c 

FILL  IN  LATER 

7V 

c 

75 

c 

76 

c 

READ  IN  RECORDS  5  AND  <i  FROM  FRAME  "IF HE" 

77* 

c 

NRITE  THESE  THO  RECORDS  TO  DISK  AND  PHASE  CORRECT  THEM 

75 

c 

79 

MODE^l 

50 

IREC’5 

57 

CALL  FFTFIL  ( IF  I LE.  MODE  .  IP.EC ,  lOPR,  XRING(  I  > ,  BUFO  > 

52 

c 

5o 

c 

CONVERT  RECORD  5  TO  36  BIT  FORMAT 

5¥ 

c 

• 

5 

ilRITE<4.600i 

56 

oOO 

FORMAT ("  STARTING  TO  CONVERT  TO  36  BIT" > 

57 

CALL  SPCONVfSUFO  ,  BUFO  ,  20  4S .) 

55 

c 

59 

c 

TRANSFER  RECORD  TO  DISK  DRIVE 

90 

c 

97 

MRITE<4,alO) 

92 

olO 

FORMAK"  TRANSFERING  TO  DISK") 

9J 

IDATA*1 

9V 

ICYLO^O 

95 

lER^O 

96 

CALL  FFTRK<IDATA,IC7LD,1,BUF0,IER> 

97 

c 

95 

c 

REPEAT  FOR  RECORD  NUMBER  o 

99 

c 

00 

IREC=6 

fOPTKHN  pfitvafirtw 


101 

203 

203 

104 

105 

C 

106 

c 

107 

c 

106 

c 

109 

c 

110 

c 

111 

c 

113 
123 

114 

12^ 

630 

126 

227 

C 

116 

C 

219 

C 

120 

121 

630 

222 

C 

223 

c*»** 

124 

c 

125 

226 

227 

126. 

129 

130 

700 

131 

232 

233 
134 

235 

236 

710 

2  37 

C 

236 

C 

239 

C 

2  40 
141 

c 

143 

c 

143 

c**** 

144 

c 

145 

146 

147 

750 

146 

2  49 

C 

ISO 

c**** 

CaCL  FFTFIi  <  IFIL  Z .  HOD£  ,  IP.iC ,  lOPP ,  FP.IHG  .  GUFO) 

CALL  SPCONV  (  SUFO ,  SUFO ,  3046-  > 

ICVLTI^l 

CALL  FFTPK(  IDATA,ICYLD,I,BLlFO.I£P) 

BATA  IS  NOH  ON  DISK 

CALL  PHAS£  C0P.R£CT10H  PROGRAM  TO  £L1MINAT£  PHASE  SHIFT  IN  DATA 

4K  DATA  POINTS  TO  PHASE  CORRECT  THEREFORE  H*I2 
FIVE  PASSES  OVER  DATA  THEREFORE  NPASS  *5 


H^IS 

HPASS^S 

NRITE<4,630^ 

FORMATf"  S£ BINNING  PHASE  CORRECTION" ) 
CALL  PHAS£C  :N  ,NP.ASS,IPASS) 

PHASE  CORRECTION  COMPLETE 

NRITE<4,6S0) 

FORMAK"  PHASE  CORRECTION  COMPLETE") 

ASSEMBLING  REAL  PART  IN  ARRAY  HBUFF(I> 

CALL  FFTRK(IDATA,0,0,BUF0,IERy 
NJ’O 

DO  700  10*0 , I02S 
HSUfF(lJ)*SUFFER<NJ> 

NJ*NJ*-2 

CONTINUE 

CALL  FFTRKi IDATA.J.O^BUFO.IER) 

NJ*0 

DO  710  IJ*2024 .3047 
HSUFF( IJ )*BUFFER<NJ > 

NJ^NJ*3 

CONTINUE 

BRING  IN  HAMMING  TABLE  VALUES 
READ( 6 >BUFO .BUFFER 


READING  IN  EXPECTED  PEAK  SAMPLE  NUMBERS 

MRIT£<  4.I0O’ 

MRITE<4,750) 

FORMAK •  TOTAL  NUMBER  OF  LINES  YOU  ORE  LOOKING  FOR  -  '.f> 
R£AD<4.Z20HTL 

MRITE  TO  SEARCH  FILL  TOTAL  NUMBER  OF  LINES 


FORTRAN  ItJ  PROSRRM 


15: 

2  5-2 

C 

2  53 
254 

2  55 
256 

C 

257 

158 

514 

259 

755 

260 

262 

7  aO 

2o2 

163 

516 

2o4 

C 

2a5 

C***-» 

2a6 

167 

2aS 

169 

c 

270 

770 

172 

900 

272 

C 

275 

c**»* 

174 

175 

176 

c  ■ 

177 

27$ 

7S0 

279 

790 

2$0 

C 

1S2 

c 

1$2 

2S3 

C 

2$4 

235 

1$6 

793 

1$/ 

2SS 

c 

1S9 

796 

190 

2000 

191 

C 

192 

c**** 

293 

194 

c 

295 

296 

297 

c 

>IR1T£<  S.T‘^'0)ITL 

no  2000  IC=1.1TL 
HRITE<4,260> 

HRni<4,S14> 

fORHAT<'  INPUT  NRV£NUMB£R  IN  CM- 2  FOR  LINE  '.*) 

R£AI><4,526>NRUN 

NRITE<4,760>IC 

FORMAT!^  INPUT  SAMPLE  NUMBER  OF'. 14.’  LINE  =  ’,*> 

R£A0<4.S20>ISN 

F0RMAT<B25.6) 

NRITE<  4 , 260> 

OUTPUT  TO  SCREEN  POINTS  AROUND  EXPECTED  SAMPLE  HUMBER 

no  POO  IN*ISN-  20.ISN-t-20 
CALL  APODKIH.'JLj 
NP.ITE(4,770)IN.UL 
FORMATS  ’.la.'  ’ ,E2i.a> 

CONTINUE 

READ  ACTUAL  PEAK  SAMPLE  NUMBER  FROM  SCREEN  AND  INPUT  TO  DISH  FILE 

NRIT£<4.260) 

URIT£<4..7S0> 

FORMATC  INPUT  PEAK  SAMPLE  HUMBER  ’.*) 

R£AD<  4,790) IRSN 
FORMAT (16) 

IS  IT  CORRECT? 

NRITE( 4 ,793) 

FORMAT<’  HAS  THIS  CORRECT  (O^N.l^T)  ’ .*) 

READ(4,790)KC 
IF (KC.EO.O)  50  TO  755 

HRIT£<  5.796) IRSN.HAUN 
FOP.MAT(  16,B2S.6) 

CONTINUE 

REMIND  TAPE 

CALL  FFTFIL(0, IDRU . 2 . I OPR .XRING, BUFO) 


E 


239 


c  rOfiT'flflw  m  pRo&naH 


1  C  FORT RON  lU  PROGROM 

i  c 

3  C  PILE  NONE  =  HOMTBL 

4  C 

5  C  MRITTEN  BY  PORRIS  HEOL 

6  C  UTOH  STOTE  UHIVERSITV  HOY  22, 1904 

7  C 

S  C  PROBROH  TO  COMPUTE  THE  UOLUE  Of  THE  H0HH1H6  OPOBIZING 

9  C  FUNCTION,  UHICH  IS  NORNOLIZEB  FOR  UNIT  OREO. 

10  C  VOLUES  COHUPTEB  FOR  POSITIVE  OFFSET  ONLY. 

11  C  ORROY  OF  UOLUES  ORE  COMPUTED  FOR  0.01  SONPLE  NUMBER 

12  C  RESOULTION.  STORED  IN  ON  FILE  COLLED  HOMUL  OH  DISK 

13  C 

14  COMMON  BUFO, BUFFER! 2047} 

15  C 

16  DO  IS  1=0,2047 

17  BUFFER! I >=0.0 

IS  IS  CONTINUE 

19  DEL=0.0 

20  BUFF£R!0)=0.S3S36 

21  PI=S. 14159 

22  DO  10  1=1,500 

23  DEL*DEL*0.01 

24  fSTM=0.S3S36*! !SIN!PX»D£L) >/!PI»DEL>> 

25  SCTM=!SIN!P1*!D£L-1.0))>/!PI*!D£L-1.0>> 

26  TDTM*!SlN!PI*!D£L4l.O)>  >/!PI*!DEL*-l  .0>  ) 

27  SUFFER!! )*FSTM*!0. 46164/2. 0>*!SCTM*TDTM> 

2$  10  CONTINUE 

29  C 

30  C  NRITE  OUT  TOBLE 

31  C 

32  i4RIT£!S)BUF0, SUFFER 

33  END 

34 


Kl  Iv.  *.*  fc* 

0|  »0  K  o  </  "J  0  U  <  i»l 


ORTRON  IV  SUSROLniNE 


FORTRRN  :v  SUFROUTIN£ 

FILE  HRHE  -  M&J 

URJTTEH  BY  PARRIS  NEAL 

UTAH  STATE  UNIVERSITY  JULY  IS,J9$4 

PART  OF  “TPAUB"  THIS  IS  THE  3  POINT  RUNNING  AVERAGING  NINVON 
TO  AVERGE  TENPS  NEIGHTEB  BY  THEIR  STANBARD  DEVIATION 

SUBROUTINE  AVGT 

CONNON  IH1(5>  ,IH-2(5)  ,  IH3<  S>  ,  IL  1<  4  >  .  ILH  4)  .1L3<  4)  .A01<4)  ,A02i  4)  , 
A03( 4y .SAU 4> ,SA2< 4>,SA3<4i. TPi<4l . TP2i 4 > , TP3( 4 > , STK 4 > , STS ( 4 > . 
ST3( 4) ,TAVG<  4) .SAVG<  4^ 

BO  100  :^i,4 

TN-(TP:<I)/<STJ<I)**2)}*<TP2<I>/(ST2<I>**2>>4<TP3(I>/<STS<I)»*i)) 
TB-( 1 .0/ (STK I >**2> >*< 1 .0/(ST2< I >**2) >*< 1 .0/ <STS( I >**2) ) 
TAVG(I)=rTN/TB 
SAVG(I)^I.O/<SQRT<TD)) 

continue 

RETURN 


FORTRAN  SUBPOUT  I  HE 

FORTRAN  SUBROUTINE 

FILE  NAME  =  RERAN 

URITTEN  BY  PARRIS  HEAL 

UTAH  STATE  UHIOERSITY  JULY  JS,I9S4 

PART  OF  “TPAU6"  READS  FRAME  OF  TEMPERATURE  DATA 

SUBROUTINE  RF RAM( IH , IL ,A0 ,SA , TP, ST) 

DIMENSION  IH(  5)  ,IL^4i  ,/)0(4  >  .SAi  4)  .TP(  4)  .ST(  4> 

READiS.lOO^IHi : ' .IH<2> ,IH<Sy .IH(S> 

>  FORMAT i SI4) 

DO  200  1=1.4 

READ(S.1S0>IL(  I  >  .A0(  I  )  .SA(  IKTP<I>,ST<I) 

<  CONTINUE 

>  FORMAT^  I4,4EJS.u,> 


RE  TURN 


FORTRfiN  /V  SUBROUTINE 

FORTRAN  lU  SUBROUTINE 

FILE  NAME  -  UFRAM 

NRITTEN  BY  PARRIS  NEAL 

UTAH  STATE  UHIOERSITV  JULY  13,1964 

PART  OF  "TPAUB"  MRITES  A  FRAME  OF  AVERAGED  TEMP  DATA 

SUBROUTINE  MFRAM( IH , IL ,AO,SA, TP,ST> 

DIMENSION  IH<5>,IL<4> ,A0<4y .SA< 4 ) , TP( 4 ) , ST ( 4 > 

NRIT£(6. ICO  tIH( 1 > , IH<S) .IH(3) . IH<  4> . IH^ 5) 

FORMAT( 514  > 

DO  SOO  1=1.4 

MF.ITE  (6 .150>IL(  I )  ,A0<  I  >  ,SA<  I)  .TP(  I )  .ST<  I  > 

CONTINUE 

F0RMAT(I4.4E15.o) 


RE  TURN 


iff  C»J  (*4  KJr  ►  ^ «'  ^  ^4  **4  ^»4  ^  ^ 

►•>  ►^  o  >£.  c>  N  Ov  m  A  «►!>>*  ^  o  'i>  «&  '•j  0  ui  1*4  ►■>  K*  o  '0  N  o-  t*i  <  c^  »•:, 


FORTRfiN  lU  SUBROUTINE 


1  C  FORTRRN  10  SUBROUTINE 

C 

C  FILE  NOME  =  HOOE 

C 

C  MRITTEN  By  PORRIS  NEOL 

C  UTOH  STOTE  UNIOERSITV  JULY  1S,J9S4 

C 

C  PORT  OF  "TPOOB"  MOOES  TEMP  FROMS  3  OND  2  TO  2  OND  1 

C 

SUBROUTINE  MOOE 
C 

COMMON  IH1<S>  ,1H2(5)  ,1H3<S>  .1L2(4>  ,IL2(  ■»>  ,IL3(-i>  .001(4)  ,002(4) 

1  003(4) ,S01( 4) ,S02( 4) ,S03( 4) ,TP1( 4) ,TP2( 4) , 703(4) .ST  I ( 4 ) ,ST2( 4 ) 

2  ST3(4) .T00G(4),S00G(4) 

C 

DO  200  1=2.5 
IH2 ( I )=IH2( I ) 

1H2( I )=IN3( I ) 

200  CONTINUE 

C 

DO  200  1=2,4 
002(1)=002(1) 

002(I)=003( 1) 

S02(1)=S02(I) 

S02( I )=S03( I ) 

TP2(I)=TP2(I) 

TP2(1)=TP3(1) 

ST2(I)=ST2(I) 

ST2(I)=ST3(I) 

200  CONTINUE 
C 


RETURN 


Cj  Uj  U|  b|  U|  0|  G|  C»|  1-^ 


/FORTRON  lU  SINGLE  PRECISION  CONVERSION 

1  /FORTRRN  IV  SINGLE  PRECISION  CONVERSION 

2  / 

3  /FILE  NOME  IS  'SPCONV.RR’ 

4  / 

5  /GENE  R.  HRRE 

6  /ELECTRO-DYNRNICS  LRSORRTORIES 
T  /UTAH  STATE  UNIVERSITY 

S  /LOGAN,  UTAH  04322 
p  /3i  GECEHSER  1P77 
10  / 

11  /THIS  SUBROUTINE  CONVERTS  12  BIT  SINGLE  PRECISION 

12  /INTEGERS  (AS  LOADED  FROM  MAG  TAPE,  FOR  INSTANCE) 

13  /INTO  FORTRAN  IV  INTEGERS.  IT  ASSUMES  THAT  THE 

14  /ARRAY  OF  INTO  NHICH  THE  VALUES  HERE  LOADED  IS 

15  /DIMENSIONED  AT  LEAST  N/3  (ROUNDED  UP)  NHERE 
l<i  /N  IS  THE  HUMBER  OF  12  BIT  DATA  NORDS  TO  BE 

17  /CONVERTED.  THE  OUTPUT  ARRAY  MUST  BE  DIMENSIONED 
IS  /AT  least  N.  the  INPUT  AND  OUTPUT  ARRAYS  MAY  BE 

19  /THE  SAME  ARRAY. 

20  / 

21  J CALLING  SEQUENCE 

22  / 

23  /CALL  SPCONVdXIN.IXOUT.N) 

24  /  IXIN  »  INPUT  DATA  ARRAY 

25  /  IXOUT  *  OUTPUT  DATA  ARRAY 

26  /  N  •  HUMBER  OF  12  BIT  DATA  HORDS 

27  / 


2 '3 

SECT 

SPCONV 

‘39 

JA 

0ST 

30 

0XR. 

ORG 

.*10 

Ji' 

TEXT 

*SPCONV* 

33 

»RET. 

SETX 

0XR 

33 

SETB 

0BASE 

34 

JA 

.  *3 

35 

0BASE . 

ORG 

.  *6 

Jo 

IXIN. 

ORG 

.  *3 

jr 

IXOUT , 

ORG 

.  *3 

30 

N. 

ORG 

.*3 

39 

IXI. 

ORG 

.  *3 

40 

IXO. 

ORG 

.*3 

41 

ORG 

0BASE*3O 

42 

FNOP 

43 

JA 

0RET 

44 

FNOP 

45 

0GOBAK, 

OsO 

46 

0ARGS . 

ORG 

.*3 

47 

I. 

ORG 

.*0003 

4$ 

ITl , 

ORG 

.*0001 

4^f 

IT2. 

ORG 

.  *0001 

50 

ITS. 

ORG 

.*0001 

/FORTRaN  lU  S1H6LE  PRECISION  CONOER 


51 

J. 

ORB 

.  y-OOOS 

52 

»TMP, 

ORB 

.-f-OOll 

53 

*LIT, 

0000 

5j 

0000 

55 

0000 

5<S 

0001 

ST 

2000 

59 

0000 

59 

0002 

3000 

•il 

0000 

o2 

0003 

63 

3000 

o4 

0000 

o5 

0LBL=. 

o6 

■ 

ORB 

0LBL 

67 

»RTN, 

BRSE 

0BRSE 

69 

JR 

0BOBRK 

69 

*ST, 

STRRTV 

70 

■  0210 

71 

FSTR 

0BOBRK , 0 

72 

0200 

73 

SC  TX 

0XR 

74 

3ETB 

0BRSE 

75 

LDX 

0,1 

76 

FSTR 

0BRSE 

/  i 

FSTR 

0RRBS 

79 

FLORZ 

0BRS£,1* 

79 

FSTR 

IXIN 

JO 

FLORZ 

0BRSE,1* 

91 

FSTR 

IX  OUT 

92 

FLORZ 

0BRSE.1* 

93 

FSTR 

N 

94 

/ 

95 

/INITIRLIZE 

96 

/ 

97 

STRRTF 

99 

FLORZ 

N 

99 

FSTR 

I 

90 

FLOR 

0L I T*0006 

91 

FHULZ 

N 

92 

FSTR 

J 

93 

/ 

94 

/HRIN 

LOOP 

95 

/ 

96 

010, 

FLOR 

I 

97 

FSUB 

0L  IT *-0006 

■  FORTRAN  HJ  SINGLE  PRECISION  CONOEP.SIOH 


101 

FOOD 

IX  IN 

lOS 

FSTF) 

I XI 

103 

STRRTF 

104 

FLDAZ 

1X1 

105 

FSTR 

ITl 

106 

FLDR 

J 

107 

FSUS 

»L  It  4-0006 

lOS 

FSTR 

J 

109 

JLT 

020 

110 

/ 

111 

/FIRST 

UORD 

112 

/ 

113 

RLN 

0 

214 

STAR  TO 

US 

FADO 

IXOUT 

116 

FSTA 

1X0 

117 

ST  ART F 

US 

SETX 

ITS 

119 

XT  A 

0 

120 

SETX 

0XR 

121 

FSTAl 

1X0 

122 

FLOA 

J 

123 

FSUS 

0LIT4OOO6 

224 

FSTA 

J 

12S 

JLT 

020 

126 

/ 

127 

/SECOND 

•  HORO 

12S 

/ 

129 

ALH 

0 

1 30 

starts 

131 

FAOD 

IXOUT 

132 

FSTA 

1X0 

2  33 

ST ART F 

134 

SE  TX 

IT2 

135 

XT  A 

0 

4  J  o 

SETX 

0XR 

137 

FSTAZ 

IXO 

13S 

FLOA 

J 

239 

FSUS 

0L  1  T4-0006 

140 

FSTA 

J 

141 

JLT 

020 

2  42 

/ 

143 

/THIRD 

HORO 

144 

/ 

145 

ALH 

0 

146 

STARTS 

147 

FADS 

IXOUT 

14S 

FSTA 

1X0 

140 

STARTF 

ISO 

SETX 

IT  I 

•_  *. 


Cl  <  U  -<!»  4^  4^  -<  ■<»  (»4  K4  A«  k.  k<  is*  K«  K« 

O  O  •»  "I  O'.  CH  ^  C*|  Ki  O  C)  'i  0  tl  ■•*  *►!  ».>  *^4  O  *  ■»  N  0^  «  A  *»4  »■*  *^  O  &  NO  U*  <  K' 


!<£B£C  9000  fOfiTPan  10  TPP£  C-UBfiOUTlN£ 


A'fSfC  9000  FOP.TP.PH  10  TPP£  SUBPOUT  1N£ 

F1L£  HaM£  IS  -TOPE' 

6£N£  P.  HOPE 

ELECTPICPL  EN61NEEP1H&  BEPPPTPENT 
UTOH  STATE  UHIUEPSITV 
LOGAN,  UTAH  $4222 
17  JULT  1979 

CALLING  SEQUENCE 

CAL  L  TAPE  (■  UNI  T  ,  COM  AND  ,  BUFFEP  ,  COUN  T  > 

UNIT  =  TAPE  UNIT  TO  BE  USED 

0  «  PHYSICAL  UNIT  0  IN  9T  MOOE 

1  -  PHYSICAL  UNIT  O  IN  7T  MODE 

2  »  PHYSICAL  UNIT  1  IN  91  MODE 

3  =  PHYSICAL  UNIT  1  IN  7T  MODE 

CONAND  -  TAPE  COMMAND 

0  *  NO  OPEPATION 

1  -  NO  OPEPATION 

2  '  REMIND 

3  =  SET  OFFLINE 

4  *  SPACE  FOPMAPD  N  FILES 

5  *  SPACE  FOPMAPD  N  RECORDS 

6  »  READ  FOPMAPD  1  RECORD 

7  ^  R£AD  FOPMAPD  1  RECORD 
$  »  SPACE  REVERSE  N  FILES 

9  *  SPACE  REVERSE  H  RECORDS 

10  »  SPACE  REVERSE  N  RECORDS  <£DIT> 

11  *  READ  REVERSE 

12  *  MRITE  file  HARK 

13  =  ERASE  3  INCH  GAP 

14  =  MRITE  RECORD  (EDIT' 

15  -  MRITE  RECORD 
BUFFER  =  BUFFER  ARRAY  FOP  DATA 

COUNT  »  MORD/RECORD  COUNT.  N0TE--F0R  READ 
AND  MRITE  COMMANDS.  THIS  IS  THE  12 
BIT  MORD  COUNT. 

REVISION  HISTORY 


17  JUL  1979 
14  AUG  19$2 


GENE  A.  MARE 
GENE  A.  MARE 


INITIAL  VERSION 

CHANGED  CALL  TO  TAPES(-1,... 

TO  HOLD  IN  MODULE  UNTIL  DONE 


SUBROUTINE  TAPE ( UNI T .COMAND .BUFF ER .COUNT > 
GET  THE  ARGUMENTS  IN  LOCAL  VARIABLES 


XEBEC  9000  FORTROH  IV  TFIPE  SUBROUTINE 


I 

tv 

kV 


I.-' 

I 

K.- 

►  -■ 

t':*; 

ii 

t:-- 


SI 

c 

5^' 

1UN17=UN1T 

53 

ICOH*COMRHV 

54 

JC0M‘IC0»X2 

55 

K0UHT^C0UHT*0.5 

56 

1-1 

57 

c 

59 

c*** 

BO  IT 

59 

c 

60 

10 

CONTINUE 

61 

CALL  TRP£3< 1 ,IUN1T  ICON .SUFFER .KOUNT > 

62 

c 

63 

c*** 

CHECK  FOR  RERO/NRITE  OPERATIONS 

64 

c 

6  5 

IF(<JCON.EO.3).0P..{JC0N.EO.7y'>GO  TO  20 

66 

SO  TO  30 

67 

c 

60 

c**» 

CHECK  FOR  ERRORS  kP.EAB/NRITE  ONLY) 

69 

c 

70 

20 

CONTINUE 

71 

CALL  TAPE 3(-l.KRES, KCOH , KBUF . KCN T > 

79 

INTN-<KRES/256>-<KRES/S12>*2 

73 

1F<  INTN.NE.OGO  TO  20 

74 

KRES-KREe/4096 

7  5 

KRE6-<KRE6/64>*iKRES/3-<KRES/16>»2) 

76 

IF<KRE6.E0.0.0R.I.0T.9)S0  TO  30 

77 

c 

79 

c»** 

BACKSPACE  AHB  TRY  AGAIN  UP  TO  10  TINES 

79 

c 

90 

CALL  TAPE3( 1 .lUNlT .9,BUFFER.1> 

$1 

1-1*1 

32 

GO  TO  10 

33 

c 

34 

c*** 

BONE 

e 

c 

•;« 

JO 

CONTINUE 

•if  1 

RETURN 

0  $ 

ENB 

09 


c 

FORTRAN  10  PROGRAM 

: 

C 

FORTRAN  20  PROGRAM 

2 

C 

3 

C 

FILE  name  =  TPAOG 

4 

C 

5 

C 

NRITTEN  BY  PARRIS  NEAL 

<i 

C 

UTAH  STATE  UNIOERSITY  JULY  13,1934 

7 

C 

■3 

c 

PROGRAM  READS  TEMPERATURE  FILES  PRODUCED  BY  "BLTZMD’ 

c 

AND  AOERAGES  THEM  OOER  A  3  POINT  UINDOM  UEIGHTED  BY 

10 

c 

EACH  TEMP'S  STD  DEO 

11 

c 

1 

COMMON  INKS)  ,IH2<5)  ,  IH3  (  5>  .  ILl  <  4  )  ,  IL2<  4  >  ,  IL3  (  4  >  .A01(4>  .A02<4> 

1 3 

1 

A03<4) ,SA1<4> ,SA2(4y ,3A3(4> . TP1<4>, TP2(4> , TP3(4) ,STU 4) ,ST2<4> 

14 

^  r 

1 

ST3<4)  ,  7A‘0G(4)  .SA0G<  4  > 

16 

P.EAD(  5. 100>  INFMS 

17 

100 

FORMAT! 14^ 

IS 

MR1TE(6, 100) INFMS 

19 

c 

20 

CALL  RFRAHdHl  ,IL1  ,A01  ,SA2  ,TP1  ,ST1> 

■21 

CALL  RFRAM(1H2,IL2,A02.SA2, TP2,ST2) 

c 

23 

c 

24 

DO  200  1*1,4 

25 

TH*<  TPKDKSTl  ( 1)**2)  >*<  TP2(I )/<ST2<I  >**2)  ) 

26 

TD*<  1 .0/<STl(l)**2))4-(1.0/(S72(l>»»2>> 

27 

TA0G(1)*TN/TD 

2S 

SAOG< I )*1.0/(S0RT(TD) > 

29 

200 

CONTINUE 

30 

c 

31 

CALL  MFRAM< IHl ,IL1, AOl ,SA1 , TAOG , SAOG) 

32 

c 

33 

500 

CONTINUE 

1 

34 

CALL  CHKEOF(E) 

35 

READ! S,120)1H3< 1 ) . 1H3<  2) ,IN3( 3  > .IN3i 4> , IH3( 5> 

3a 

IF!E.NE.O)  GOTO  600 

37 

120 

FORMAT (514) 

33 

DO  130  10=1,4 

39 

RCAD(5,140)IL3!IJ),A03i IJ) ,SA3! IS) , TP3! IJ ) ,ST3 ! IJ ) 

40 

130 

CONTINUE 

41 

140 

FORMAT! 14,4E15.6) 

42 

C 

43 

CALL  AOGT 

44 

C 

45 

CAL L  MFRAM! IH2,1L2. A02 , SA2 , TAOG , SAOG ) 

46 

CALL  MOOE 

47 

c 

4S 

GOTO  500 

49 

aOO 

CONTINUE 

Ml  Ml 


/FOfiTRaN  ;v  FFT  KUOS  HOOULE 


t 

/FOK.KKN  10  FFT 

RFC 5  NODULE 

s 

/6Z.Z  f) 

.  N'2Z 

4 

/£ . ZM42 

XU/.  £.• 

>.IZ2.).‘  DZ0'24-Z.4 

S 

/U4.'t  S 

4/4Z  U.)6Z23>46 

o 

/L/' ‘ 

U4f (  S4322 

/ 

/30  JS. 

t 

Z  19$  2 

4 

/CKLLIHB  SEQUENCE 

JO 

/ 

11 

/ 

CALL  FFTRK(UN1T ,CYLnP..P.DNT ,X,£RAOP> 

12 

/ 

13 

/ 

UNIT  *  PHYSICAL  DISK  UNIT  NUMBER  (i 

14 

/ 

CYLDR  =  CYLINDER  TO  BE  TRANSFERED 

IS 

/ 

RDNT  =  READ/NRITE  SELECT  <0  -  READ 

16 

/ 

X  ^  BUFFER  USED  <2049  F4  NORDS> 

IT 

/ 

ERROR  =  RKOS  STATUS  REGISTER 

IS 

/ 

19 

/ 

20 

/«£C/7S; 

ON  NJ STORY 

21 

/ 

•><> 

/JO  J5. 

Z  1992 

GZ.Z  A.  N.'2Z  I.)4>/,  6Z23)/. 

25 

/ 

24 

/ 

25 

/DEFINITIONS 

26 

/ 

27 

DSNP^ 

6741 

/DISK  SKIP  ON  FLAG 

2S 

OCLP.* 

6742 

/DISK  CLEAR 

29 

DLBB* 

67 43 

/LOAD  ADDRESS  AND  GO 

30 

DLCfl‘ 

6744 

/LOAD  CURRENT  ADDRESS 

31 

DFIST-r 

6745 

/READ  STATUS 

32 

DLDC^ 

67 46 

/LOAD  COMMAND 

33 

BSN^ 

7002 

/BYTE  SNAP 

34 

/ 

3  5 

/FFTP.K 

ENTRY 

36 

/ 

37 

BASE  0 

3S 

SECT  FFTRK 

39 

SETX  CON  /SET  BASE  ADDRESS 

40 

STARTD 

/GET  UNIT  NUMBER 

41 

FLDAZ  0 

.1 

42 

FSTA  3 

43 

ST ART F 

44 

FLDAZ  3 

45 

ALN  0 

46 

STARTD 

47 

FSTA  UNIT 

4S 

FLDAZ  0 

,2  /GET  CYLINDER  NUMBER 

49 

FSTA  3 

50 

STARTF 

;'02.> 

-■  UKITE) 


/foktrun  iu  fft  rkos  module 


51 

FLDRZ  3 

52 

RLH  0 

53 

STRRTD 

5d 

FSTR  C'/LDR 

55 

FLDRZ  0.3 

,GET  RERD/MRITE 

56 

FSTR  3 

57 

STARTF 

58 

FLDRZ  3 

59 

RLN  0 

60 

STRRTD 

61 

FSTR  RKMT 

62 

FLDRZ  0.4 

,3ET  BUFFER  ADDRESS 

63 

FSTR  3 

64 

FSTR  BUFH 

STRRTF 

IZifi 

TRRP4  RK05S 

/BO  TO  S-MODE  SUBROUTINE 

67 

STRRTD 

/RETURN  STATUS  REGISTER 

63 

FLDRZ  0.5 

6R 

FSTR  3 

70 

STRRTF 

71 

fldr  error 

f 

FNORM 

7Z 

FSTRZ  3 

74 

FLDR  30 

75 

3RC 

/RETURN 

7  6 

/ 

7  / 

IS- MODE 

RK05  MODULE 

78 

/ 

79 

COHMZ  RK05M 

SO 

ORB  .*200 

31 

/ 

32 

/RKOS  CONTROL  SUBROUT 

INE 

S3 

/ 

34 

RH05S. 

0 

3  5 

/ 

/INITIOL  SETUP 

/ 

33 

CLL  CLR  IRC 

/DO  A  POMER  CLEAR 

39 

DCLR 

OO 

CLR  IRC  RTR 

/MAKE  LO  PART  OF  READ  RLL/MRITE  ALL 

91 

TRD  RDMT 

/GET  HRITE  BIT 

92 

RHD  K4001 

/MAKE  SURE  ITS  OK 

93 

CLL  RTR 

/AND  MOVE  IT  INTO  POSITION 

94 

DCR  COM 

/IN  THE  COMMAND  MORD 

95 

TRD  C'/LN 

/NON  SET  THE  CYLINDER  NUMBER 

96 

RHD  N377 

/RESTRICT  IT  TO  3  BITS 

97 

CLL  RRL 

/AND  PUT  IT  IN  the  PROPER  POSITION 

93 

RTL 

6  O 

RTL 

100 

DCR  C'/LN 

For  the  disk  address  register  acv> 

/FORT Ran  10  FFT  RKOS  MOHLUE 


201 

TAD  UNUM 

102 

AND  K3 

20S 

SNA 

104 

JMP  RKRET 

105 

RAL 

106 

TAD  COM 

107 

DCA  COM 

lOS 

TAD  M30 

109 

DCA  SLKCNT 

110 

DCA  ERP.EG 

111 

lOF 

112 

/ 

113 

/MAIN 

DATA  TP.AHSFEl 

114 

/ 

215 

RLOOP. 

DCLP. 

1  !•/, 

TAD  BUFL 

117 

DLCA 

110 

TAD  BUFH 

119 

AND  K7 

120 

CLL  RAL 

121 

RTL 

122 

TAD  COM 

123 

DLDC 

124 

TAD  CYLN 

125 

DLA6 

126 

ISZ  CYLN 

127 

CLL  CLA 

12S 

TAD  BUFL 

129 

TAB  K400 

130 

DCA  BUFL 

131 

RAL 

132 

TAD  BUFH 

133 

DCA  SUFH 

134 

RHAIT, 

DSKP 

1 3  5 

JMP  CNTLC 

1 36 

DRST 

137 

AND  K3777 

1 30 

32  A 

139 

JMP  RERR 

140 

ISZ  BLKCNT 

142 

JMP  RLOOP 

142 

RMAITl 

,  DRST 

143 

SMA  CLA 

144 

JMP  RMAITl 

245 

ISZ  BELAY 

146 

JMP  ,-l 

147 

RKRET. 

CDF  CIF  0 

140 

ION 

149 

JMPl  P.K05S 

150 

DELAY , 

0 

UOU  SET  THE  UNIT  HUMBER 
■UNITS  0-S  ONLY 
/  :;>  NOT  ALLON  UNIT  0 
TO  PROTECT  OPERATING  SYSTEM 
.PUT  IT  TOGETHER  NITH  HI  CYLINBER  SIT 
/THEN  INCLUDE  IT 
/IN  THE  DISK  COMMAND 
/SET  THE  BLOCK  COUNT 
/TO  ■24C103  SLOCKS  (6144  NORDS) 

/.RESET  THE  ERROR  MORD 
. SORRY,  NO  INTERRUPTS  ALLONED 


/CLEAR  STATUS 

/SET  LO  ORDER  SUFFER  ADDRESS 
/AND  LOAD  INTO  CURRENT  ADDRESS 
/SET  HI  ORDER  SUFFER  ADDRESS 
/LIMIT  TO  FIELDS  0-7 

/then  shift  it 

/AND  PUT  IN  COMMAND  MORD 
/LOAD  THE  COMMAND  REGISTER 
/SET  THE  CURRENT  CYLINDER  NUMBER 
/LOAD  IT  AND  GO 

/UPDATE  CYLINDER  MHILE  MAI  TING 
/MUST  CLEAR  LINK  AN'/UAV 
/SET  lO  order  suffer  ADDRESS 
/AND  SUMP  IT  BY  TMO  MEMORY  PAGES 

/NON  GET  THE  OOERFLOM 
/AND  ADD  IT 

/TO  THE  HI  ORDER  SUFFER  ADDRESS 
/IS  THE  DISK  DONE' 

/CHECK  OUT  CONTROL  C  MHIlE  MATTING 
/CONE  —  GET  THE  STATUS  REGISTER 
/REMOVE  COMPLETION  FLAG 
ANY  ERRORS' 

/■■ES — GO  CHECK  THEM  OUT 

/HAVE  ME  TRAHSFERED  ALL  THE  BLOCKS'’ 

/NO — DO  ANOTHER 

/MAIT  FOR  THE  COMPLETION  FLAG 


/TURN  THE  INTERRUPTS  BACK  ON 
/A^L  DONE 


In  >  'O  0  K  0^  O  ’'I  ri  M  >  'O  0  lx  *>  O  *'♦  M  fn  >  %■>  '/)  K  <•>  •>■  O  <4  >o  >  0  »>  <•'/  '5  O  •^i  M 


/FOF.TP.dN  ;y  i'/^r  fiA'05  MODULI 


1  52 

tRP^OR  Cn£Cr\ 

P.LP.P..  DCP  tPREG 
TAD  CP.P.E6 
RND  K^OS 
SNR  CLP 
JMP  PKP.£T 
DCLP 
STL  RTL 
DCLP. 

DSKP 
JMP  .-1 
DCLP. 

DP.ST 
SIR  CLR 
JMP  .-5 
JMP  P.KP.ET 

/ 

, CONTROL  C  CHECH 
/ 

CMTLC,  KP.S 

PHD  H277 
TRD  MS 
SNR  CLR 

KSF 

JMP  P.HR1T 
JMP  P.KP.ET 

/PRP.RMETEP.S 

/ 

COM,  0 

ONE .  2 

7  HO,  I 

7HRE .  S 

FOUR.  N 

FIVE,  t 

UNIT,  0 

UNUM .  O 

C'.'LDP..  0 

C'/LN,  0 

RUNT,  0 

P.DNT,  0 

29 J  DUFH.  0 
205  BUFL,  0 
29-:,  ERROR.  i7 

2  97  0 

2  90  EP.P.EG.  0 
2-'9  BLHCNT  .  0 


SPUE  THE  STATUS 
■DO  NE  NEED 
/TO  RECRLIBRRTE? 

/NO — STOP  Transfer  and  return 

/START  THE  RECRLIDP.RTE--CLEAR  DRIVE 

/AC  *•  2 

/RECALIBRATE 

/NRIT 

/UNTIL  DONE 

/CHEAP  STATUS  REGISTER 
/NRIT  FOP.  STATUS 

/still  doing  RECALIBRATE 
/NON  STOP  TRANSFER  AND  RETURN 


/READ  KEYBOARD  STATIC 
/REMOVE  PAP.ITY  BIT 
/IS  IT  A  '■C? 

/IT  IS  A  —  IS  THE  FLAG  UP'' 

/NO  —  KEEP  ON  TRUCKING 

/A  'C  —  LET  FRTS  TAKE  CARE  OF  IT 


A 

>. 

UJ 

r 

,■% 

s" 

s* 

.V 

niES 

£C  ‘-'OOO 

FOPTPPN  m  S-MOBE  E-UBPOUr I tIES 

,v 

m  * 

/:{EBEC 

9000  FOPTPPN  m  i-NOBE  SUBPOUT^HES 

3 

/ 

/FILE 

NPHE  IS  ’TPPES.PP’ 

% 

4 

/ 

y 

5 

/  QEHE 

P .  HAA£ 

6 

/ELECIRICF)L  EneiHEEP.INO  BEPPPTMENT 

/' 

/UTPH 

Srflr£  UNtUEPSIT'/ 

$ 

/LOGPN 

,  y7fl//  S432S 

9 

/1 2  JULV  ly'TF) 

:o 

/ 

r 

11 

/CPLilHS  SEQUENCE 

12 

/ 

13 

! 

CALL  TPPE&i  1 , 7i/N7,  ;£/:/«,  7Bi/f  .  7CN7  ) 

14 

/ 

7  =  7/)P£  OPEPPTION  MODE 

1  e 

/ 

I  LINT  -  TAPE  UNIT  NUMBEP 

y. 

1^ 

/ 

ICON  -  TAPE  CONNPNB 

IT 

/ 

ISLIf  =  TAPE  BUFFEP  NAME 

1& 

19 

/ 

/ 

ICNT  =  NOP.B/PECOPD  COUNT 

m 

20 

/ 

CAL L  TAPES ( 0 , lUNT . IF  NT , PPBN .PTFS} 

■/ 

21 

/ 

0  -  TAPE  PAPAHETEP  SET  NOBE 

•>2 

/ 

lUNT  =  TAPE  UNIT  HiJNBEP 

23 

/ 

IFHT  »  TAPE  UNIT  FOPNAT 

24 

/ 

PPBN  -  TAPE  UNIT  PAPIT'/~BENSITy 

£ 

/ 

PTFS  -  TAPE  UNIT  P£AD  THPESHOLB- F AST  SAP 

26 

/ 

27 

/ 

CALL  TAPES/- 1  ,IPeS,TCON.JBUF  .SCNT.) 

t. 

/ 

-7  '  TAPE  PESISTEP  PEAB  MOVE 

2  9 

/ 

IPES  '  TAPE  UNIT  EPPOP  PEGIST EP.-OPEPAT INS  STATUS 

30 

/ 

TCOH  •=■  TAPE  COHHANB  UOPB 

31 

/ 

JSUF  »  TAPE  BUFFEP  ABBP.ESS 

.V 

3^ 

/ 

JCHT  =  PESIBUAL  UOPB/P.ECOPB  COUNT 

V 
•  ■ 

33 

/' 

,  4 

3  4 

/lUNT 

TAPE  UNIT  SALLIES 

3  5 

/ 

0  -  PHYSICAL  UNIT  0  IN  ’■IT  MOVE 

% 

J-j- 

/ 

1  -  PHYSICAL  UNIT  0  TH  77  HOBE 

37 

/ 

S  =  PHYSICAL  UNIT  1  IN  ■■■■T  MODE 

30 

/ 

j  -  PHYSICAL  UNIT  J  IN  77  PODE 

39 

/ 

% 

40 

/icon 

TAPE  CONN AND  VALUES 

41 

/ 

0  »  NO  OPEPATION 

42 

/ 

I  ^  NO  OPEPATION 

43 

/ 

2  =  PEMINB 

r 

44 

/ 

3  *  SET  OFFLINE 

{^■a 

45 

/ 

4  -  SPACE  FOPNAPB  N  FILES 

4ty 

/ 

5  »  SPACE  FOPNAPB  H  PECOPBS 

47 

/ 

6  =  PEAB  FOPNAPB  1  PECOPB 

40 

/ 

7  ‘  PEAB  FOPNAPB  1  PECOPB 

4^' 

/ 

S  '  SPACE  F.E'-JERSE  N  FILES 

50 

9  =  SPACE  P.EUEPSE  N  P.ECOPBS 

js 

ft 

260 


/>i£B£C  -^000  FURTRON  I'.l  H-HODE  SUBROUTINES 


51  / 


10  =  SPACE  REVERSE  N  RECORDS  (EDIT) 


5S  / 

5S  / 

54  / 

55  / 

56  / 

57  / 

5S  /IFMT 

59  / 

60  / 
ol  / 

62  / 

63  / 


11  =  READ  REVERSE 

12  =  NRITE  FILE  NARK 

13  =  ERASE  3  INCH  GAP 

14  -  NRITE  RECORD  ^£DIT> 

15  •  URITE  RECORD 

TAPE  FORMAT  MODE  CONTROL 
0  =  UNPACKED  DATA  FORMAT 

1  •»  TEST  DATA  FORMAT 

2  =■  PACKED  FORMAT  A 

3  -  PACKED  FORMAT  S 


64  /PRDN 

65  / 

66  / 

67  / 

6S  / 

69  / 

70  /RTFS 

71  / 

72  / 

73  / 

74  / 

75  / 

76  /IRES 

77  / 

7S  / 

79  / 

■30  /  TCOM 

31  / 

32  HOT  INSTRUCTION  DEFINITIONS 


33 

34 

/ 

SKNB^ 

63  54 

/SKIP 

ON  CONTROLLER  NOT  BUSY 

0  5 

LDMAa 

6335 

/LOAD 

MEMORY  ADDRESS 

36 

LDMF‘ 

6333 

/LOAD 

MEMORY  FIELD 

37 

LHNC^ 

6311 

/LOAD 

HI  MORD  COUNT 

■33 

LDNC  = 

6327 

/LOAD 

LO  NORD  COUNT 

■39 

LDCM^ 

6325 

/LOAD 

COMMAND  REGISTER  AND  GO 

90 

RDNC^ 

6347 

/READ 

RESIDUAL  NORD  COUNT 

91 

RHNC‘’ 

6312 

/READ 

HI  RESIDUAL  NORD  COUNT 

92 

RDST* 

6343 

/READ 

STATUS  REGISTER 

93 

RDES^ 

6345 

/READ 

ERROR  REGISTER 

94 

DSIN* 

6351 

/DISABE  TAPE  INTERRUPT 

95 

BSN^ 

7002 

/BYTE 

■SNAP  ■ 

96  / 

97  / 

93  /REVISION  HISTORY 
69  / 

100  /12  JUL  1976  5£«£  INITIAL  VERSION 


TAPE  PARITY-DENSIT/  CONTROL 
0  =  EVEN  PARITY-LON  DENSITY ( NRZI > 

1  =  EVEN  PARITY-HIGH  DENSITY  (PE) 

2  -  ODD  PARITY-LON  DENSITY  (NRZI) 

3  =  ODD  PARITY-HIGH  DENSITY  (PE) 

TAPE  READ  THRESHOLD-FAST  SAP  CONTROL 
0  =  NORMAL  THRESHOLD-NO  FAST  SAP 

1  •=•  NORMAL  THRESHOLD-FAST  SAP  SET 

2  »  LON  THRESHOLD-NO  FAST  SAP 

3  LON  THRESHOLD- FAST  GAP  SET 

TAPE  REBISTERS— LON  ORDER  12  BITS  CONTAIN  THE  OPERATING 
REGISTER.  THE  HIGH  ORDER  12  BITS  CONTAIN  THE  ERROR 
REGISTER. 

TAPE  COMMAND  REGISTER  CONTAINED  IN  LON  ORDER  12  BITS. 
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aEB 

£C  9000  FORTRflH  10 

S-MOBE  SUBROUTINES 

01 

/ 

OH 

/ 

03 

/TRPES  ENTRY  (fPP 

INTERFRCE) 

OJ 

/ 

05 

BRSE  0 

06 

■  SECT  TRPES 

07 

SETX  HOBE 

OS 

STRRTD 

/GET  MODE  ORLUE 

09 

PLDRZ  0.1 

10 

FSTR  3 

11 

STRRTF 

12 

FLBRZ  3 

13 

RTX  0 

/STORE  HOBE  ORLUE 

14 

JLT  Ml 

/GO  TO  HOBE  =  -1  SECTION 

1 5 

STRRTB 

/GET  REHRINBER  OF  RRGUMENT. 

16 

FLBRZ  0,2 

/GET  NORB  2 

17 

FSTR  3 

IS 

STRRTF 

19 

FLBRZ  3 

20 

RLH  0 

21 

STRRTB 

22 

FSTR  MB2H 

23 

FLBRZ  0,3 

/GET  NORB  3 

24 

FSTR  3 

25 

STRRTF 

26 

FLBRZ  3 

27 

RLH  0 

2S 

STRRTB 

29 

FSTR  NOSH 

30 

FLBRZ  0.5 

/GET  NORB  5 

31 

FSTR  3 

32 

STRRTF 

33 

FLBRZ  3 

34 

RLH  0 

3  5 

STRRTB 

36 

FSTR  NOSH 

37 

FLBRZ  0.4 

/GET  NORB  4  RBBRESS 

3S 

FSTR  3 

39 

STRRTF 

40 

XTR  0 

/GET  HOBE  RGRIH 

41 

JEO  HZ 

/GO  TO  HOBE  »  0  SECTION 

42 

/ 

43 

YTRPE  CONTROL  SECT 

ION  <H0BE  -  1) 

44 

/ 

45 

STRRTB 

46 

FLBR  3 

47 

FSTR  JSUFH 

4S 

STRRTF 

4^1 

TRRP4  PBPT 

50 

FLBR  30 

/XBlSeC  i‘000  FORTRdN  lU  ■'■-MOVE  SLUSPOU TINES 


2z: 

JRC 

iZS 

/ 

1S3 

/TAPE  PPRPNETER  INIT 

lALIZATlON  SECTION  '.MODE  =  0> 

IS-t 

/ 

1S5 

»Z.  ELDRZ  3 

/GET  NORD  4  VALUE  <JBUF> 

156 

ALH  0 

157 

STRRTD 

«  «*•;< 

FSTR  MD4H 

159. 

STRRTF 

160 

TRRP4  PDPZ 

/GO  TO  ZERO  S-MODE  SECTION 

161 

FLDR  30 

lai 

JRC 

loS 

/ 

164 

/TAPE  RESISTER  READ 

SECTION  (MODE  ■»  -I) 

165 

/ 

2o6 

Ml,  TRAP 4  PnPMl 

/SO  TO  S-MODE  SECTION  TO  RE 

167 

STARTD 

16$ 

FLDAZ  0,2 

/RETURN  NORD  2  (IRES) 

169 

FSTA  3 

170 

ST ART F 

17 1 

FLO A  M02 

172 

FNORM 

17S 

FSTAZ  3 

174 

STARTD 

•?75 

FLDAZ  0,3 

/RETURN  NORD  3  (TCOM) 

176 

FSTA  3 

177 

ST ART F 

17S 

FLDA  1403 

179 

FHORM 

ISO 

FSTAZ  3 

ISl 

STARTD 

1 S2 

FLDAZ  0,4 

/RETURN  NORD  4  (JDUF) 

1  S3 

FSTA  3 

1$4 

STAPTF 

1S5 

FLDA  JDUF 

1 S6 

FNORM 

1 S7 

FSTAZ  3 

ISS 

STARTD 

1S9 

FLDAZ  0,5 

/RETURN  NORD  5  (JCNT) 

190 

FSTA  3 

191 

STARTF 

192 

FLDA  UD5 

193 

FNORM 

194 

FSTAZ  3 

195 

FLDA  30 

196 

ZAC 

197 

/ 

19S 

/e-MODE  TAPE  control 

SECTION 

199 

/ 

1:00 


COMHZ  RliOSM 


'XEBEC  yOOO 

FORTRPH  10  S- 

NODE  SUBROUTINES 

?Oi 

ORG  .*400 

?02 

/ 

?OJ 

/$-fiOV£ 

TAPE  CONTROL 

(NODE  =  1> 

?0-/ 

/ 

?05 

PBPT, 

0 

>06 

CLA 

i 

>or 

TAD  ND2L 

/GET  THIS  UNIT  NOP  CONNANB  POINTER 

i 

;os 

AND  K3 

1  i 

109 

TAD  TN0P*1 

>10 

DCA  TNOP 

i 

111 

TAD  MD2L 

/GET  THIS  UNIT  LAST  CONNAND  POINTER 

,  i 

112 

AND  K3 

113 

TAB  TC0M*1 

i 

111 

DCA  TCOM 

i 

115 

lOF 

/NO  INTERRUPTS 

•  - 

>16 

SKNB 

/ISSUE  HOP  CONNAND  UHEN  CONTROLLER  READ 

1  - 

117 

JMP  .-1 

1  : 

110 

TADZ  TNOP 

‘  i 

119 

LDCM 

t  - 

120 

CLA 

>21 

RDST 

/NAIT  UNTIL  THIS  UNIT  IS  ON  LINE.  READY 

122 

AND  K7000 

/AND  NOT  RENINDIHG 

123 

TAD  N6000 

1  ^ 

121 

SZA  CLA 

1 

125 

JNP  .-4 

1  ^ 

126 

TAB  JSUfL 

/LOAD  LO  BUFFER  ADDRESS 

"'iT 

LDNA 

‘  i 

120 

CLA 

129 

TAD  JSUFH 

/LOAD  HI  BUFFER  ADDRESS 

130 

AND  NIT 

/ 

131 

LDNF 

4 

132 

CLA 

1  ^ 

133 

TAD  ND5H 

/LOAD  HI  HORD  COUNT 

1  ^ 

131 

AND  K17 

>35 

LHNC 

136 

CLA 

137 

TAD  ND5L 

/LOAD  LO  HORD  COUNT 

133 

LDNC 

i 

139 

CLA 

4l 

140 

TAD  HD3L 

/FORN  CONNAND  HORD 

5  - 

141 

AND  K17 

1  i 

142 

SSN 

1 

4 

143 

CLL  RTL 

i 

144 

TADZ  TNOP 

145 

LDCN 

/LOAD  CONNAND  HORD  AND  GO 

146 

DCAZ  TCON 

/SAVE  CONNAND  HORD 

i 

147 

SNNB 

148 

JNP  .-1 

i 

>49 

CIF  CDF  0 

:i 

150 

ION 

/TURN  THE  INTERRUPTS  BACH  ON 

i£F£C  <■000  FORTRAN  lU  i-MOD£  SUBROUTIHtS 


Cl 

JMPZ  PDPT 

'  C2 

/ 

'S3 

/ -3- MODE 

PARAMETER  i 

’S4 

/ 

55 

PDPI, 

0 

'56 

CLA 

'S7 

TAD  ND2L 

!5e 

AND  K3 

'5? 

TAD  TN0P*1 

>6,0 

DC  A  TNOP 

>61 

TAD  ND3L 

>62 

AND  K3 

'.63 

DCA  ND3H 

>64 

TAD  ND4L 

>65 

AND  K3 

!66 

DCA  ND4H 

>6T 

TAD  NDSL 

>60 

AND  K3 

!6? 

DCA  MD5H 

>70 

TAD  NDSL 

!71 

CLL  RTL 

J72 

TAD  ND4H 

173 

RTL 

>74 

TAD  ND5H 

>75 

RTL 

>,76 

TAD  ND2L 

>.77 

DC At  TNOP 

\70 

CIF  CDF  0 

>.7P 

JHPt  PDPI 

'.SO 

/ 

>31 

/ 3- MODE 

RESISTER  S. 

yF!2 

/ 

>33 

PDPMl . 

0 

>34 

lOF 

>  V  5 

SKNB 

':36 

3MP  .-1 

'  3  , 

RDNC 

233 

DCA  NDSL 

>36 

RHNC 

>60 

AND  K17 

261 

DCA  HD5H 

262 

TADt  TCOM 

>63 

DCA  NDSL 

264 

DCA  ND3H 

265 

RDST 

296 

DCA  ND2L 

267 

RDES 

263 

DCA  ND2H 

CIF  CDF  0 

300 

ION 

:  VONE 

'TIRLIIRTION  (MODE  ^  0* 

/CALCULATE  TAPE  NOP  ADDRESS 

/MASK  FORMAT  MORD 

/MASK  PARlTy-DENSITV  NORD 

/MASK  THRESHOLD-SAP  MORD 

/BUILD  NOP  COMMAND 


/■STORE  NOP  COMMAND 
/DONE 

riON  (MODE  =■'/•> 

/TURN  OFF  interrupts 

/MAIT  until  controller  NOT  BU 

/SET  LO  RESIDUAL  NORD  COUNT 

/SET  HI  RESIDUAL  NORD  COUNT 

/SET  last  unit  tape  COMMAND 

/SET  STATUS  RESISTER 
/SET  error  resister 


■'TURN  OH  interrupts 


■HESEC  000 

FORTRRN  10  ■ 

J'O/ 

JMFZ  POPMl 

j02 

/ 

lOS 

/PRRRNE 

TERS 

so  4 

/ 

305 

MODE , 

0 

306 

ONE , 

1 

30T 

7  NO, 

2 

SOS 

THRE. 

3 

309 

FOUR, 

4 

310 

FIVE, 

5 

31 1 

HD2, 

27 

312 

ND2N, 

0 

313 

ND2L  . 

0 

314 

ND3. 

l‘7 

315 

ND3H, 

0 

316 

ND3L  , 

0 

317 

NB4, 

27 

31S 

ND4N, 

0 

319 

ND4L  , 

0 

320 

ND5, 

27 

321 

NOSH, 

0 

322 

ND5L  , 

0 

323 

JSUF , 

27 

524 

JSUFH, 

0 

325 

JBUFL , 

0 

326 

THOP, 

RDDR  TNOPO 

327 

TNOPO, 

240 

320 

THOPl, 

241 

329 

TH0P2 , 

242 

330 

TN0P3, 

243 

331 

TCOM, 

RDBR  TCOMO 

332 

TCOMO . 

0 

333 

TCOHl , 

0 

33  4 

TC0M2, 

0 

i3  S 

TCOM 3, 

0 

j  Jo 

K3, 

3 

J  J7* 

K17 . 

17 

339 

K7000, 

7000 

339 

M6000, 

-6000 

C£  SLIBROUriNES. 

.  oONk 

/ILINT/IUHT/lREe 

/ icon/ IFHT/TCOH 

/ISUF/PRDN/JBUf 

/ICNT/RTF6/JCNT 

/HOP  JNSTRUCTION  POINTER 
/244  FRS7  SOP 
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